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Rotary-Vane Attenuator With an Optical Readout 


W. E. Little*, W. Larson*, and B. J. Kinder* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colorado 80302 


(May 13, 1970) 


A precision x-band rotary-vane attenuator has been developed that follows its theoretical cosine 
squared law of attenuation to within 0.002 dB to a vaiue of 20 dB. The best available attenuation 
measurement techniques verify this relationship to within the accuracy of the measurements. From 
20 to 50 dB, the deviation from the theoretical value slowly increases and is frequency dependent. 
At present, three attenuators of this type have been constructed by NBS. Two of these attenuators 
have excellent attenuation versus frequency characteristics in the 20 to 50 dB range (0.08 dB devia- 
tion at 50 dB) while the third has a marked degree of frequency sensitivity (0.21 dB at 50 dB). 


Key words: Attenuator; microwave; precision; rotary vane. 


1. Introduction 


The principle of the rotary-vane attenuator [1,2]! has 
been known for a number of years, and relatively good 
attenuators have been commercially available for some 
time. A considerable amount of analysis [3,4,5,6,7] of 
this type of commercially available attenuator has been 
made, and some work has been directed toward the devel- 
opment of more precise attenuators [8,9]. Most of the 
major problems have been identified and analyzed. 

The steps taken in the work described in this paper 
were to eliminate or reduce the major limiting errors. The 
results achieved using this approach were very gratifying. 
The evaluation of the three constructed attenuators has 
shown that the property of attenuation difference up to 20 
dB is improved an order of magnitude over that reported 
previously. However, two of the attenuators retain their 
excellence above 20 dB, while the third shows some fre- 
quency dependence at higher values of attenuation. A 
description of the attenuators and the results of the meas- 
urement on two of them by the i-f substitution [10], d—c 
substitution [12], and modulated subcarrier [13] tech- 
niques of attenuation measurement are given. 

The “true” rotary vane value is estimated by averaging 
the results from the three measurement techniques. Meas- 
urement results from the third attenuator by the sub- 
carrier technique are also given. 


2. Theory 


As is well known, the ideal behavior of the rotary-vane 
attenuator is given by 
* National Bureau of Standards, Electromagnetics Division, Boulder, Colo. 80302 


1 Figures in brackets indicate the literature references at the end of this 
paper. 
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where E£;, indicates the strength of the input wave and 
E out is the strength of the output wave. 
In decibels, this is 


A = —40 logio (cos y). 


Thus, ideally, the output wave is a simple function of 
the angle of rotation of the center section of the atten- 
uator. 

There are many problems encountered in the actual 
construction of this type attenuator which has many 
sources of error. The errors fall roughly into five cate- 
gories: (1) read-out resolution and resettability; (2) 
basic mechanical precision such as vane alinement, ro- 
tating section eccentricity, resolution of vane angle set- 
ting, etc.; (3) external leakage, primarily from the rota- 
ting joints; (4) internal leakage (lack of absorption of 
the parallel component); and (5) internal reflections. 
Mariner [8] gives adequate analysis of categories 2, 4, 
and 5. James [9] also considers these errors as well as 
category 1 for the particular type tape readout which he 
was using. Larson [7] further considers the vane aline- 
ment error. The errors caused by category 3 have not 
been considered and in fact have not been readily observ- 
able in the attenuators that have been tested to date, 
primarily because there was not sufficient resolution to see 
the results of the leakage or the problem was not identi- 
fied. In the development of the attenuators described here, 
this problem was identified and had to be corrected for 
before substantial improvement could be accomplished. 








3. Description of the Attenuator 


3.1 Readout 


One of the major limitations of the design of this type 
of attenuator has been the technique used to determine 
the actual vane angle setting. In the commercial low- 
resolution (and some recent high-resolution) instruments, 
the technique is to simply use a gear reduction system 
with all the associated mechanical problems. These errors 
have been analyzed by Larson [4,6] and it is quite con- 
clusive that the technique rapidly becomes the limiting 
factor in the design of an ultra-precision attenuator. James 
[9] used a tape-dial indicator arrangement but, again, it 
is not suitable for ultra-precision attenuators. (See fig. 1.) 

To achieve the desired accuracy, the readout has to be 
a direct-reading device attached directly to the rotating 
cylinder. Fortunately, there are commercially available 
engraved glass scales with an optical readout that have the 
desired accuracy and resolution. Figure 2 shows a sketch 
of one of the scales mounted on the cylinder of the atten- 
uator’s rotating center section. The scale has an absolute 
accuracy of + 1” of arc and a resolution of the same 
amount; hence, the total uncertainty is + 2” of arc. Fig- 
ure 3 shows the attenuation deviation in dB from nominal 
for the variation of + 2” of arc versus attenuator dial 
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Ficure 1. Basic components of rotary-vane attenuator. 
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Ficure 2. Engraved glass scale mounted on rotating cylinder. 


setting in decibels. If this were the only source of uncer- 
tainty, the values of the uncertainty in attenuation would 
be about + 0.00025 dB and 0.003 dB at the dial setting 
of 10 and 50 dB respectively. 

A sketch of the dual microscope readout is illustrated 
in figure 4. Note the two viewing apertures are located 
180 degrees apart in order to reduce the eccentricity 
error. 


3.2 Mechanical Precision 


The final quality of the attenuator, even with the en- 
graved scales, is primarily a function of the precision and 
accuracy with which it is assembled. Larson [7] shows 
that the vane alinement is quite critical and a relatively 
small misalinement can cause a considerable error. Fur- 
ther, in order to utilize effectively the + 1” of arc resolu- 
tion of the engraved scale, rotational precision is a must. 

The required precision was obtained by carefully ma- 
chining and strengthening the attenuator housing as 
shown in figure 5. In some of the later testing of atten- 
uator No. 1, it was observed that even with the improved 
housing there was a nonrepeatability of setting of ap- 
proximately 10 »B. The cause of this problem was found 
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Ficure 3. Attenuator resolution. 
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Ficure 4. Dual microscope mounted in attenuator housing. 





Ficure 5. Modified attenuator housing. 


to be the forces exerted on the rotating cylinder by the 
drive gear even though the finest quality triple A bearings 
had been used in the rotating joints. The problem was 
eliminated by driving the cylinder through a bellows as 
shown in figure 6. This addition was not made to atten- 
uators No. 2 and 3. 


3.3 External Leakage 


The effects of external leakage are quite subtle and in 
many cases are very difficult to identify. In the early de- 
sign of the attenuators, this leakage caused nonrepeata- 
bility problems at the lower values of attenuation and in- 
stability at the higher values. The primary source of the 
leakage signal was determined to be the rotating joints. 
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Ficure 6. Bellows drive of attenuator No. 1. 


The addition of chokes alone did not provide sufficient 
attenuation to stop the leakage. In some cases, the leaked 
power is down only 40 dB from the incident power in the 
waveguide and it is rarely down more than 60 dB. In 
the attenuators described here, the leakage was reduced to 
more than 120 dB below the incident power by liberal 
use of microwave-absorbing material in the design of the 
housing and rotating joints as shown in figure 7. 


3.4 Internal Leakage 


Internal leakage is the result of insufficient attenuation 
of the parallel component of the wave (the component 
that is in the plane of the resistive card). The card ab- 
sorption was checked by measuring the maximum atten- 
uation that could be achieve with a particular card. For 
commercially available cards, this ranged from 110 dB 
to 135 dB at 9 GHz. The results of this leakage component 
is a deviation from the cos*y law at the higher values 
of attenuation. Both Mariner [8] and James [9] analyzed 
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Ficure 7. Addition of microwave absorbing material 
around rotating joints. 6 








this error. The actual effect of this leakage for the atten- 
uators is shown in figures 8, 9, and 10. Some work was 
done on the design of different types of cards but no im- 
provement was achieved, 
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Ficure 9. Frequency sensitivity of attenuator No. 2. 
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Ficure 10. Frequency sensitivity of attenuator No, 3. 





4. Results 


In order to test the attenuators over the desired 50 dB 
range, several techniques of attenuation measurement 
were used. The d—c substitution [12] system is the most 
accurate technique of attenuation measurement available 
in the range of 0.01 to approximately 20 dB. The modu- 
lated subcarrier [13,14] technique is quite accurate in the 
range of 0.1 to 20 dB while the i-f [15] substitution tech- 
nique is quite good in the range of 1 to 50 dB. All the 
techniques are comparable over parts of the overlapping 
range up to 20 dB. Above 20 dB, the d—c substitution 
technique loses resolution and the particular subcarrier 
system used for the measurements had not been carefully 


evaluated above 20 dB at the time of the measurements. 
Attenuators No. 2 and No. 3 were calibrated by all 
three techniques of attenuation measurement at several 
values between 0.5 and 20 dB at 9 GHz as shown in table 
1. (Values above 20 dB were obtained by the i-f tech- 


TaBLE 1. Comparison of calibration data for the same attenuator 
by three different methods 
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Average 
Nominal measured limits to Standard 
value, value, systematic deviation, 
dB dB error, 
Attenuator No. 2. 

0.5000 d—c 00.5006 0.0006 .0001 
i-f 00.4990 .0150 .1000 

sc 00.4996 .0050 .0005 

1.0000 d—c 01.0018 .0013 .0006 
i-f 01.0002 .0150 .0023 

sc 00.9997 .005 .0016 

3.0000 d—c 03.0016 .0047 .0001 
i-f 03.0000 .0150 .0023 

sc 02.9993 .0050 .0014 

6.0000 d—c 06.0021 .0049 .0001 
i-f 05.9986 .0150 .0028 

sc 05.9984 .0050 .0025 

10.0000 d—c 10.0020 .0053 .0055 
i-f 09.9990 .0150 .0013 

sc 09.9973 .0050 .0021 

20.0000 d—c 20.0021 .0100 0042 
i-f 19.9994 .0150 .0013 

sc 19.9946 .0100 .0023 

Attenuator No. 3. 

0.5000 d—c 0.5002 0.0006 0.0004 
i-f .4988 .0150 .0010 

sc .5007 .005 .0002 

1.0000 d—c 1.0003 .0013 .0002 
i-f 0.9992 .0150 .0010 

sc 1.0003 .0050 0011 

3.0000 d—c 3.0007 0047 .0002 
i-f 3.0010 .0150 .0010 

sc 3.0003 .0050 0015 

6.0000 d—c 6.0008 .0049 .0005 
i-f 6.0000 .0150 .0023 

sc 6.0007 .0050 .0027 

10.000 d—c 9.9996 .0053 .0016 
i-f 9.9992 .0150 .0010 

sc 9.9994 .0050 .0025 

20.0000 d—c 20.0069 0100 .0014 
i-f 19.9964 0150 .0013 

sc 20.0000 0100 .0045 
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nique only and are shown in table 3). The estimated sys- 
tematic error for each system is the summation of the 
estimated limits to all of the known system errors (this 
gives the estimated maximum possible systematic error) 
and the standard deviation is determined from the meas- 
urement data. 

The estimated true values for the attenuators given in 
table 2 are the average of the values obtained from the 
three different systems. 


TABLE 2. Estimated true value for each nominal value determined 
from the average of the three different methods 














Estimated 
Nominal true 
value, value, 
dB dB 
Attenuator No. 2 
0.5000 00.4997 
1.0000 01.0006 
3.0000 03.0003 
6.0000 05.9997 
10.0000 09.9996 
20.0000 19.9987 
Attenuator No. 3. 
0.5000 00.4999 
1.0000 00.9992 
3.0000 03.0005 
6.0000 06.0004 
10.0000 10.0003 
20.0000 20.0018 





TaBLe 3. Measured values—i—f substitution technique 
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Nominal Soom limits to Standard 
value, a systematic deviation, 
dB 4B’ error, dB 
dB 
Attenuator No. 2 
30.000 30.001 +0.018 0.001 
40.000 40.000 + .028 .003 
50.000 49.968 + .085 .004 
Attenuator No. 3 
30.000 29.995 +0.018 0.003 
40.000 39.982 + .028 .003 
50.000 49.945 + .085 .003 





5. Conclusion 


The development work outlined in this report confirms 
that over a limited range a precision-built, rotary-vane 
attenuator has excellent properties for a primary reference 
standard of attenuation difference. These properties are 
nearly ideal up to the attenuation value of 20 dB. Above 
this value of attenuation the rotary-vane attenuator is an 
excellent working standard. (The attenuator needs cali- 
brating because of the frequency-sensitivity problem.) 
Solution of the remaining two outstanding problems 
(frequency sensitivity and internal leakage) will estab- 


lish it as a primary reference standard over a much larger 
attenuation range. 

The attenuator should find wide application for the 
evaluation of systems where it is not feasible to use sec- 
tion-by-section evaluation techniques. For example, in 
some of the dual channel i-f attenuation-measuring sys- 
tems, the standard piston atienuator cannot be indepen- 
dently evaluated. In addition to being used in the straight 
rf substitution technique, this rotary-vane attenuator can 
be used in the same fashion that the 30 MHz waveguide- 
below-cutoff atenuator is presently used. For the higher 
frequencies (millimeter waves), a frequency down con- 
verter can be used to obtain the proper frequency and for 
the lower frequencies, an upconverter can be employed. 

Another quite useful area for this rotary-vane attenua- 
tor is in the measurement of very small values of attenua- 
tion. To the authors’ knowledge, there is no other way to 
achieve a change of attenuation of a few hundredths or 
tenths of the decibel to the accuracy and resolution that 
can be achieved with the ultra-precision rotary-vane 
attenuator. 
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Design and Operational Characteristics of a High-Speed 
(Millisecond) System for the Measurement of Thermophysical 
Properties at High Temperatures: 


Ared Cezairliyan 
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Design and constructional details of a high-speed system for the measurement of selected 
thermophysical properties of electrical conductors at temperatures above 2000 K in experiments of 
subsecond duration are described. Operational characteristics of such a system are given. Various 
phenomena that affect the design and successful operation of the system are discussed and, when- 
ever possible, quantitative results are given. Certain experimental checks to assess the operation of 


the system are also described. 


Key words: High-speed measurements; high temperatures; thermodynamics; thermophysical prop- 


erties. 


1. Introduction 


Interest in thermophysical properties at temperatures 
above the limit of accurate steady-state experiments has 
necessitated the development of high-speed measurement 
techniques. Most of the difficulties inherent in steady- 
state methods are eliminated in high-speed experiments. 
However, new problems arise due to the nature of rapid 
experiments, some of which increase rapidly with in- 
creasing speed. This suggests that an optimized system 
has to be designed for most accurate measurements. 

The objective of this writing is to describe the details 
of a high-speed system for the measurement of selected 
thermophysical and related properties (heat capacity, 
electrical resistivity, hemispherical total and normal spec- 
tral emittances, melting point, etc.) of electrical conduc- 
tors at temperatures above 2000 K in experiments of sub- 
second duration. Operational characteristics of such a sys- 
tem under actual experimental conditions are presented. 
Various phenomena that affect the design and successful 
operation of the system are discussed in the context of 
high-speed experiments and, whenever possible, quantita- 
tive results are given. Certain experimentai checks to 
assess the operation of the system are described. 

The method employed in the high-speed measurements 
is based upon rapid resistive heating of the specimen by 
a single pulse (subsecond duration) of direct current and 
measuring the pertinent quantities with millisecond reso- 
lution. Formulation of the relationships between various 
experimental quantities and properties is given in another 
publication [1]? 


* This work was supported in part by the Propulsion Division of the U. S. 
Air Force Office of Scientific Research under Contract ISSA-70-0002. 

1 Figures in brackets indicate the literature references at the end of this 
paper. 


2. Measurement System 


The high-speed measurement system consists of an elec- 
tric power pulsing circuit and associated measuring and 
control circuits. A functional diagram of the complete 
system and a photograph of a portion of the system are 
presented in figures 1 and 2, respectively. 


2.1. Pulsing Circuit 


The pulse circuit includes the specimen in series with 
a battery bank, a standard resistance, a variable resist- 
ance, and a switching system. The battery bank consists 
of 14 series-connected 2 V batteries each having ap- 
proximately 1100 A-—h capacity. A standard resistance 
(0.001 © shunt) made of manganin strip is used to meas- 
ure the pulse current flowing through the specimen. A 
variable resistance made of water-cooled Inconel tubes 
(total resistance — 30 mQ) enables control of the heat- 
ing rate of the specimen and the shape of the current 
pulse. The switching system consists of two series-con- 
nected, fast-acting switches. The second switch is used as 
a backup in the event the first one fails to open at the 
end of the heating period. The switches are operated by a 
series of pulse control units. 

The pulsing circuit is a simple RL circuit whose dy- 
namic characteristics, after closing of the switch, may be 
expressed by 


E—L% =i (R, +R) (1) 
dt 
where 
R, = resistance of the specimen 
R, = resistance of the circuit with the exception of the 
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Ficure 1. Functional diagram of the complete high-speed measurement system. 


L = self-inductance of the complete circuit 
E = battery potential 
i >= instantaneous current 
¢ == thee: 
The solution of eq (1) is 
: a 5 oy — gee 
t= (l—e ) (2) 


where R = resistance of the complete circuit (R =R, + 
R.). Equation (2), which gives the increase of current as 
a function of time, is obtained on the assumption that the 
quantities E, L, and R are constant. In an actual experi- 
ment, this assumption is satisfied approximately for E 
and L, however, R increases with time as the specimen is 
heated. Therefore, the current versus time relationship is 
of a distorted exponential form. The circuit’s “time con- 
stant,” which is expressed by the quantity L/R, decreases 
with increasing R as the specimen is heated. 

Energy stored in the magnetic field of an inductor, with 
no ferromagnetic material in its vicinity, is given by 


Li? 
y -—. 3 
5 (3) 


It may be noted that stored energy depends on the instan- 
taneous value of the current and not on the current his- 
tory. In a typical high-speed measurement system, in 
which inductance of the specimen is of the order of 10-7 
H, and current is the order of 10° A, magnetic energy is 
of the order of 0.1 J. This is small compared to imparted 
energy, which is the order of 10* J. For a complete cycle 
of operation, consisting of closing the switch and after a 





Ficure 2. The experiment chamber and the high-speed pyrometer. 


finite period opening the switch, the stored energy in the 
specimen is zero. 


2.2. Measuring and Control Circuits 


Voltage signals from both the standard resistance and 
the specimen are sent to the recording system via differ- 
ential amplifiers. A potentiometric system is used to 
calibrate pulse-voltage and pulse-current measuring cir- 
cuits including the differential amplifiers and the digital 











recording system. A Kelvin bridge is used to measure 
resistance of the specimen at ambient temperature before 
and after a pulse experiment. A 400 Hz synchronous 
motor used in the high-speed pyrometer provides the time 
base for the data acquisition system. Timing of various 
events, such as closing of the switch and triggering of 
various electronic equipment, during a pulse experiment 
is achieved by a series of time-delay units. Triggering of 
the recording system and opening of the switch at the end 
of the desired heating period are automatically activated 
when the specimen temperature reaches preadjusted initial 
and final values. 


2.3. Specimen 


The specimen is a tube of the following nominal dimen- 
sions: length = 4 in (102 mm), outside diameter = 0.25 
in (6.3 mm), wall thickness = 0.02 in (0.5 mm). A small 
rectangular hole (1 0.5 mm) is fabricated in the wall 
at the middle of the specimen to approximate blackbody 
conditions. In order to compensate for the cross-sectional 
nonuniformity created by the hole, a portion from the rest 
of the specimen is removed by grinding its surface. The 
criterion is that the cross section of the removed portion 
is equal to that of the opening created by the hole. From 
geometrical considerations this may be expressed as 


R? [{ 76 
Di ae [eee 955 
a 5 (4 sin ) (4) 
where 
a = width of hole 
b = depth of hole (thickness of tubular specimen) 
R = outer radius of specimen 
‘= 


angle subtended by the are corresponding to the 
ground flat. 


After solving the above transcendental equation for the 
angle @ by iteration, the pertinent dimensions as shown 
in figure 3 may be obtained from the following relations: 


L=2R dine 


; (5) 
C = R cos = (6) 








Ficure 3. Specimen cross section at the plane of the flat. 


A photograph of the specimen is shown in figure 4. 

The rectangular sighting hole in the specimen wall is 
fabricated 0.8 mm offcenter to improve the blackbody 
quality. For the specimen geometry given in an earlier 
section, blackbody quality of the hole is estimated using 
De Vos’ [10] method. Computations are made assuming 
four different specimen inner surface reflectances and for 
two specimen lengths, one corresponding to the “effec- 
tive” length and the other to the actual length. The results 
(fig. 5) do not show any appreciable dependence of the 
blackbody quality on length and reflectance. 


2.4. Experiment Chamber 


The experiment chamber, shown schematically in fig- 
ure 6, contains the specimen, the clamping electrodes, 
an expansion joint, potential probes, and other auxiliary 
components. 


Ficure 4. Photograph of the specimen. 
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Ficure 5. Blackbody quality of a specimen computed for 
different inner wall reflectances and lengths. 
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Ficure 6. Experiment chamber. 


The specimen is mounted vertically 6 mm offcenter 
with respect to the axis of the chamber to reduce the 
effect of internal reflections. The inner wall of the 
chamber is coated with nonreflecting paint. The chamber 
wall, as well as the specimen clamps, are water cooled. 
Thermocouples are connected (electrically insulated) to 
the two end clamps to measure the specimen temperature 
before each pulse experiment. An expansion joint allows 





the expansion of the specimen in the downward direction. 
The effect of specimen evaporation on the window trans- 
mission is minimized by placing a flat partial shield in 
front of the specimen. 

The potential probes are knife-edges made of the speci- 
men material and are placed at a distance approximately 
13 mm from the end clamps. The knife-edges define an 
“effective” portion of the specimen, which should be free 
of axial temperature gradients for the duration of the 
experiment. A photograph of the specimen with the clamps 
and potential probes is shown in figure 7. 

The chamber is designed for conducting experiments 
with the specimen either in vacuum or in a controlled 
atmosphere. 


2.5. High-Speed Pyrometer 


Temperature of the specimen is measured with a high- 
speed photoelectric pyrometer, which permits 1200 evalu- 
ations of the specimen temperature per second. The 
pyrometer alternately passes precisely timed samples of 
radiance from the specimen and a reference source 
through an interference filter (wavelength 650 nm, band- 
width 10 nm) to a photomultiplier. During each exposure 
the photomultiplier output is integrated and is recorded. 
Successive exposures to the reference source, which is a 
calibrated, gas-filled tungsten-filament lamp, are taken 
through a sequence of three different optical attenuators 
mounted on a rotating disk, resulting in a staircase of 
reference exposures with approximately 50 percent at- 





Ficure 7. Photograph of the specimen, 
clamps, and potential. probes. 
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tenuation for each step. This scheme is used for measure- 
ments up to 2500 K, the limit of reliable operation of 
gas-filled lamps. For measurements at higher temperatures 


calibrated optical attenuators are placed in the path of 


the radiation from the specimen. The details regarding 
the construction and operation of the pyrometer are given 
in the literature [2]. 


2.6. High-Speed Digital Data Acquisition System 


The high-speed digital data acquisition system consists 
of a multiplexer, analog-to-digital converter, and a core 
memory together with control and interfacing equipment. 
All signals are brought to the multiplexer through differ- 
ential amplifiers in order to avoid inaccuracies arising 
from common ground points. The multiplexed signals go 
to the analog-to-digital converter, which has a full scale 
reading of + 10 V and a full-scale resolution of one part 
in 8192 (8192 = 2!*). Digital output from the converter 
consists of 13 binary bits plus a sign bit. This output is 
stored in a core memory having a capacity of 2048 words 
of sixteen bits each. The data acquisition system is capable 
of recording a set of signals corresponding to temperature, 
voltage, and current approximately every 0.4 ms. At the 
end of the pulse experiment, information stored in the 
memory is retrieved in the form of numeric printing and 
punched tape using a teletypewriter. Since the laboratory 
has access to a time-shared computer, it is possible to 
unload the memory directly to the computer, bypassing 
the intermediate stage of punching paper tape. A manu- 
script regarding the constructional and operational details 
of the high-speed digital data acquisition system is in 
preparation [3]. Oscilloscopes are used only to monitor 
the general pattern of the experimental results, and to 
detect any anomalies. 


3. Operational Characteristics 


Highest measurement precision is achieved when the 
experimental quantities are recorded close to the full- 
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Ficure 8. Oscilloscope trace photographs of 
typical voltage and current pulses. 


Equivalence of each major division is: time = 50 ms, voltage = 2 V, 
current = 1000 A. 





scale value of the data acquisition system. However, since 
quantities vary with time it is not possible to realize this 
completely. The recording of current can be optimized by 
adjusting the series resistance, which minimizes this vari- 
ation. Also, the proper adjustment of the series resistance 
assures the quasi-linear variation of both voltage and 
current pulses. Oscilloscope trace photographs of typical 
voltage and current pulses are shown in figure 8. Radiance 
from the blackbody hole in the specimen, as well as that 
of the reference lamp, as seen by the pyrometer as a func- 
tion of time is shown in figure 9. Dots forming the hori- 
zontal lines correspond to the three reference radiances. 
Temperature of the specimen is computed from the inter- 
polation of its radiance between adjacent reference 
radiances. 

In the computation of properties third-degree poly- 
nomial functions for voltage, current, and temperature as 
functions of time are obtained by least squares approxi- 
mation of the data. Typical standard deviation of an 
individual point from the pertinent function is 0.02 per- 
cent for voltage, 0.03 percent for current, and 0.5 K for 
temperature. 

Typical results on imparted power to the specimen, the 
heating rate of the specimen and the ratio of heating to 
cooling rate (M) for a molybdenum specimen computed 
from measurements of voltage, current, and temperature 
are shown in figures 10, 11, and 12, respectively. They 
represent a total of 40 experiments that were performed 
at two different speeds, designated as “fast” and “slow,” 
over the temperature range 1900 to 2800 K. The smooth 
curves shown in the figures are obtained by fitting the 
computed points to third-degree polynomial functions 
using least squares approximation. The standard deviation 
of an individual point from the pertinent function is 0.5 
percent for power, 0.7 percent for heating rate, and 0.8 
percent for the ratio of heating to cooling rate. The small 
values of the standard deviations are an indication of the 
high degree of reproducibility of the current pulses and 
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Ficure 9. Oscilloscope trace photograph of 
radiance of rapidly heating specimen. 
Dots forming the Yong horizontal lines correspond to radiances from the 
reference source. Equivalence of each major division is: time = 20 ms, 
radiance = arbitrary unit. 
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Ficure 10. Variation of the power imparted to a molybdenum 
specimen as a function of temperature. 





o “N we 
T | i T | 
y 
QP 
yn 
4 
eee eee ee dieereneen 





Ta 
E 
=x 
w = 

a 
ae 
9 

254+ : 

4 | 

¥ | 

. 4 

on: SLOW \ 

| 

i 7 








1800 2000 2200 2400 


TEMPERATURE , K 


2600 2800 


Ficure 11. Variation of the heating rate of a molybdenum 
specimen as a function of temperature. 


the entire measurement system. Reproducibility of the 
pulses was also determined by performing a total of 12 
experiments at 5 min intervals. The maximum difference 
in the measured current values was less than 0.02 percent. 

Although heating rate of the specimen depends, among 
other factors, on the power supply and the external circuit 
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Ficure 12. Variation of the quantity M (ratio of heating to cooling 
rate) of a molybdenum specimen as a function of temperature. 


characteristics, the cooling rate following the opening of 
the switch depends on the specimen material and its 
geometry. 

Power balance during cooling yields the following 
expression for the cooling rate of the specimen: 


dT a . 
(4 ), a 


B= eA / Cyl 

= effective surface area of specimen 
heat capacity 

weight of effective specimen 
specimen temperature 
hemispherical total emittance 
Stephan Boltzmann constant 
(5.6697 « 10° W m= K~*) 


cooling rate. 


(7) 
(8) 


where 


qa salem 


(dT /dt),= 


The cooling rate as a function of temperature for a 
typical molybdenum specimen is shown in figure 13. 

In order to obtain an expression for the specimen tem- 
perature during cooling as a function of time, eq (7) is 
rearranged and integrated. The resultant relation is 


1 1/3 
3Br+_1. 
ea 


r= 
(9) 
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Ficure 13. Variation of the cooling rate of a molybdenum 
specimen as a function of temperature. 
where 
t — time 
T= specimen temperature at the start of cooling 
(att —0). 


Variation of the temperature of a molybdenum speci- 
men during cooling as a fungtion of time corresponding 
to different initial temperatures is shown in figure 14. 

In preliminary calculations and also in computations 
related to error analysis it is necessary to know the quan- 
tity P,/P;, where P, is power loss from the specimen due 
to thermal radiation and P; is power input to the speci- 
men. From power balances for heating and cooling peri- 
ods the following relation is obtained: 
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Ficure 14. Variation of the temperature of a molybdenum speci- 
men during cooling from various initial temperatures. 


where 
(dT /dt) }, 
(dT /dt),.° 


It may be seen that for the steady-state case (d7/dt), = 
0, which implies that M = 0, and thus P, = P;. 


M= (11) 


4. Consideration of Various Phenomena 


In this section various physical and electrical phenom- 
ena that affect the design and operation of high-speed 
thermodynamic measurement systems are presented. 


4.1. Axial Temperature Distribution 


Heat transfer from specimen to clamps causes the 
establishment of temperature gradients in the specimen. 
However, since the voltage probes are 12.7 mm away 
from the clamps, sharp temperature gradients near the 
clamps do not affect the measurements. The magnitude 
of the axial temperature gradient in the “effective” speci- 
men is estimated by solving the transient heat conduction 





P, 1 equation assuming constant properties and zero radial 
P,- 14+M (10) heat transfer [4]. The series solution is given by the 
; following equation: 
co 
_ ( Ppl? a By (—1)" [| 4s olan Ur] 12 
= ($4) i (2n + 1)8 “°° 2 —- 40 (3) 
n=0 
where A= thermal conductivity 
i= electric current a= thermal diffusivity. 
p= electrical resistivity 
A = cross-sectional area of specimen If radial heat transfer is considered, the actual axial tem- 
l= half length of “effective” specimen perature gradient will be somewhat smaller than that 
x= axial distance from midpoint of specimen given by eq (12). 
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The temperature gradient between the midpoint of the 
molybdenum specimen and various points along the speci- 
men is calculated using eq (12) (with n= 30) for os 
various heating pulses. The normalized results are shown 
in figure 15. Computations for heating pulses of 350 and - 
700 ms in duration indicate that specimen temperature 
at the plane of the potential probes is approximately 99.9 
and 99 percent of its midpoint value, respectively. The 
respective average temperature of the “effective” specimen i 
is 99.99 and 99.9 percent of its midpoint value. 

To ascertain the temperature uniformity, a photograph ¢— TusuaR _, 
of a molybdenum specimen was taken during the last part SPECIMEN 
of the heating period. This photograph is shown in figure 
16 in conjunction with the schematic of the specimen. 
No variations in exposure are detected for the portion 
between the potential probes. The region of sharp tem- 
perature gradients was studied photomicrographically. A 
photomicrograph of a molybdenum specimen, obtained 
after heating it to 2800 K for several times with pulses 
ranging from 0.4 to 0.8 s is shown in figure 17. Pertinent 
locations on the photomicrograph are identified with re- 
spect to the specimen. Two distinctly different grain 
structures may be seen. The portions of the specimen that 
were in the clamp or close to it have retained the original 
“drawn” pattern, while the portion that was heated to a 
high temperatures shows a considerable grain growth. 
It is interesting to note a clear boundary between the 
two patterns, which correspond to a plane (P1) approxi- 
mately 3 mm from the plane (P2) of the end clamp. > 
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4.2. Radial Temperature Distribution cusmP 








There is no exact solution for the time-dependent radial 
temperature distribution in conductors carrying pulse 
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Ficure 16. Photograph of a molybdenum specimen 
during pulse heating to 2500 K. 


0.8}- currents with radiative heat transfer from their surfaces. 


7 
| 
e 
| Even the numerical methods become very complicated 
7 because of the nonlinearity of the differential equation. 
In pulse experiments where heat loss from the surface 
06 4 is small compared to power imparted to the specimen, the 
| radial temperature drop (temperature difference between 
4 inner and outer surfaces for tubular conductors) may be 
| estimated using the following approximate relation: 
al 
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Koh 
Toh 





04 AT, =N aft, (13) 
eer J 
= TIME 

= HALF LENGTH 


= THERMAL DIFFUSIVITY 


where 


AT, = temperature drop under transient conditions 
| AT, = temperature drop under steady-state conditions 
| N = ratio of heat loss from the surface to imparted 
: power under pulse conditions. 
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The solution of the steady-state radial heat conduction 
; - ah ae T3 "7 equation, assuming constant properties, for radial tem- 





x perature drop is: 
roeoT* 9 9 9 r 
Ficure 15. Axial temperature distribution in a resistively AT, = ase — 17°) — 2r? In (=)| (14) 
self-heated specimen under transient conditions. A(To — ) nm 
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Ficure 17. Photomicrograph of a molybdenum specimen 
after pulse heating to 2800 K for several times. 
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i-* (2 + 26, + 6m) —e 


where 
T) = outer radius 
r; = inner radius 
o = Stephan-Boltzmann constant 
« = hemispherical total emittance 
d = thermal conductivity. 


Computations on a molybdenum specimen indicate 
that radial temperature drop at 2800 K may be approxi- 
mately 0.3 K for a heating pulse of 0.4 s. 


4.3 Evaporation 


The rate of evaporation m(T) from metallic surfaces 
may be expressed by 


m(T) = Ae ®/t (15) 


where, 7 is temperature in K, A is a constant for a given 
surface, and E is equal to the energy of evaporation in 
the units of the Boltzmann constant. 

Since in high-speed experiments the heating period is 
short compared to the cooling period, one may assume 
that most of the evaporation takes place during cooling 
of the specimen. Total evaporation during cooling becomes 


to lz 
u= |" m(T)ae= | Ae *®/Tdt. 
th ty 


After performing a time-to-temperature transformation 
using eq (9) and integrating eq (16), one obtains 


(16) 


—# is 
(2+ 26,4 e)| 





where 
E 
On = = 
Tn 
E 
6, = — 


B = quantity expressed by eq (8) 


and 7, and 7; are initial and final temperatures, respec- 
tively. 

Computations on a molybelenum specimen indicate that 
the total weight loss due to evaporation during cooling 
from 2800 K is approximately 0.003 percent of the speci- 
men weight. 


4.4. Inductive Effects 


Using equations given in the literature [5], self- 
inductance of the standard resistance and the specimen 
were computed to be 0.13 and 0.018 uH, respectively. 

The maximum rate of change of the electric current 
during measurements in a typical high-speed experiment 
is approximately 3000 A s"‘. Then, the possible error 
(due to self-inductance) in potential measurements across 
the standard resistance and the specimen is 0.4 mV and 


0.054 mV, respectively. Both of these values are smaller 
than the resolution of the recording system (1 mV). 

In order to check the magnitude of mutual-induction 
between the main circuit and the signal circuits, the fol- 
lowing experiment was performed. Potential connections 
to the standard resistance and the specimen were dis- 
connected and were shorted duplicating, as much as 
possible, the geometrical configuration that the leads have 
during an actual high-speed experiment. A heavy pulse 
current was allowed to flow through the main circuit and 
the signals were recorded. The output of the recording 
system did not show any indication of voltage signals 
(within 1 mV resolution) that could have been generated 
as the result of mutual inductive coupling between the 
main circuit and the signal circuits. 


4.5. Skin Effect 


The skin effect, which is caused by changing currents 
in a circuit, alters the temperature distribution in the 
specimen. 

For cylindrical conductors the skin effect is a function 
of the following parameter: 


x = 2rav/ (2pf/10°%p) (18) 





where 


a = radius of the conductor 

f = frequency of current 

p = electrical resistivity of the specimen 
# = permeability of the specimen. 


For nonmagnetic materials, » = 1. 
In the case of tubular (thin-wall) specimens, the skin 
effect is a function of the following parameter: 


B= 2 r/aV/2 (19) 


where 


x = quantity defined by eq (18) 
a = outer radius of the specimen 
7 = thickness of the specimen. 


In general, the contribution of the skin effect in a 
specimen is expressed in terms of the ratio of its electrical 
resistance at the frequency in question to that under d-c 
conditions. Results on the resistance ratio in terms of the 
parameter x and £ are given in tabular form in the 
literature [5]. Computations indicate that for the speci- 
men geometry and the type of pulses used in the present 
measurement system, the contribution of the skin effect 
is less than 0.01 percent. 


4.6. Magnetic Forces 


Magnetic force (per unit length) on a tubular specimen 
resulting from current flowing through it is expressed as 


r=($) [HE 
—\3b/ La+m)? 


where i =current, m — a/b, b = outside radius, a = 
inside radius. The corresponding magnetic pressure is 
2i? 


given as 
bl ) | m | 
” 3arb* (1+ m)? 


Excessive magnetic forces, which act towards the center 
of the specimen, tend to collapse it. The collapsing pres- 
sure for thin-wall tubes is 

ya —m)* 


(20) 


(21) 





2E 


p.=(2, 


=, (22) 


where E = modulus of elasticity and » = Poisson’s ratio. 
The maximum allowable current is obtained by equat- 

ing eqs (21) and (22), which yields 
‘| (23) 
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1+ 2m 
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4.7. Thermionic Emission 
The contribution of thermionic emission is manifested 
in the following two forms: (1) energy loss from the 
specimen as a result of energy carried away by emitted 
electrons, (2) establishment of undesirable current paths 
around the specimen. 
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Dependence of thermionic emission current on tem- 
perature is given by 


eo. 


I=AT*e ™ (24) 


where 


A = thermionic emission constant 
¢ = thermionic work function 

e = electronic charge 

k = the Boltzmann constant. 


Considering the values A = 55 & 10* A m~*K~ and 
¢@ = 4.15 eV [6,7] for molybdenum, thermionic emission 
current is computed to be approximately 7 & 1075 
A mm ~ at 2000 K and 0.15 A mm~ at 2800 K. 

It may be noted that computed values are for the case 
where electrons are continuously removed from the sur- 
face, which requires a collector plate at positive potential. 
In the present system, the chamber wall is electrically 
neutral and the only electrical potential gradient that 
exists is that across the specimen, which is less than 10 V. 
Thus, in the present experiments thermionic emission is 
most probably space charge limited; its contribution 
therefore would be much less than above computed values. 


4.8. Thermoelectric Effects 


Thermoelectric effects play an important role in experi- 
ments where temperature gradients exist and where cur- 
rents flow through interfaces of dissimilar conductors. 
The Peltier effect between the specimen and clamps and 
the Thomson effect in the specimen where sharp tem- 
perature gradients exist are appreciable. However, they 
do not affect the measurements in the present experiments 
since their contribution is negligible in the “effective” 
specimen. To further reduce the possible contribution of 
thermoelectric effects, the knife-edges used for potential 
measurements are made of the same material as that of 
the specimen. In the high-speed experiments, the record 
of potential readings at high temperatures after the open- 
ing of the main switch shows “zero” for the potential 
between the probes. This indicates that there are no 
measurable erroneous signals due to the interaction (in 
thermoelectric sense) of the probes with the specimen. 


4.9. Thermodynamic Equilibrium 


Whether a specimen is under thermodynamic equilib- 
rium while measurements are taken during rapid heating 
and cooling periods may be determined by considering 
the relaxation times of various rate processes. 

Total relaxation time that must be considered for 
thermodynamic equilibrium calculations may be expressed 
as 


1 1 

where each 7; represents a different mechanism of re- 
laxation. The various relaxation mechanisms are the 
results of different scattering processes. The main three 
relaxation processes are due to electron-phonon, electron- 


imperfection, and phonon-phonon interactions. 











In electrical conductors, electron-phonon interaction is 
the predominant one in determining the overall relaxation 
time. For estimative purposes, the relaxation time due to 
this may be expressed as 


om 


ne* 


where e = electronic charge, m = electronic mass, n = 
number of electrons per unit volume, and o = electrical 
conductivity. 

Relaxation time due to electron-phonon interactions for 
metals at room temperature is of the order of 1071s. 
Since relaxation time is directly proportional to electrical 
conductivity, and since the latter decreases with increasing 
temperature, it takes a longer time to establish thermal 
equilibrium at high temperatures. The variation of elec- 
trical conductivity in metallic elements between room 
temperature and their melting points is within a factor of 
100; therefore, it can be seen that thermal equilibrium 
in metals takes place in times of the order of 107s. 

Information concerning relaxation times for vacancy 
equilibrium is not readily obtainable; but since the 
mechanism is related to that of diffusion, the relaxation 
time is an exponential function of temperature. The mag- 
nitude of vacancy relaxation time as observed in quench- 
ing experiments at lower temperatures and extrapolated 
to higher temperatures may be of the order of a milli- 
second for refractory metals above 2000 K. 


Figs 


(26) 


5. Experimental Checks 


The following sections outline the procedures and the 
results of selected experiments, which were performed to 
assess the accuracy of measured quantities under various 
operational conditions. 


5.1. Resistance Comparison 


The objective was to determine the value of a test 
resistance by measuring the current flowing through it 
and the voltage drop across its terminals under transient 
conditions and comparing this measurement with its re- 
sistance under steady-state conditions as determined by 
the potentiometric method. 

In order to eliminate any doubt of the possibility of 
the effect of induced voltages on electrical signals, sepa- 
rate experiments were conducted using two geometrically 
different test resistances. One of the resistances (0.001 0) 
was a flat manganian strip identical to that used for pulse 
current measurements. The other was an Inconel tube 
with a nominal resistance of 0.003 0. The results on both 
test resistances do not indicate any dependence attribut- 
able to geometry. In these experiments pulse currents in 
the range 800 to 2100 A were used. The results do not 
indicate any current dependence. The results on resistance 
comparison and standard deviation of the measurements 
are shown in figures 18 and 19, respectively. It may be 
seen that, in general, resistance measurements by pulse 
and steady-state techniques are in agreement within 0.03 
percent. 
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Ficure 18. Resistance difference of a test resistor measured 
under steady-state and pulse conditions. 
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Ficure 19. Standard deviation of resistance measurements 
under pulse conditions. 


5.2. Effect of Different Heating Rates 


Several phenomena depend upon the heating rate of 
the specimen. Heating rate is related, among other factors, 
to the electric current that flows through the specimen. 
In high-speed experiments, the rate of change of the cur- 
rent creates a number of problems such as induced volt- 
ages in measuring circuits, varying forces in the circuit 
components, etc. In addition to these extrinsic effects, 
there is the consideration of the specimen’s attaining 
thermodynamic equilibrium at the time of measurements. 
To demonstrate partially that the above are negligible, 
in the present system under present operating conditions, 
the following experiments were conducted. 








A tungsten strip (3.2 mm wide, 0.25 mm thick, 100 mm 
long) was placed in the high-speed measurement system. 
Two voltage probes were welded 20 mm apart at the 
center portion of the strip. Two pulse experiments were 
conducted in which the strip was heated from room tem- 
perature to 2200 K (surface brightness temperature) in 
approximately 2000 and 250 ms, respectively. During the 
experiments, quantities required for the determination 
of voltage, current, and surface brightness temperatures 
were recorded. Comparison of the resistances obtained 
from the two (slow and fast) experiments corresponding 


to the same brightness temperature yields the following 
average value: 


R(slow heating) / R(fast heating) — 1.0003. 


The resultant 0.03 percent deviation is within the uncer- 
tainty of resistance measurements. Thus, it may be con- 
cluded that approximately eightfold variation in the heat- 
ing rate does not introduce any significant differences in 
the measured properties. 


5.3. Melting Point of Platinum 


As a partial check on the accuracy of temperature 
measurements, experiments were conducted to obtain the 
melting point of platinum. A high-purity (99.999 percent) 
platinum wire (diameter — 0.001 in, 0.025 mm) was 
wound (half turn) around a thin-wall molybdenum tube 
(similar to regular specimens) at the plane of the black- 
body sighting hole. Two pulleys supported the wire. A 
0.5 g weight was attached to each end of the wire to 
assure good contact between the wire and the tube. The 
ends of the platinum wire were connected to a special 
circuit, which indicated the time of the opening of the 
circuit resulting from the melting of platinum. The break 
in the circuit took place at a point of contact with the 
tube as the latter was pulse heated. Duration of an experi- 
ment (heating of the specimen from room temperature 
to its melting point) was approximately 0.8 s. Tempera- 
ture of the tube corresponding to the time of the break 
in the platinum wire gave the melting point of platinum. 

The results of three experiments gave an average value 
of 2044.2 K for the melting point of platinum on the 
1968 International Practical Temperature Scale [8] with 
an average deviation from the mean of the three amount- 
ing to 1.1 K. This compares favorably with the accepted 
value of 2045 K [8]. Details regarding the experiments, 
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interpretation of results, and errors are given in another 
publication [9]. 


6. Conclusions 


Estimates indicate that the imprecision and inaccuracy 
of measured quantities in high-speed experiments using 
the present system are: 0.02 and 0.05 percent for voltage, 
0.03 and 0.06 percent for current, and 0.5 K and 4 K for 
temperature at 2000 K. 

The high-speed system described in this paper has been 
used successfully to measure heat capacity, electrical re- 
sistivity, hemispherical total and normal spectral emit- 
tances of molybdenum [1] and tantalum [11] from 1900 
to near their respective melting points. The system was 
also used to obtain the melting point of molybdenum [12] 
and the variation of normal spectral emittance of tantalum 
during melting [13]. The results indicate that at 2000 K 
accuracy of measured properties by the system described 
here is about the same as that obtained by conventional 
methods. At temperatures above 2500 K, it surpasses any 
other method now known. Estimates of the imprecision 
and inaccuracy of measured properties are given in the 
above cited references. 
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Chromatographic paper, saturated with ammonia, is used to extract the residual thiosulfate from 
the gelatin layer of processed film. The paper is treated in a silver nitrate solution and fixed in an 
ammonia-sodium chloride solution. Thiosulfate, if present, forms silver sulfide and the thiosulfate 
is determined by measuring the transmission density of the paper darkened by silver sulfide. The 
paper extraction method is simple, rapid, requires only a densitometer as special equipment, and 
is very sensitive especially when the transmission density is measured for two layers of the paper. 
No filters are required. The test readily reveals the uneven distribution of residual thiosulfate on 
the film. By selective removal of thiosulfate from the image silver by ammonia and potassium 
bromide solutions, it was shown that a small amount of thiosulfate was absorbed on the image silver 


immediately after processing. 


Key words: Archival film; hypo; microfilm; permanent record film. 


1. Introduction 


Residual thiosulfate is readily determined in processed 
photographic papers by measuring the transmission den- 
sity of the silver sulfide formed in situ by the reaction of 
the silver ion with the residual thiosulfate [1].2 However, 
the transmission density of the silver sulfide in the gelatin 
layer of film is much lower for corresponding thiosulfate 
concentrations, making it difficult to use the density of 
the silver sulfide as a measure of residual thiosulfate con- 
centration in film [2, 3]. 

It has been found that chromatographic paper wetted 
with an ammonia solution will absorb the residual thio- 
sulfate from the gelatin layer of the film. Such paper was 
treated in silver nitrate solution forming silver sulfide in 
situ and the excess silver ion was removed by treating the 
paper in an ammonia solution. The transmission density 
of the silver sulfide was large enough to determine 
thiosulfate concentrations as low as 0.15 y»g/cm*. The 
film was not destroyed in the test. 


2. Experimental Procedure 


Whatman No. 3 MM chromatographic paper* was 
selected for this study. The paper, which was 0.03 cm 
thick, withstood the treatments in the solutions and drying 
without distortion. A 2-percent ammonia solution was 
found to be very efficient in extracting the residual thio- 
sulfate from the gelatin. The ammonia solution should 
either be used in a hood or covered during the test. 

A strip of film 10 to 12 cm long and 1.6 to 2 cm wide 
was taped, emulsion side up, at the two ends to blotter 





1 Retired. Present address: 2911 Dawson Ave., Silver Spring, Md. 20902. 

2 Figures in brackets indicate the literature references at the end of this paper. 

3 Commercial materials and equipment are identified in this paper in order to 
specify the experimental procedure adequately. Such identification does not imply 
recommendation nor endorsement by the National Bureau of Standards. 
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paper resting on a level glass plate. The chromatographic 
paper, cut the width of the film sample and of a length 
to fit between the tape strips at the ends of the film, was 
soaked 3 min in about 100 ml of a 2-percent ammonia 
solution. Then the paper was grasped with tweezers at one 
end, drained briefly, and the other end brought in contact 
with the film, allowing the paper sample to roll down onto 
the film to prevent air bells. The paper was rubbed once 
gently with the back of the tweezers. The blotter paper 
absorbed any excess ammonia solution at the edges of the 
film. After 5 min the paper was removed and bathed 
4 min in 100 ml of a solution containing 10 g of silver 
nitrate and 30 ml of glacial acetic acid per liter. A control 
strip was made by treating a sample of the chromato- 
graphic paper in the 2-percent ammonia solution 3 min 
before the test sample was removed from the film and 
the two paper samples were processed together. Next, the 
papers were placed in crystallizing dishes and treated 3 
min in each of two solutions containing about 100 ml of 
5-percent ammonia solution and 1 g of sodium chloride 
per liter, to remove the excess silver ion. The solutions 
were stirred during treatments. Then the papers were 
removed from the ammonia solution with tweezers and 
placed on blotting paper but not blotted. The paper was 
turned over three or four times during a 10-min period 
and allowed to dry in the room atmosphere. The trans- 
mission density was measured for the control and test 
papers. The effective optical density of the silver sulfide 
was determined by subtracting the density of the control 
paper from that of the test paper. The corresponding 
residual thiosulfate concentration of the film was deter- 
mined from a calibration curve. 

It was found that 3 to 5 percent ammonia solutions 
containing 1 to 5 g of sodium chloride per liter may be 








used to fix the paper samples. The sodium chloride stabi- 
lized any trace of ammonia-silver complex left in the paper. 

The chromatographic paper was sufficiently uniform in 
density that only an occasional control run would be 
necessary in routine testing. 

A pink dye was absorbed by the paper from some films 
during the extraction operation. The dye was bleached 
out by exposing the dried paper sample for 1 h to a 20-W, 
20-in daylight fluorescent lamp at a distance of 1 in. 

The transmission density measurements reported in this 
paper were made with an Ansco-Sweet Densitometer, 
Model 11. In this paper, the concentration of residual 
thiosulfate is given as micrograms per square centimeter. 
In most of the literature, the concentration is reported 
in milligrams or micrograms per square inch. One micro- 
gram per square inch is equivalent to 0.155 yg/cm*. 


3. Evaluation of the Paper Extraction Method 


In evaluating the paper extraction method for the de- 
termination of residual thiosulfate in processed film, six 
microfilms, two motion picture films, two x-ray films, one 
aerial film, one print film, and one professional film were 
tested. The films were made by five manufacturers. The 
unexposed film samples were treated 4 min in a solution 
containing 90 g sodium sulfite, 52 g sodium carbonate 
monohydrate, and 5 g of potassium bromide per liter. 
This solution is typical of developers but without a de- 
veloping agent. The developing agent was eliminated to 
avoid any reduction of silver in the gelatin layer. Then 
the film samples were fixed in an acid hardening bath 
and washed in chlorine-free tap water as described in a 
previous publication [4]. The excess water was removed 
from the film samples with wiping tissue. The washing 
time controlled the residual thiosulfate concentration in 
the films. 

For comparison, the residual thiosulfate concentrations 
in the processed film samples were determined by the 
paper extraction method and the modified Crabtree-Ross 
method [5]. Samples for the two methods were taken 
adjacent to each other on the film. The thiosulfate concen- 
trations ranged from 0.15 to 16 ug/cm* of sodium thio- 
sulfate. The two methods were in general agreement. 
However, the paper extraction method revealed the fact 
that the residual thiosulfate is not always evenly dis- 
tributed in the gelatin layer, there being higher and lower 
values for residual thiosulfate concentration than the 
average value found by liquid extraction using the modi- 
fied Crabtree-Ross method. 

Tests were made to determine if any thiosulfate re- 
mained in the gelatin layer of the film after the first 
extraction. A second extraction was made with the paper 
extraction method immediately after the first extraction. 
The effective optical density of the silver sulfide was 
determined for two layers of the test paper and the resi- 
dual thiosulfate concentration was read from the calibra- 
tion curve shown in the next section. The results are 
tabulated in table 1. For 11 films, no thiosulfate remained 
where there had been a low original thiosulfate concentra- 
tion and only a negligible amount remained where it had 
been high. Two films with a thick layer of gelatin retained 
about 5 percent in the low thiosulfate concentrations and 
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Taste 1. Residual thiosulfate found in a 
second ammonia extraction of processed film. 





Approximate micrograms of residual 
thiosulfate per centimeter originally 
found in the film after processing 


15 3 8 16 





Micrograms of thiosulfate per 
square centimeter found in the 
second ammonia extraction 








Six microfilms 

Two x-ray films 

One motion picture film 0 0 
One aerial film 

One print film 


0 to 0.15 0 to 0.7 


One motion picture film ” _ 
One professional film } 0 to 0.07 0.15 0.7 15 








10 percent in the high. A third extraction detected no 
thiosulfate. The presence of thiosulfate in the gelatin 
layer after the first, second, or third extractions was also 
determined by processing the film, after extraction, b 
the procedure used by Crabtree, Eaton, and Muehler a 
for determining the thiosulfate in photographic papers. 
Using this method, one can detect 0.08 pg/cm? of thio- 
sulfate as silver sulfide by measuring the density of two 
layers of film with a Kodak 18A filter. No thiosulfate was 
detected in the film with low original concentrations and 
the amount found where concentrations had been higher 
correlated with that found in the second paper extraction. 

The effect of a high concentration of silver in the fixing 
bath on the paper extraction of thiosulfate was investi- 
gated. Film was fixed in a fixing bath containing 4 g of 
silver per liter, the silver being added as silver bromide. 
The silver thiosulfate complex left in the film did not 
affect the extraction of the thiosulfate. 

Tests were made to determine if any residual develop- 
ing agents remained in processed film as these would 
form reduced silver in the test paper in the silver nitrate 
solution. Unexposed samples of microfilm, motion picture 
film, and x-ray film were normally processed and washed 
free of thiosulfate. The samples were treated in a bath 
containing 10 g sodium sulfite, 5 g boric acid, 10 g 
potassium aluminum sulfate, 10 mi glacial acetic acid, 
0.05 g metol, and 1 g hydroquinone per liter and washed 
30 s in running tap water. The above bath contained 
chemical compounds found in a used fixing bath in addi- 
tion to thiosulfate. When the paper extraction test was 
applied, there was no reduction of silver in the paper. The 
developing agents were rapidly removed from the films 
and therefore would not interfere with the paper extrac- 
tion method. 

Paper extraction does not remove all of the thiosulfate 
from areas on the film containing a large amount of 
image silver. One microfilm was flashed and developed 
to a density of 4.8 and washed 30 min. The paper extrac- 
tion method detected no thiosulfate in either the image- 
free area or the image area, whereas the modified 
Crabtree-Ross method gave negative results for the clear 
area but 1.2 »g/cm* of thiosulfate in the image area. A 
sample of the high-density area was treated 15 min in a 
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5-percent potassium bromide solution, washed, and 
bleached with an untreated sample in a dichromate solu- 
tion. The effective densities of the remaining silver sulfide 
measured for two layers with a Kodak 94. filter were 0.78 
for the untreated sample and 0.72 for the one treated in 
potassium bromide solution. The silver sulfide density of 
0.72 corresponds to the silver sulfide formed in the image 
during fixation and the difference of 0.06 density cor- 
responds to thiosulfate absorbed on the image silver. 
These tests, made immediately after the film was washed, 
showed the absorbed thiosulfate was not removed by 
ammonia but was by potassium bromide. The amount of 
thiosulfate absorbed is very small. The above results 
indicate that the paper extraction method should be ap- 
plied to an image-free area of the film or one not having 
more than 0.1 density of silver. The absorbed thiosulfate 
is unstable and reacts rapidly at high relative humidity 
with the silver in the image to form silver sulfide. 


4. Calibration Curve 


A calibration curve was prepared using the Whatman 
3 MM chromatographic paper. A sample of the paper, 
cut 10 by 20 cm, was treated 5 min, with agitation, in 
two liters of a sodium thiosulfate solution in a photo- 
graphic tray. Solutions containing 0.050, 0.065, 0.075, 
0.10, 0.15, 0.2, 0.3, 0.4, and 0.8 g of Na2S203°5H20 per 
liter were used. The paper sample was removed from the 
thiosulfate solution and blotted immediately between clean 
photographic blotters with pressure, blotted a second 
time between clean blotters, and allowed to dry flat on 
clean blotters, turning the sample over every 3 or 4 min 
for the first 30 min to obtain a more even distribution 
of the sodium thiosulfate in the paper. 

The dried paper sample was cut into two parts, 5 by 
10 cm, and one part was treated, together with a control 
sample of the chromatographic paper containing no thio- 
sulfate, 4 min in a 400-ml solution containing 10 g of 
silver nitrate and 30 ml of glacial acetic acid; fixed 3 min 
in each of two 400-ml solutions containing 3-percent 
ammonia and 5 g of sodium chloride per liter; and dried 
flat on blotter paper. The transmission densities of the 
paper sample containing silver sulfide and the paper 
control were measured. The density of the silver sulfide 
was obtained by subtracting the density of the control 
from that of the paper sample containing the silver sulfide. 
A single layer of the chromatographic paper had a density 
of 1.19 + 0.01, two layers had 1.47 + 0.01 with an 
occasional deviation of 0.02. Treatment in the processing 
solutions did not change the density. 

Four 214-cm squares were cut from the second part of 
the sample mentioned above, adjacent to the first part 
and about 5 cm from the ends, and the residual thiosulfate 
was determined by the modified Crabtree-Ross test 
method [5]. In this test method 5 ml of a 5-percent 
solution of mercuric chloride was added, with rapid stir- 
ring, to a vial containing 5 ml of a 5-percent solution of 
potassium bromide and 100 ug of gelatin. The thiosulfate 
reacted with the mercuric chloride to form a turbidity 
which was determined with a nephelometer. The thiosul- 
fate concentration was determined from a calibration 
curve. Four squares were used in the test for sodium 
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thiosulfate concentrations up to 2.5, 2 squares for 5, and 
1 square for 10 »g/cm*. The 5-percent solutions of potas- 
sium bromide and mercuric chloride were added from 
burets. The samples of paper were treated 5 min in a 
50-ml- beaker with 5 ml of the potassium bromide solu- 
tion, while agitating with a thin glass rod, and were 
transferred to a 10-ml graduate. The samples were treated 
a second time with 5 ml of the potassium bromide solution 
to remove any trace of thiosulfate remaining in the paper 
and the solution transferred to the same graduate. After 
the paper samples were removed from the beaker the 
volume of solution in the 10-ml graduate was made up to 
10 ml with the 5-percent potassium bromide solution to 
replace the solution absorbed by the paper and the 10-ml 
solution was transferred to the beaker in which the extrac- 
tion was made. Four drops of a 0.1 percent gelatin solu- 
tion and then 10 ml of the mercuric chloride solution 
were added to the beaker with rapid stirring. Ten ml of 
the solution in the beaker were transferred to a vial and 
allowed to stand 30 min and the turbidity was determined 
with a Coleman No. 9 nephelometer. The vial was handled 
gently to avoid stirring up a few fibers of the paper that 
had settled to the bottom. The concentration of the sodium 
thiosulfate was determined from a calibrated curve for 
known concentrations of sodium thiosulfate in the potas- 
sium bromide—mercuric chloride solution. This value was 
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Ficure 1. Silver sulfide calibration curve showing the transmission 
density of silver sulfide formed in the Whatman No. 3 MM 
chromatographic paper by the reaction of sodium thiosulfate 
with silver nitrate. 








multiplied by 2 and divided by the area in square centi- 
meters tested to obtain the concentration of the sodium 
thiosulfate in micrograms per square centimeter in the 
paper. The transmission densities of the silver sulfide for 
one layer and two layers of the paper were plotted against 
the corresponding sodium thiosulfate concentrations to 
obtain the calibration curves shown in figure 1 for the 
Whatman No. 3 MM chromatographic paper only. Cali- 
bration curves must be prepared for each type of paper 
used in extracting the thiosulfate from the gelatin layer 
of film. The calibration procedure is very simple if a 
nephelometer is available for measuring turbidity. 


5. Summary 


The new paper extraction method is simple, rapid, and 
requires no light filters. Only a densitometer is needed as 
special equipment. The method is nondestructive as tests 
may be made on roll or sheet film which, after drying, 
may be returned to storage. The residual thiosulfate may 
be extracted in the field and the paper may be forwarded 
for final processing to a laboratory which has a calibra- 
tion curve for the type of paper used in the extraction. 
The test should be made in an image-free area of the 
film or one not having a density greater than 0.1. The 
test is quite sensitive, especially if the silver sulfide of two 
layers of the paper is measured. The test demonstrated 





that the residual thiosulfate is not always evenly distrib- 
uted in the gelatin layer. Residual sodium thiosulfate as 
low as 0.15 ug/cm* can be detected if the density of the 
silver sulfide for two layers of paper is measured. This 
concentration is visually discernible for one layer viewed 
with a control. For this reason the paper extraction test 
is an excellent qualitative test for the detection of residual 
thiosulfate in processed film. The test method will deter- 
mine trithionate and tetrathionate as thiosulfate. This is 
advantageous since these compounds also react with the 
image silver to form silver sulfide. 


6. References 


{1] Crabtree, J. I., Eaton, G. T., and Muehler, L. E., A review of 
hypo testing methods, J. Soc. Motion Picture Engrs. 42, 
34-57 (1944). 

[2] Mattey, D. A., and Henn, R. W., Determination of thiosulfate 
and thionates in film with silver nitrate, Photo. Science and 
Engr. 10, No. 4, 202-208 (1966). 

[3] Pope, C. L, A simplified method for determining residual 
thiosulfate in processed microfilm, Photo. Science and 
Engr. 13, No. 5, 278-279 (1969). 

[4] Pope, C. I, Blemish formation in processed microfilm, J. 
Res. Nat. Bur. Stand. (U.S.), 72A (Phys. and Chem.), 
No. 3, 251-259 (1968). 

[5] Pope, C. L, Determination of residual thiosulfate in pro- 
cessed film, J. Res. Nat. Bur. Stand. (U.S.), 67C (Eng. 
and Instr.), No. 3, 237-245 (1963). 


(Paper 75C1-309) 


22 


—a2 








JOURNAL OF RESEARCH of the National Bureau of Standards — C. Engineering and Instrumentation 
Vol. 75C, No. 1, January-March 1971 


Scanning Densitometer for Continuous Recording of 
Spectral Transmission Density at Low Spatial Contrast 


W. L. McLaughlin, M. Rosenstein,’ E. K. Hussman, and J. J. Lantz, Jr. 


Center for Radiation Research, National Bureau of Standards, Washington, D.C. 20234 
(October 6, 1970) 


A simple photoelectric transmission densitometer can be converted in the laboratory into an 
instrument for measuring spatial variations in spectral density of small areas of thin, nonturbid 
films. This is accomplished by adding relatively inexpensive components that are readily available 
and easily assembled. The basic instrument can be one of a number of commercial densitometers, 
consisting of stabilized electronics, an adjustable light source, and a photosensitive probe which 
has a wide range and measures spectral transmission density. For continuous scanning of relatively 
small sample areas, the followng additions are needed: (1) narrow-bandpass filters; (2) inter- 
changeable field illuminating slits; (3) a short-focal-length illuminating-condenser lens which com- 
promises between a relatively large numerical aperture and a reasonable depth of focus; (4) a 
motor-driven stage with adjustable speeds; (5) a light diffuser and light guide; and (6) a con- 
tinuous, relatively fast-response data-recording system (such as a fast-response strip-chart recorder) . 
This instrument has been assembled and tested especially for measuring spatial variations of 
spectral density produced in irradiated photochromic films; and also for calibrating a train of 
irradiated film samples, in terms of transmission density at a given wavelength as a function of 
irradiance, radiation absorbed dose, or photochemical reactivity. 


Key words: Densitometer; microdensitometer; radiochromic dye films; spectral transmission density. 





1. Introduction 


With the increasing use of thin gel, polymer, and 
photochromic films in chromatography and ionizing radi- 
ation applications, such as obtaining beam profiles, reg- 
istering chemical-reactivity, depth-dose, or penetration 
distributions, monitoring high-intensity radiation fields, 
and monitoring beam uniformity in radiation processing 
of industrial products, there is a need for continuous 
optical-density scanning at given wavelengths across rela- 
tively small areas of these films. A commercial micro- 
densitometer is suitable for very high resolution work, 
such as determining granularity of images, modulation 
transfer functions of imaging systems, or acutance of 
photographic images [1—4],? but is expensive and difficult 
to use for large volumes of routine work. The aim of the 
present work is to show how a simple, inexpensive photo- 
metric apparatus can be assembled and used to measure 
spectral transmission densities [5] of relatively nonscat- 
tering thin-film images with moderate spatial resolution. 

Most microdensitometers have two objective lens-slit 
systems: an illuminating and a scanning system. The 
illuminating system is situated between the light source 
and the sample. The scanning objective lens focuses a 
certain area of the narrow sample image on a reading slit 
in front of the photocathode. These components are 


* This work was partly supported by the Division of Isotopes Development, U.S. 
AEC 


1 Bureau of Radiological Health, Environmental Health Service, Rockville, Md. 
and Graduate Student, Laboratory for Radiation and Polymer Science, Univer- 


sity of Maryland, College Park, Md. Guest worker at the National Bureau of 
Standards. 


2 Figures in brackets indicate the literature references at the end of this paper. 
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arranged so as to reduce errors due to flare from the 
sample, diffraction from the apertures caused by partial 
coherence in the illumination, and film granularity [6, 7]. 
Great pains have to be taken to improve performance by 
adjusting the components so as to minimize angular and 
focal misalinements and non-uniform photocathode sen- 
sitivities [8]. The slit images falling on the photosensor 
are narrow in order to reduce adjacency effects [9]. 
The instrument described here is not a conventional 
microdensitometer, in that a scanning lens-aperture sys- 
tem has not been used; rather the slit-image of narrow- 
band light focused on a small area of the relatively non- 
scattering photochromic film is transmitted according to 
the spectral density of that area, and then is diffused and 
channeled efficiently to the photomultiplier. A schematic 
chart displaying the main components is shown in figure 
1. Relatively turbid samples, such as developed silver 
halide emulsions, would cause excessive light-scattering 
in such an arrangement, so that a scanning objective lens 
would be needed to focus the small portion of sample 
image through a scanning slit on the photocathode. The 
instrument diagrammed in figure 1 should be restricted 
to use with relatively nonturbid, nonreticulated plastic 
films, photochromic layers, and chromatographic gels. 
The light passes through a multiple filter system to a slit 
whose image is reduced and focused on the sample. 
Essentially all the light transmitted by the sample is first 
diffused and then collected by the photomultiplier. The 
photomultiplier current (proportional to spectral density) 
can be read relative to the lateral film position by an 
electrometer, an electrometer connected to a strip-chart 
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Ficure 1. Schematic chart of main components of densitometer 
scanning slit and light collection system. 

















recorder, or, in an expanded system, an analog-to-digital 
converter may be used and the resulting pulses fed into 
a tape-recorder readout. The light-collecting system is 
similar to microscope systems used in recording histo- 
photometers for integral density scanning of small bio- 
logical specimens or for microscopic spectrophotometry, 
in that essentially all transmitted light from a small sample 
area is collected and measured [10—12]. A removable thin 
focusing screen can be placed into the plane of the sample 
for sharp focusing of the slit image while viewing with 
an 8-power lens. 

The main limitation of the nonfocusing diffuse-light 
gathering system is that only moderate spatial resolution 
is possible, even with relatively nonscattering media. The 
main advantages are that problems of diffraction, critical 
slit alinement, focal adjustment, and nonuniformity of 
photocathode sensitivity are less severe than for a micro- 
densitometer. Even though granularity noise in turbid 
materials, such as developed silver halide emulsions, would 
be averaged somewhat by the present system, the scat- 
tering of light by the same would result in flare and 
limited definition of the measured image. 


2. Instrumental Details 
The basic instrument, adapted for the purposes de- 


scribed here, can be one of a number of photoelectric 
densitometers with a wide response range, adjustable 
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lamp brightness, and stabilized electronics.? The adapta- 
tion for continuous recording of spectral densities of 
small areas can be accomplished at relatively low cost 
without disturbing the primary function or opera*ion of 
the original instrument. If necessary, the instrument can 
easily be returned to its original configuration to measure 
macroscopic diffuse optical densities of large specimens. 

Figures 2 and 3 are photographs of details of the added 
components. Figure 2 shows the drive, gear, and counter 
combination which is powered by a variable d-c power 
supply. Inexpensive permanent-magnet reversible-type 
motors with integral reduction gearing drive the worm- 
gear screwjack. Limit switches are mounted to the body 
of the jack to prevent over-travel. At full voltage, the film 





Ficure 2. Close-up photograph of drive, gear and 
counter assembly and d-c power supply. 





Ficure 3. Close-up photograph of scanning densitometer 
as arranged during operation. 


3 For example, the apparatus used in this application was a Welsh Densichron 
Model 400 Photometer and Densichron 3853D Transmission light source. (This 
equipment is identified here in order to specify the instrumental details ade- 
quately. Such identification does not imply recommendation nor endorsement by 
the National Bureau of Standards.) 














holder directly connected to the jack moves laterally in 
either direction across the photometric aperture at 2 
mm/s. Slower speeds variable down to 8 pm/s are 
possible by lowering the voltage with a rheostat and by 
using reduction gears. In order to monitor the relative 
sample position, a calibrated revolution counter is used. 
These parts provide satisfactory reproducibility and con- 
trol of the position of samples, without requiring more 
expensive synchronous motors or stepper motors and 
gear trains. 

Figure 3 shows the scanning densitometer in its operat- 
ing mode. The sample holder is continuously movable 
from left to right or vice versa in a horizontal plane across 
the slit. A vernier scale allows positioning of the sample 
perpendicular to the scanning direction. The vertical 
position of this whole assembly is also adjustable for use 
with different illumination lenses which focus the illumi- 
nation slit images in a given horizontal plane. The sample 
holder has a variable opening and can hold films or 
plates as large as 20 cm long by 3.8 cm wide. 

For the purpose of carrying out many simple spatial 
density gradient readings, an inexpensive drive-strip chart 
recorder system is satisfactory. If a more sophisticated 
readout is required, particularly in the case of large- 
volume data acquisition, an analog-to-digital converter 
can be used. In the latter case, a ramp voltage can be 
supplied via a servo loop to drive the motor in stages and 
thus give a record of optical density as a function of 
position. 

Since one of the main purposes of the instrument is 
to make spectral-density measurements at optical wave- 
lengths near the absorption maxima of particular photo- 
chromic and plastic films, special consideration has been 
given to providing intense essentially monochromatic light 
sources corresponding to the wavelengths of the absorp- 
tion peaks of interest. Interference filters which have been 
used successfully in spectrophotometry [13] were fitted 
to the instrument for this purpose. In this case, the optical 
density measured is standard spectral density [5, 14, 15]. 
Table 1 lists the filters used, according to wavelength 
passed, their characteristics in terms of maximum percent 
transmittance, and spectral bandpass full width at half 
maximum. Also listed are the plastic or dye films for 
which the change in absorption characteristics due to ir- 
radiation is most pronounced at the indicated wavelength. 


TaBLE 1. Interference filter characteristics and corresponding 
radiochromic film type measured at a particular wavelength 








Transmittance Maximum 








Band width Radiochromic 
wavelength percent at half max. film type 
(nm) transmittance (nm) 
401 44.7 10 polyvinylchloride 
436 55.5 10 methyl yellow 
malachite green 
559 59.4 10 pararosaniline 
new fuchsin 
601 58.5 10 hexahydroxyethyl] violet 
crystal violet 
620 60.0 10 malachite green 
640 80.0 10 blue cellophane 


red Perspex 








Several interchangeable slits have been made by an 
electrical-discharge machining process [16] in aluminum 
and blackened to reduce reflection from the edges. A list 
of these slits and their effective image size at the sample 
plane, when a 4-power, 3.0 numerical aperture condenser 
lens is used as an illumination condenser system, is given 
in table 2, 


TABLE 2. Effective image sizes of slits 
used with scanning densitometer 








Illuminating slit 


Effective image 
(width X height, mm) 


(width X height, mm) 





0.068 x 7.0 0.018 x 1.84 
-133 X 1.0 035 X 0.26 
25 X 2.0 .066 X 0.53 
50 X 3.0 13 X 0.79 

10 X75 26 X 1.97 





The lens is an achromatic, spherical-aberration-corrected 
short focal-length (5.5 mm) lens with an aperture vari- 
able from {/1.8 to {/22. The lens is positioned at a dis- 
tance from the slit somewhat greater than twice the focal 
length and the distance between the lens and sample falls 
between one and two focal lengths. In this geometry, the 
dimensions of the slit image focused on the sample are 
between one-third and one-fourth of the slit dimensions. 
It should be observed from the combination of the 
smallest effective slit width (18 um) and the numerical 
aperture of the lens at full aperture that diffraction at the 
wavelengths of interest is not a handicap. A slit of this 
size is still large enough to produce adjacency effects 
when used for viewing samples with high spatial fre- 
quency; however, in most images encountered in the 
present work, with the exception of edge images, only 
relatively small variations in optical density occur over 
a distance the size of the slit width. 


3. Instrument Performance 


It was noted earlier that the scanning densitometer is 
not suitable for evaluation of turbid photographic films 
with high spatial contrast due to excessive light-scattering 
effects. An example of this limitation is demonstrated in 
figure 4, where a knife-edge image made on photographic 
film is shown: (1) as traced by a commercial micro- 
densitometer, and (2) superimposed, as traced by the 
scanning densitometer. Due to the procedure for calibrat- 
ing the instrument, data from the scanning densitometer 
were originally obtained in terms of transmittance. For 
purposes of comparison with data in figures 4 and 5 
from the microdensitometer, these data have been con- 
verted to optical density, where 


Optical Density = — log (Transmittance). 


Both traces were made with 601 + 10 nm (full width at 
half maximum) light provided by an interference filter; 
the effective slit of the microdensitometer was 16 um 
800 um and the effective slit image of the scanning densi- 
tometer was 18 um X 1,840 um. An objective measure 
of sharpness for sigmoid-shaped traces is acutance, which 
is calculated from the relationship [17]: 





S* (AD/Ax)? 
i=] 


Acutance = ey, = 4) 


where, 

(AD/Ax)i = density change per corresponding trace 
distance (microns) for each of n equally 
spaced trace distances between D, and Dp. 

D.4,D,; = optical density at points on knife-edge 
trace where AD/Ax is 0.005 »m™. 


Acutance for the microdensitometric trace of the photo- 
graphic knife edge, using n — 10, was a factor of about 
10 higher than for the scanning densitometer. In addition 
to loss of resolution in the trace of the scanning densitom- 
eter, a second degradation of the knife-edge trace is ob- 
served due to flare from the illumination of the knife edge. 
Flare [7, 17] arises from light reflected into the emulsion 
when a bright region surrounds the image and is within 
the field of the detector. This is the case when a photo- 
graphic knife edge is viewed by the scanning densitometer. 
The result is a reduction in optical density due to col- 
lection of extraneous light in the detector. Flare is pri- 
marily a function of the material being viewed; it is not 
a serious problem with the nonturbid films for which 
the scanning densitometer was designed. 

Calculations have been reported of modulation transfer 
functions (MTF’s) for microdensitometer systems from 
knife-edge traces based on Fourier-transform techniques 
and convolution methods [17-19]. These computational 
techniques, with the aid of high-speed computers, can 
generate separate MTF’s for the photographic system and 
the microdensitometer from knife-edge scans and provide 
other basic information required for the evaluation of 
image quality. Such an elaborate analysis was not deemed 
appropriate for the scanning densitometer because the 
instrument is not proposed for high-resolution work, but 
rather as a mechanical apparatus for the purpose of driv- 
ing a scanning slit across a nonturbid film having rela- 
tively low spatial contrast. 
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Ficure 4. Comparison of photographic knife edge traces for com- 
mercial microdensitometer (dotted line) and scanning densi- 
tometer (solid line). 


Vertical marks on the traces indicate the points between which acutance was 
determined (AD/Ax = 0.005 um~1). 


The most practical measure of the performance of the 
scanning densitometer is by comparison of optical density 
traces of nonturbid films with traces produced by a com- 
mercial microdensitometer. A typical comparison is shown 
in figure 5, where traces of an exposed radiochromic dye 
film, which has recorded the spatial-distribution of dose 
(proportional to optical density) with depth in a medium, 
are drawn on the same scale. The important feature of 
the comparison is the overall agreement in the optical- 
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Ficure 5. Comparison between spatial optical density tracings of an irradiated radiochromic-dye film 
for commercial microdensitometer and scanning densitometer. 
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Ficure 6. Scanning densitometer trace of a series of irradiated radiochromic dye films. 


(1), (2), and (3) are calibration films exposed to a uniform radiation field, and (4) is the depth-dose response of a 
film inserted longitudinally in a 4 mm thick phantom. The steep vertical lines represent the edges of the individual films. 


density trend between the two traces. As expected, the film 
edge is more clearly defined in the commercial micro- 
densitometer trace. The sharp spikes evidenced there are 
responses to small scratches and dust particles which can- 
not be resolved by the scanning densitometer. It is also 
observed that the problem of flare noticed with the pho- 
tographic knife edge is not as pronounced at the edges 
in the traces of these radiochromic dye films. Except for 
the improvement in determining the position of the edge 
of the film, no great advantage would be gained in using 
the more sophisticated microdensitometer, since the in- 
strument resolution required to register gradual spatial 
variations in the dye films is available from the scanning 
densitometer. As an example of a practical application of 
this instrument, figure 6 shows a trace of a series of 
irradiated radiochromic-dye films read in tandem on the 
densitometer. 


4. Summary 


A simple, inexpensive, scanning densitometer, with 
effective viewing areas as small as (0.035 < 0.26)mm* 
or (0.018 < 1.84) mm?, has been designed and assembled. 
It is based on existing components, namely a commercial 
photoelectric densitometer, an illuminating optical-slit 
system with various wavelengths of monochromatic light, 
a variable-speed sample driving mechanism and table, and 
a continuous fast-response readout system. Experiments 
have shown that images of relatively low spatial contrast 
on thin nonturbid photochromic films can be scanned 
competitively with more expensive and complicated micro- 
densitometers. Advantages are that alinement, diffraction, 
focal adjustment, and photocathode non-uniformity are 
not severe problems with the instrument described. Dis- 
advantages are that high spatial frequency and sharp 
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knife-edge images on turbid media cannot be handled 
adequately with this system. 
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A rotating attenuator for the generation of subsecond duration sawtooth shape radiance pulses 
is described. The attenuator disk is in the form of a cam. The geometry of the opening of a 
diaphragm used as the aperture stop determines the shape of the radiance pulses. Radiance is 
determined from the measurements with a high-speed photoelectric pyrometer. Recording of signals 
is made with a high-speed digital data acquisition system. The combined measuring and recording 
systems have a full-scale signal resolution of approximately one part in 8000 and a time resolution 
of 0.4 ms. Two different diaphragms are used in this study yielding radiance pulses (20 to 150 ms 
long) with linear and quadratic rise. The standard deviation of the experimental points from the 
pertinent functions describing the radiance variation is less than 0.5 percent. 


Key words: High-speed methods; light attenuator; photometry; radiance modulation; radiation. 


1. Introduction 


A considerable amount of work has been done to 
develop various techniques to modulate radiance beams 
in the visible as well as in the ultraviolet and infrared 
regions. Most of these techniques were confined to gen- 
erating either repetitive rectangular pulses, or short- 
duration fast-risetime single pulses. A review of various 
methods for modulating radiance pulses was given by 
Jones [1]. Optical modulation techniques in the infrared 
region were discussed by Hudson [2]. 

Increasing interest in transient techniques for the 
measurement of properties of matter at high temperatures 
and related studies has necessitated the development of 
programmable optical attenuators. Such attenuators are 
needed to simulate (using steady-state radiance sources) 
the rapid heating or cooling of the specimen and also to 
check dynamic characteristics of radiation detectors and 
of high-speed temperature measurement systems in general. 


2. Method and Apparatus 


In this study an optical system, including a rotating 
optical attenuator, is described that is capable of generat- 
ing repetitive sawtooth shape (either linear or quadratic) 
radiance pulses, ranging from 20 ms to several seconds 
in duration. The optical system (fig. 1) consists of a radi- 
ance source (tungsten filament lamp), a rotating disk 
(cam with helical circumference), a special diaphragm, 
a high-speed pyrometer, and other pertinent optical com- 
ponents. The lenses 01 and 02 shown in figure 1 are 63 
mm in diameter with an aperture of {/2.5. The circum- 
ference of the disk was machined in such a way that its 

* This work was supported in part by the Propulsion Division of the U.S. Air 


Force Office of Scientific Research under contract ISSA-69-0001. 
1 Figures in brackets indicate the literature references at the end of this paper. 
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radius r was proportional to the angle from a fixed 
reference point according to the following relation: 


(1) 


where ro is the maximum radius, a is the angle, and k is 
the proportionality constant. The maximum and minimum 
radii of the disk were 7 in (177.8 mm) and 6 in (152.4 
mm), respectively, yielding a maximum difference of 1 in 
(25.4 mm). The thickness of the disk was 0.12 in 
(3 mm). 

The aperture stop of the optical system was defined by 
the circumference of the disk and the diaphragm placed 
in the collimated portion of the beam. As the disk rotated, 
the effective cross-sectional area of the beam changed, 
which, in turn, generated sawtooth shape pulses. The 
shape of the diaphragm opening determined the shape 
of the radiance pulse. For example, a rectangular opening 
as seen in figure 2 yields a pulse with linear rise, while 
a triangular opening gives a pulse with quadratic rise. 
Rotation of the disk was synchronized with other measur- 
ing equipment with a signal obtained from a radiation 
detector. The radiation detector PM (fig. 1) generated a 
signal each time a 0.2 mm in diameter hole H (fig. 2) 
fabricated in the disk was in line with the detector and 
the light source L (fig. 1). Relative radiance of the source 
after attenuation was measured with a high-speed photo- 
electric pyrometer [3]. The electrical signals from the 
pyrometer were recorded with a high-speed digital data 
acquisition system [4]. The combined measuring and 
recording system has a full-scale signal resolution of 
approximately one part in 8000 and a time resolution 
of 0.4 ms. 

The photomultiplier in the pyrometer was alternately 
exposed to the radiance from the tungsten filament lamp 
and to the radiance from a reference lamp. This scheme 


r=ro — ka 
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Ficure 1. Functional diagram of the optical system for the genera- 
tion of subsecond duration sawtooth shape radiance pulses. 


eliminated errors that result from photomultiplier fatigue, 
etc. 
Pulse shape as a function of time for various dia- 
phragm geometries is derived in the following paragraphs. 
The effective cross-sectional area of the light beam, 
when the disk circumference is at point x on the dia- 
phragm, may be expressed as 


x 
A, =|, y dx 


where x = ro — r and y is the width (height) of the 
diaphragm at x (fig. 2). 


(2) 


For rectangular geometry, y = constant = a, thus 


Ae = a. (3) 
For triangular geometry, y = bx, thus 
A, = bx?/2. (4) 


When the disk rotates «t a constant speed, the angle of 
rotation is proportional to time ¢, that is 


a= Kt. (5) 
From eq (1) and the definition of x, it follows that 
Xx=ro —r= ka. (6) 
Combination of eqs (5) and (6) yields 
<= hits. (7) 


The radiance L as seen by the pyrometer at a given 
instant is 


L= i A. (8) 
Combination of eqs (3), (4), (7), and (8) yields the 
following: 
For rectangular geometry 


L = Ke; (9) 
For triangular geometry 
L=K’? (10) 


where K and K’ are combination of constants. Equations 
(9) and (10) represent radiances that are linear and 
quadratic functions of time, respectively. 





Ficure 2. The attenuator disk and the diaphragm opening. 
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3. Experimental Results 


In order to assess the operational characteristics of the 
system, experiments were performed using diaphragms 
of rectangular and triangular openings. These experi- 
ments are summarized in the following paragraphs. 

Linear radiance pulses: A total of 15 experiments were 
conducted at five different attenuator speeds yielding 
radiance pulses (increasing linearly with time) in the 
range 20 to 150 ms. Additionally, four experiments were 
performed with the direction of rotation of the attenuator 
reversed, This scheme generated radiance pulses that 
decreased linearly with time. 

Quadratic radiance pulses: A total of four experiments 
were conducted at two different attenuator speeds. An 
additional experiment was performed with the direction 
of rotation of the attenuator reversed. 

Oscilloscope trace photographs of representative radi- 
ance pulses are shown in figure 3. The radiances as de- 
tected by the high-speed pyrometer and recorded with 
the digital data acquisition system were fitted to linear 
and quadratic functions, depending on the pulse shape, 
using the least squares method. A summary of the experi- 
mental parameters and results is given in table 1. It may 
be seen that average standard deviation (individual point) 
of the results of both pulse shapes is approximately 0.6 
percent corresponding to data covering a radiance ratio 
of over 10. 

The above value for the standard deviation reflects 
departure from ideal behavior that may be due to several 
factors, such as random noise in radiance samples and in 
the pyrometer operation, nonuniformity in the lenses and 
in the optical alinement, nonuniformity in the attenuator 
circumference, variation in the attenuator speed, etc. In 
order to separate the contribution of the pyrometer and 
the random noise inherent in radiance measurements 
from the attenuator and the associated optical system, 
two separate experiments were conducted under steady- 
state conditions. In one of the experiments, the attenuator 
was positioned to provide the maximum opening of the 
diaphragm as seen by the pyrometer. Then the pyrometer 
response to the radiance source was recorded. The stand- 
ard deviation of the results was 0.2 percent. A similar 
experiment was performed with the attenuator positioned 
to provide radiance equivalent to that of the lowest level 
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Ficure 3. Oscilloscope trace photographs of linear (upper) 
and quadratic (lower) radiance pulses. 


Horizontal scale: 10 ms (upper), 20 ms (lower) per major division. Vertical 
scale: arbitrary units. Dots forming the long horizontal lines correspond to radi- 
ances from the reference source. 


used in the dynamic experiments. The standard deviation 
of the results of the pyrometer output was 0.8 percent. 
Assuming an approximately linear variation of random- 
ness between the two extremes, a value of 0.5 percent 
for the standard deviation contributed by the pyrometer 
and radiance measurements in general is obtained. Thus, 
if it is assumed that total standard deviation is the square 
root of the sum of the squares of individual standard 
deviations, the standard deviation corresponding to the 
nonuniformities of the attenuator and the associated opti- 
cal components becomes approximately 0.3 percent. This 
value is comparable to the uncertainty with which the 
pyrometer linearity could be determined in separate 
experiments, 

In order to ascertain the fact that the value obtained 
for the standard deviation was the result of random de- 
viations in the data, the data on radiance were also fitted 


TABLE 1. Summary of experimental parameters and results 








Direction of 








~~ Pulse oo Number of Average Average total 
Item sala radiance length : : : points per radiance standard 

memes variation (ms) ins icienmaans experiment ratio deviation, % 
1 Linear Increasing 23 3 20 10.3 0.6 
2 Linear Increasing 30 3 26 12.1 “a 
3 Linear Increasing 39 3 33 9.0 6 
4 Linear Increasing 59 3 50 118 A" 
5 Linear Increasing 144 3 123 11.3 6 
6 Linear Decreasing 23 1 20 10.0 7 
7 Linear Decreasing 30 1 26 9.5 my | 
8 Linear Decreasing 39 1 33 10.4 6 
9 Linear Decreasing 59 1 50 10.6 a 
10 Quadratic Increasing 60 3 39 11.0 6 
ll Quadratic Increasing 40 1 26 10.1 5 
12 Quadratic Decreasing 60 1 39 11.1 6 
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to functions of one and two higher degrees than the 
pertinent one. No improvement in the standard deviation 
was observed, indicating that there was no detectable 
systematic curvature imposed on the proper function. 


4. Conclusion 


The results of the experiments described above have 
shown the feasibility of generating subsecond duration 
pulses of prescribed shapes (linear and quadratic) with 
a standard deviation of less than 0.5 percent over a radi- 
ance ratio of 10. The method allows one to insert ap- 
propriate diaphragms in the optical system to obtain 
radiance pulses of various shapes. For example, in addi- 
tion to linear and quadratic, it is possible to generate 
pulses of cubic, exponential, and other shapes. 
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The scheme of generating radiance pulses allows one 
to detect the linearity of radiation detectors under dy- 
namic conditions, and also provides the experimenter 
with means of simulating high-speed heating or cooling 
experiments. 
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A modified transient pulse technique for measuring shear and longitudinal sound speeds in 
solids has been devised. The technique is described and evaluated for both precision and accuracy 
on a variety of solids. This evaluation has been done by experiments in which the constancy of 
sound speed with path length is used as the prime test for accuracy. A number of variables such 
as transducer frequency, transducer size, choice of different coincidence points, etc., have been 
examined. The transducer bond problem is examined and overcome by the use of a time delay 
bond. Indicated accuracy of the method is found to be better than one part in 10*. Advantageous 
sample geometries for accuracy testing are discussed. 


Key words: Acoustical pulse; longitudinal and shear sound waves; solids sound speed accuracy; 


speed of sound; transient pulse-acoustical, 


1. Introduction 


For the purpose of this discussion, measurement meth- 
ods for sound speed can be divided into two groups: 
(1) equilibrium-monochromatic sinusoidal methods (ES), 
and (2) nonequilibrium sharp-transient pulse methods 
(NET). For the ES methods, error evaluation depends 
upon their described character. The specification of the 
particular equilibrium condition is important. Many vari- 
ations of the ES method are discussed in the extensive 
works of H. J. McSkimin and others [1—5].! But in the 
NET methods the predominant notion is that the first 
sound arrival proceeds by the shortest path and travels 
with the sound speed of the infinite medium. The nearer 
to the first arrival that the measurement is done, the less 
significant are the various end effects and unwanted 
boundary reflections. 

Because the NET method performs poorly for materials 
that have frequency dispersion or high losses, its use has 
been limited in recent years. The monochromatic ES 
systems are more applicable to these cases. Despite these 
limitations, various NET methods have been used suc- 
cessfully in this laboratory [6-8]. Recently, the latter 
has been applied to the measurement of speed of sound 
in solids with surprisingly good results. Modifications of an 
older technique have largely overcome the remnant sound 
energy problem which occurs with continuous pulsing. 

The method uses a pair of very short pulses, each ap- 
proximating a delta function. The pulses are separated 
by a time interval which is precisely adjusted to cor- 
respond to the round-trip time of flight for sound through 
the sample. The signal pair is followed by a dead time 
which is sufficiently long enough for the residual sound 
to have died out. During this dead time, the transducer 
is used as a receiver in such a manner that the n™ return 





1 Figures in brackets indicate the literature references at the end of this paper. 
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of the second pulse is superimposed on the (n + 1)™ 
return of the first pulse. The earliest recognizable feature 
of these returning signals is used to make this coincidence. 
This is related to the most direct path for sound travel. 
This contrasts sharply with the ES methods. In general, 
the phase lag of the acoustical signal is measured by 
comparison to the infinitely long wave train from which 
the initial pulses are chopped. These ES experiments are 
essentially monochromatic, as the propagated pulses are 
made up of many wavelengths in signals of constant 
amplitude. 

The method is discussed with special emphasis on pre- 
cision and accuracy as determined empirically. The uni- 
formity of the sound speed and the similarity between 
the two pulses used for coincidence were used as tests 
for accuracy. The effects of different variables on the 
sound speed have been examined and are noted. Some 
guidelines (shape and material) for the choice of a suit- 
able sample for testing the accuracy of the method have 
been developed. 


2. Experimental Features 


This transient video pulse technique is used for sound 
speed measurements in short fat solid samples (3.5 cm 
long X 3.8 cm diam) and can be used either single 
ended or with separate sender and receiver. The single 
ended way produces the best results and is used for this 
study. The transit time for a round trip is the principal 
measurement made. 

With continuous pulsing at intervals comparable to the 
transit time of the sample, remnant sound will build up 
in small samples of high Q and coincidences for timing 
the transit time will be impaired. The remnant sound 
problem is handled by the use of a special sequence of 
driving pulses. The transducer is driven by two very short 
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Ficure 1. Experimental setup schematic. 


(<30 ns) 50-V pulses, each approximating a delta func- 
tion, which are separated by a known and adjustable time 
period, r. These pulses are followed by a long dead time 
(2" — 1)r and then the full sequence is repeated. The 
period r is determined by the frequency of the control 
oscillator (B) as measured by the frequency counter (G) 
(see-fig. 1). This drives the two-shot apparatus (A) which 
produces the pulse sequence. 

‘The«dead time generated by the two-shot apparatus 
can be adjustable to any value (2" — 1)7 where n is an 
integer ranging from 1 to 8. Such adjustments allow for 
variations in sample losses by permitting the dead time 
to be adjusted so that, at the beginning of the next 
sequence, the left-over residual sound is negligible. Most 
of these measurements were made with n = 8 and with 
a dead time of 255 r. The two-shot apparatus is a parallel 
set of 1 MHz binary chains having an adjustable phase 
difference for pulses traveling down the two chains [9] 
(see fig. 2). A provision which allows interchanging the 
firing order of the respective chains is used to cancel 
slight electronic phase differences between the outputs of 
the two chains. This function is indicated by the “A 
before B” switch on figure 2. 

The pulses are then amplified and fed to a blocking 
oscillator which drives the transducer. The transducer, 
acting as a highly damped resonant sender, sends the 
sequence of sound pulses into the sample (H), figure 1. 
In this system a single pulse is reflected back and forth 
in the sample creating a whole family of successive 
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Ficure 2. Two-shot pulse generator — a block diagram. 
A before B switch sets the chain firing order, namely A before B or B before A. 
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acoustic reflections: first round trip reflection, second 
round trip reflection, etc. In this way the two pulses as 
generated by the apparatus form two overlapping families 
of acoustical reflections. The superposition of all the 
acoustical reflections is received by the transducer. The 
electrical output of the transducer is amplified by a low 
noise amplifier ((D) of fig. 1) and then the coincidences 
are displayed on a delayed sweep 100 MHz scope (E). 

Coincidences are set between the two overlapped fam- 
ilies of reflected pulses. With the control oscillator, the 
period r is adjusted until the second reflection generated 
by the first pulse of the pair coincides with the first re- 
flection generated by the second pulse. Specifically this 
coincidence is made by aligning the first half cycles of 
the two reflections. Usually the combination of the two 
first half cycles is maximized to indicate alinement. This 
represents a coincidence made at a time between 14 to % 
of the resonant period of the transducer after the receipt 
of the very first arrival (earlier coincidences can be 
obtained photographically). The definition of the coinci- 
dence setting is improved by a low frequency phase modu- 
lation of the control oscillator (fig. 1 [F] modulation 
oscillator). 

Because the two-shot apparatus affords a very low 
pulse rate, coincidences are especially easy to set even 
for materials of the highest Q. This condition is due 
primarily to the absence of driving energy on the trans- 
ducers at the time of coincidence and to the lack of 
acoustical reverberation from the previous pulses afforded 
by the long dead time. 

Very thin barium titanate transducers are used for both 
longitudinal and shear wave generation and reception 
(see fig. 3). A backing mixture of epoxy resin and 
tungsten graded powders is attached to the transducers. 
The high loss characteristics of this mixture helps keep 
the Q of the transducer system below 10, which is desira- 
ble. This reduces long ringing and improves the initial 
pulse edge. 

The transducers are coupled to the sample by an 
acoustically thick joint. Joints 10 to 25 um in thickness 
are constructed of materials having an acoustic impedance 
much lower than that of either the transducer or sample. 
This helps isolate the sample from the transducer and 
the reflections of sound energy in the sample tend to take 
place at the proper surface. These bonds have the ad- 
vantage that they are simple and easy to construct. 
Compressional bonds are usually made with water or 
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Ficure 3. Typical transducer and bond construction. 
































alcohol at room temperature and 4-methyl-1-pentene [10] 
or 2-methyl butane [11] at low temperature (77 to 196 
K). A good shear coupling agent at room temperature is 
Dow Resin 276-V9 [1] (poly-methyl styrene), and at 
low temperatures is 4-methyl-1-pentene. 


3. Technique Evaluation 


Errors and accuracy have been evaluated by a number 
of different experiments for which various accessible 
variables can be controlled. Two criteria are used for 
evaluation: sound speed and round-trip transit time. 
Round-trip transit time is used primarily for studying 
small time delays associated with bond and end effects. 
The sound speed is measured whenever the “Crucial 
Test” is applied. 


3.1 The Crucial Test (Sound Speed Uniformity) 


The accurate determination of sound speed is hindered 
greatly by the lack of a uniquely definitive experiment. 
In general, the toughest convenient test for accuracy that 
is readily available is to examine the change of sound 
speed as a function of sample length. In liquids this is 
relatively simple. The physical limitations of solid sam- 
ples cause special problems. The following four are 
examples: (1) The length of a solid sample cannot be 
changed in a continuous fashion. (2) Material uniformity 
is always questionable. (3) Spurious reflections occur 
more readily in the smaller physical sizes commonly used 
in solid measurements. Finally (4) there are difficulties in 
bonding the transducer to the sample. The importance 
of the constancy of sound speed with length cannot be 
overemphasized. Although it is not a sufficient test, it 
must necessarily be applied whenever accuracy is of 
interest. 

This test was applied for three different materials. The 
first sample used was a family of aluminum gage blocks 
(5 lengths °4 in thru 4 in long X 2 in diam). The 
blocks were cut consecutively from a_ polycrystalline 
mother bar. Under the best measurement conditions, the 
longitudinal speed of sound was found constant within 
+2 X 10°‘ parts (18 MHz resonant transducer) and the 
shear speed within +1 107° parts (8.5 MHz resonant 
transducer) for the five lengths. These errors in sound 
speed for the aluminum sample could easily be explained 
by the lack of radial stress uniformity in the bar samples 
and by nonuniform sound scattering caused by the poly- 
crystalline state. 

A search for a better material was begun. Commercial 
grade fused silica was tried and found to be poor for the 
purpose. Local variations in the internal stress and density 
as seen in a polarimeter could be roughly correlated with 
variations in sound speed measured locally. Fused silica 
of more uniform optical density gave better uniformity 
in the sound speed. Finally a successful material was 
found in a borosilicate optical crown glass manufactured 
by the National Bureau of Standards under the name 
BSC 517/645. Some physical properties of this glass are 
given in the Smithsonian Physical Tables [12]. This 
material was found to have high optical and acoustical 
uniformity. At room temperature this glass has the 
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bonus property of a near zero velocity temperature 
coefficient. 

A test sample of this material was fabricated so as to 
have three different path lengths (a parallelepiped of 1 
in X 1.5 in X 2.0 in). The question of acoustical uni- 
formity is then isolated to local sound speed uniformity, 
and this has been determined to be below the threshold 
of measurement for this method (< +5 X 10° parts). 
Results for compressional waves, using a transducer 
resonant at 18 MHz on the BSC 517 sample, showed a 
barely measurable increase in speed for shorter path 
lengths [V(1 in) >V(2 in) by ~ 10° parts]. As will 
be discussed later, adjustments that reduced the bond 
error made the sound speed on this sample constant to 
the precision of the method. The same sample was used 
to test the shear mode (8.5 MHz resonant transducer) and 
the shear speed was found to be constant to within the 
precision of +5 X 10°° parts. 

Using different path lengths, another material was 
measured. With a glass optical absorption cell (nominally 
2 in by 3 in) the longitudinal (18 MHz resonant trans- 
ducer) sound speed of pure water was measured along 
two path lengths (a third length was obtained from 
another cell). The transducers were fixed to the outside 
walls of the container in the same manner as they are for 
any solid. Thus solid bonding problems are duplicated 
here although these effects are smaller than normal be- 
cause of the larger transit times involved. The temperature 
environment was well isolated and highly controlled. A 
platinum resistance thermometer was used to measure 
temperatures to a millidegree. Results on three path 
lengths (6.3, 8.9, 14.0 cm) indicate that the constancy 
of sound speed is better than the total precision (+ 5 X 
10-5 parts) and that the measured value agrees extremely 
well with the recent work of Carnvale et al. [13]. 


3.2 Effect of Diameter of Transducer 


The effects of changing transducer diameter with con- 
stant frequency were investigated. Very thin paper masks 
(25 wm) were used to restrict the radiating area of the 
transducer. Bonds were made by filling the hole in the 
mask with the bonding fluid. (It was determined that no 
sound could be transmitted through the wetted mask 
area.) Figure 4 shows how the mask hole size is related 
to the measured transit time through the sample. This 
experiment was performed on fused silica using three 
different resonant transducers. In general, the transit time 
increases (speed decreases) for smaller effective trans- 
ducer areas. This behavior might be explained by the 
notions discussed in 4.1. The abrupt change of speed from 
the unmasked value to the largest masked value is due 
primarily to the thick bond effect which will be discussed 
later. When the same experiment was performed using 
shear waves, very little change in transit time was 
measured (less than 10°* parts). 


3.3 Effect of Resonant Frequency of Transducer 


Early experiments on fused silica indicated that the 
measured sound speed tends to decrease with increasing 
resonant frequency of the transducer. A similar effect was 


found with the BSC 517/645 sample. Figure 5 shows the 
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Ficure 4. Effect of mask diameter on transit time. 


3.5-cm long fused silica sample with each of three transducers of different 
frequency. Water forms the bond. Slashes mean ‘‘no mask.” 
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Ficurs, 5. Sound speed (L-waves) in borosilicate crown glass BSC 
517/645 for three path lengths and four transducers of different 
frequency. 

Tight water bond. Flags indicate precision. 


constancy test of sound speed as performed with the three 
lengths of the borosilicate sample using four longitudinal 
transducers of different resonant frequency. The resonant 
frequency of each transducer is marked on the appro- 
priate curve. Note that as the resonant frequency in- 
creases, the speed decreases and becomes more uniform 
with path length (maximum spread for all transducer 
frequencies and sample lengths is 2 X 10°*). These 





measurements were made with a tight bond (water bonds 
of less than 10 pm). 


3.4 Earlier Coincidences 


As already mentioned in the introduction, it is felt that 
the closer to the first sound arrival that one sets the timing 
coincidences the better the resulting accuracy will be. One 
way to facilitate this is to use transducers of higher 
resonant frequency as was discussed in section 3.3. An- 
other is to use earlier recognizable features of the tran- 
sient pulse. Efforts were made to examine the effects on 
the transit time of the position of the coincidence setting 
in the received pulse. Transit times for longitudinal sound 
through a fused silica sample were measured for nine 
successive half cycles. The coincidences were set by maxi- 
mizing each half cycle starting with the first and working 
back to the ninth half cycle in the pulse. These results are 
shown in figure 6 for two different resonant frequency 
transducers. The various half cycles at which coincidences 
were made are appropriately marked and the position of 
each related to the beginning of the pulse is measured on 
the abscissa in us. This test was also performed using 
shear waves (8.5 MHz resonant transducer) with very 
similar results. 

An arbitrary scheme of obtaining information about 
coincidences earlier than the first half cycle can be car- 
ried out photographically. Photos of the front edge of the 
first half cycle were obtained from a 100 MHz fast sweep 
scope. Pictures were taken for both the first and second 
reflections separately. Amplitude measurements were 
taken from these as a function of time back to the begin- 
ning of the signal from the first half cycle maximum. A 
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Ficure 6. Transit time (L-waves) versus time from the beginning 
of the signal at which coincidence was made. 


Two transducers (18 and 7 MHz resonance frequency) on 3.5-cm long fused 
silica sample. Tight water bond. 
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Ficure 7. First quarter cycle of first and second reflections 
normalized at the first % cycle maximum. 


Ar is excess delay for points of equal amplitude (arbitrary time zero at first 


% cycle max.) L-waves in 2.5-cm long borosilicate crown glass. Transducer, 
7 MHz. Tight water bond. 


normalized amplitude graph demonstrates any shape 
change between the front edges of the first half cycle for 
the first and second reflections. 

Figure 7 is one such normalized comparison. The time 
reference origin is fixed at the first maximum of each and 
negative time proceeds forward into the pulse in the nor- 
mal way. Note that in this case the second reflection rises 
faster than the first (indicative of some geometric inter- 
ference effects). The time difference between equal ampli- 
tude points of the two pulses is plotted as Ar, the pulse 
shape parameter. This is associated with pulse shape 
change and has been used to correct the first half cycle 
coincidence data to give better values of sound speed. As 
one goes forward into the pulse the precision of Ar gets 
steadily worse making values of Ar earlier than 30 ns 
from the first half cycle peak undependable. 

Using Ar to correct the transit times, corrections can 
be found for the sound speed. These corrected transit 
times are shown by the small A data and the e data for 
the two different resonant frequency transducers on figure 
6. Note that these data tend to converge to the same value 
at the pulse beginning. Further confirmation for this arbi- 
trary correction is obtained from the favorable results of 
sound speed uniformity tests shown in figure 8. Figure 8 
is figure 4 corrected for the early pulse shape variations 
using Ar for the three different path lengths and the four 
resonant transducers on the BSC 517/645 sample. Note 
that not only is the sound speed more nearly constant for 
the three lengths, but the data for the different resonant 
frequency transducers tends to cluster around the 18MHz 
transducer value (maximum spread for all lengths and 
frequencies is 2.5 10°‘). In the case of shear waves, 
the difference in shape of the first two reflections is un- 
detectable and there is no reason for such a correction. 

The following general observations can be made about 


the behavior of Ar, pulse shape, and sound speed. 
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Ficure 8. Speed of sound (L-waves) in borosilicate crown glass 


BSC 517/645 for three path lengths and four transducers of dif- 
ferent frequency. 


The data are those of figure 4 corrected according to figure 7. 


Flags indicate 
precision, 


(1) The higher the resonant frequency of the trans- 
ducer the smaller the pulse shape factor. In fact for the 
18 MHz resonant transducer Ar is essentially zero. 

(2) The smaller Ar is (or the more nearly alike the 
first and second reflections are) the more uniform is the 
sound speed for the three path lengths. 

(3) By applying the corrections for the shape factor 
to the sound speed, the constancy is improved and the 


values of sound speed, as corrected, bunch closer to the 
highest frequency value. 


3.5 Transducer-Sample Bonds 


The effects of a thick bond on the transit time for 
longitudinal waves were studied by using a variable water 
bond. Apparatus that allowed continuous fine adjustment 
of the bond thickness while maintaining parallelism be- 
tween the transducer and sample was used. Changes in 
the transit time were measured as a function of bond 
thickness for a number of different materials, a number 
of path lengths and different frequency transducers. Fig- 
ure 9 is typical. It shows the increase in the measured 
transit time with bond thickness for three transducers of 
different resonant frequencies. The three path BSC 517 
sample was used. Measurements along these three path 
lengths give essentially the same increase in transit time, 
indicating that the bond errors are independent of sample 
length. Furthermore, using transducers of higher resonant 
frequency reduces the magnitude of the bond error. Sim- 
ilar work using the 8.5 MHz shear resonant transducer 
with a viscous resin bonding material (Dow 276-V9 
Resin) was done. A similar but much smaller increase in 
transit time was obtained (increase was ~ 1 X 1074 
parts). 
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Ficure 9. Change in transit time versus water bond thickness. 


L-waves in borosilicate crown glass BSC 517/645. Three transducers of different 
frequency. Data for three lengths on each curve. 5 is the half period of the 
transducer frequency. 


Generally as the thickness is increased to that of the 8 
bond (thickness corresponding to the time between the 
first arrival and the first maximum) both the constancy 
of the sound speed and pulse shape factor are improved. 
In fact when the sound speed is measured on the three 
length BSC 517 sample using 6 bonds, in a similar man- 
ner to figure 4, one finds that all the transducers of 
different resonant frequencies give more constant values 
of sound speed and the values obtained cluster much closer 
to each other than for a thin bond (maximum spread for 
all frequencies and lengths is 1.5 & 10~*). At the highest 
transducer frequency of 18 MHz, measurement of sound 
speed with a 5 bond gives values constant within +5 X 
10°* parts which is lower by 10~* parts than the data of 
figure 4, taken with thin bond. 


3.6 Damping Backing to Transducer 


Interference in the leading edge of the received pulses 
which occurs because of reflections off the back of the 
transducers was investigated. A number of different back- 
ing mixtures were tried in an effort to reduce transducer 
ringing. It was found that lowering the Q of the trans- 
ducer by use of a backing with higher damping improved 
its performance. Both the pulse shape and the constancy 
of sound speed were improved by a tungsten-epoxy back- 
ing which lowered the Q to less than 10. At 18 MHz pulse 
shapes were made nearly identical by the use of these 
improved backings. 


4. Discussion 


4.1 Expected Interference Errors 


It is difficult to estimate theoretically the magnitude 
and the sign of errors in sound speed unless sufficient 
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knowledge about the particular measurement is available. 
In these experiments, coincidences are made between 
corresponding points on the front edge of the first and 
second reflections. Depending on the kind of distortion 
acquired between these reflections, error of either sign 
may be expected in the transit time. 

Figure 10A represents the condition observed experi- 
mentally by this method at the lower frequencies. This 
shows the leading edges of the first and second reflections 
as they are related to the first time of arrival. The first 
reflection rises more slowly to the peak than does the 
second. The contrary situation has never been observed 
in all the experiments performed. Let the measured transit 
time, X, be given by 


X=1+ Ate — At, 


where r is the shortest time of flight, associated with first 
arrival. Because Ato, the time between the first arrival 
and the first maximum of the second reflection, is smaller 
than At,, the corresponding time for the first reflection, 
the measured transit time will be too small when coinci- 
dences of the first half cycle are used. Thus, improving 
the similarity between pulses of the first and second reflec- 
tions should cause the transit time to increase. This is 
what is observed experimentally. In fact, it occurs in con- 
junction with the improvement of sound speed constancy. 
The following simplification may help to explain the 
observed differences in the first and second pulse 
reflections. 

Assume that for the leading edge of the signal the 
sender acts more like a point source than like a plane 
wave generator. Figure 10B shows the longest and the 
shortest rays that can connect a point sender and a finite 
planar receiver for the first two reflections. The difference 
between the longest received and most direct rays is 
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Ficure 10-A. Front edge of the first and second reflections. 


Time versus amplitude. First arrival is at the origin. 
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Ficure 10-B. Extreme rays connecting a point sender to a finite 
receiver for two consecutive reflections. 


greater for the first reflection than for the second. Thus 
one might expect the delay between the first arrival and 
the initial maximum to be greater for the first reflection 
than for the second. Although this model fits the qualita- 
tive behavior, it oversimplifies a much more complicated 
geometric interference problem. It is for this and similar 
reasons that the empirical test has become the guideline 
for this experiment. 


4.2 Bond Errors 


In this method one of the biggest sources of error con- 
cerns multiple reflections between the two faces of the 
bond. When sound is sent through a bond of finite thick- 
ness some of the energy is transmitted and some internally 
reflected. In such a manner a whole family of pulses 
separated by some multiple of twice the transit time of 
the bond will be transmitted into the sample. The super- 
position of these reflections changes the shape of the 
observed pulse. The resultant shape depends on how the 
intrabond reflections add to the straight-through pulse. 
The leading edge of the pulse may be sharpened, length- 
ened, or left essentially unchanged. A similar process will 
occur at the encounter of the sound with the bond after 
each round trip. 

Not all of the encounters are of the same kind. The 
most important of these is the one after the sound has 
completed the first round trip. Part of the sound is trans- 
mitted through the bond and the rest is reflected back into 
the sample. With a thin bond, there are multiple reflec- 
tions at the faces of the bond which modify both the 
reflected and transmitted waves. Experiment indicates 
that the signal received by the transducer after one round 
trip is lengthened by more than the geometrical effect 
discussed above. On the other hand, the signal received 
after two or more round trips is modified less by the 
bond. 

To avoid the influences of these intrabond reflections, a 
thick bond is used. This “8 bond” delays the unwanted 
reflections so coincidences formed from the direct pulses 
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will be unhindered by later reflections. Experimentally, 
when the 8 bond is used, one finds that the constancy of 
sound speed is improved, that the shape differences be- 
tween first and second reflections become less and that 
the transit times are longer (measured speeds are lower) 
than for a thin bond. 


5. Summary 


The transient pulse technique of sound speed measure- 
ment has been shown to be highly accurate when used 
on nondispersive materials. Accuracies of the order of 
the precision, which is +5 & 107°, have been realized. 
The constancy of sound speed was used as the measure of 
accuracy. The accuracy of this technique is affected by 
a number of variables. The observed sound speeds are 
more nearly constant when using: (a) higher resonant 
frequency transducers, (b) earlier data in the pulse to 
set coincidences, (c) thick buffer bonds of 8 thickness, 
which remove interference caused by the intrabond re- 
flections, and (d) transducers of smaller diameter-to- 
thickness ratio. The obvious has been confirmed in that, 
as the sound speed constancy is improved, the leading 
edges of the first two reflections become more identical. 
Local acoustical uniformity of the test samples is very 
important and the NBS optical glass BSC 517 fills these 
requirements quite well. The choice of a multipath sample 
as opposed to many samples of different length removes 
much of the question of sample uniformity. 


The author is deeply indebted to Martin Greenspan 
and Frank Breckenridge for discussions about the acousti- 
cal aspects, to Carl Tschiegg for help with certain elec- 
tronic features, and to J. Wilson and C. Smith for typing. 
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Factors affecting the (unwanted) seismic response of infrasonic microphones are indicated. 
Past measurements of ground motion deduced from the radiated atmospheric sound measured with 
infrasonic microphones are reviewed, and such measurements are compared with seismometer meas- 
urements of ground motion. Seismic motions caused by the Japanese earthquake of May 1968 are 
used in this example. The seismic response of the infrasonic microphone used for this measurement 
was experimentally determined and the results are presented. A simple method of compensating for 
interfering seismic effects on microphones is described. 


Key words: Compensation; ground motion; infrasonic microphones; seismic response. 


1. Introduction 


If seismic motions caused by a distant earthquake are 
strong enough, an earthquake-associated air-pressure 
wave may be observed. These air-pressure changes are 
observed coincident with vertical motions of the earth 
in the immediate vicinity of the microphone, and are 
caused by radiation of sound into the atmosphere by the 
ground motions. Observations on such acoustical radi- 
ation have been reported previously [Cook and Young, 
1962; Donn and Posmentier, 1965]. The amplitudes of 
vertical ground motion deduced from the air pressure 
measurements have been compared with ground motions 
measured with seismometers operated in the vicinity of 
the infrasonic microphones [Cook, 1965; Donn and 
Posmentier, 1965; Golitsyn and Klyatskin, 1967]. In 
general quite good agreement has been obtained using 
the two methods. 

In addition to the infrasound caused by local ground 
motions, air pressure waves have been observed originat- 
ing from both the vicinity of the epicenter of an earth- 
quake [Benioff and Gutenberg, 1939; Mikumo, 1968] 
and from intermediate regions along the seismic wave 
propagation path [Benioff et al., 1951; Donn and Pos- 
mentier, 1965; Press and Ewing, 1951]. The Alaskan 
earthquake of 1964 provided examples of both these 
classes of earthquake-associated atmospheric sound waves 
[Donn and Posmentier, 1965]. 

An example of an atmospheric pressure disturbance 
generated by the local passage of seismic waves from 
distant seismic activity is the record shown in figure 1. 
The acoustic data from microphones in the vicinity of 
Washington, D.C., appear in the upper portion of this 
figure. The microphone array at a typical station is de- 
signed for measuring four principal characteristics of 
infrasonic waves: (1) the amplitude and waveform of 
the incident sound pressure, (2) the azimuth of arrival 
of the wave, (3) the horizontal phase velocity, and (4) 
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the dominant period of the wave. An infrasonic station 
normally has four or more microphones equipped with 
noise-reducing devices located at ground level, approxi- 
mately in the same plane and spaced about 7 to 10 km 
apart. Figure 2 shows the microphone layout for the 
Washington, D.C., station. 

Each line-microphone produces frequency-modulated 
voltages proportional to the sound pressure in the at- 
mosphere. The tones are usually transmitted by telephone 
wires to a recording location. At this central recording 
location these tones are demodulated, amplified, filtered, 
and recorded in analog form both as ink-on-translucent- 
paper traces, and on magnetic tape. The system pass-band 
is designated N7 and is shown in figure 3. 

The disturbance shown in figure 1 is related to an 
earthquake on May 16, 1968, with the epicenter under 
the sea floor 100 miles east of Hachinohe, Honshu, Japan 
(U.S. Coast and Geodetic Survey, Preliminary Determina- 
tion of Epicenter: 0 = 00" 48™ 55.4* UT, 40.84° N, 
143.22° E, h = 7 km, M, = 7.9). This disturbance was 
much smaller at Washington, D.C., than the disturbance 
following the 1964 Alaskan earthquake. It was thought 
that the smaller ground motions would make intercom- 
parison with seismic data more feasible. Note that this 
disturbance was distinguished from local pressure changes 
caused by turbulent motions of the wind, because of its 
presence at all of the microphones of a multipartite array. 

The analog acoustic record shown in figure 1 is a 
superposition of the recordings from the various micro- 
phone sites. The waveform coherence between the various 
sites is quite good. Note that the time marks are almost 
coincident at the chart speed of 34 in/min, indicating a 
high horizontal phase velocity. The disturbance was dis- 
tinguished from acoustic waves propagating in the at- 
mosphere parallel to the earth’s surface, by its higher 
horizontal phase velocity. Rayleigh waves on the earth’s 
surface have a horizontal phase velocity of about 3 km/s, 
compared with 330 m/s for acoustic waves in the 
atmosphere. 

The disturbance started at 0135 UT and was lost in 
noise at about 0214 UT. The observed periods ranged 
from 17 s to 40 s. The maximum amplitude occurred at 
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Ficure 1. Seismic disturbance of May 16, 1968. 


0153 UT and was 5.5 bars at a period of 19 s. The 
ground displacement computed from these pressure data 
was 4.2 mm. The ground motion determined from seis- 
mometers located in Washington, D.C., was 4.3 mm. The 
Georgetown University Seismometer Laboratory made the 
recordings of ground motion shown in figure 1 available. 
The short-period worldwide standard seismograph was 
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Ficure 2. Microphone site layout. 


42 


used in making the recording. The seismic data were 
corrected for system amplitude response before being 
compared with the infrasonic microphone records. Both 
measurements may be considered accurate to within 
+ 10 percent. 

R. K. Cook [1965] has shown that the air pressures for 
such a disturbance are due to an integrated effect of the 
sound radiated by Rayleigh waves over a large area, as 
opposed to a seismometer which measures motion at a 
point; hence, it seems reasonable to consider infrasonic 
microphone data in order to obtain a more representative 
picture of ground motion over an area of the earth’s 
surface. Infrasonic data will be useful for measurement 
of very large amplitude ground motions because the 
sound pressure levels will probably be larger than the 
pressure variations due to local winds. It is possible that 
for very large ground motions the dynamic range of 
sensitive seismometers will be exceeded, and then acousti- 
cal measurements can provide the desired measurements 
of the motions. But the microphones must perforce be 
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Ficure 3. Microphone system response. 














located on the surface of the earth, and they are therefore 
subjected to the latter’s seismic motions. It is therefore 
necessary to determine the effects of the microphone’s 
seismic motions on its electrical output. The effects might 
be small or negligible, but they must be known 
quantitatively. 


2. Seismic Response Considerations 


In attempting to test and evaluate the seismic response 
of infrasonic microphones several factors must be con- 
sidered. Rayleigh waves occurring at the location of the 
pressure transducer generate peak pressure changes in 
the atmosphere calculated according to the relation 


AP = pcwX 
where 
w = 2nf 
f = frequency of the wave 
X = double amplitude of ground motion 
AP = peak-to-peak pressure change in the atmosphere 
p = density of air 
c = velocity of sound in air. 


This pressure change would be measured by the micro- 
phone together with any other instrumental seismic re- 
sponse effects that might be present. A frequently used 
form of infrasonic microphone is one that uses the motion 
of a diaphragm to detect pressure changes. This is the 
type of detector that will be considered here. 

In addition to the pressure AP = pcwX, such a trans- 
ducer will respond to accelerations independently of any 
pressure changes occurring in the atmosphere. A trans- 
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Ficure 4. Schematic view of microphone. 


ducer such as that shown in figure 4 will show seismic 
response due to both the diaphragm inertia and the 
inertia of air columns in the unit. Also, when the trans- 
ducer is moved vertically through the atmosphere, a 
pressure change will be introduced whose magnitude 
will be equal to pgX. Pressure changes due to air motion 
past the microphone’s openings to the atmosphere could 
also occur. The various factors contributing in theory to 
the seismic response are summarized as follows: 


Bernoulli Pressures* 





— peX + 1/2 pV? +... 
Open Tube Diaphragm Local Closed Tube 
Air Column Effects* Inertia Radiation Air Column Effects 
..-pwhX + ©” pprass2A + pwcX wee + 1/2 pw*hX 
where: account. Effects of local sound radiation by the shaker 
g = local gravitational acceleration will be negligible. Note that the relation listed for air 


= displacement double amplitude 
h; = vertical height of microphone’s external air 
column 
ha = thickness of diaphragm 
h = length of microphone’s internal air column 
V = maximum velocity of motion 
o— 2rf 
p = density of air 
Pbrass — density of the diaphragm 
c = velocity of sound in air. 


In experimentally testing seismic response on a vibra- 
tion shaker, the effect of translating the transducer ver- 
tically through the atmosphere will have to be taken into 








* There are no contributions due to the hydrostatic equation or Bernoulli pres- 
sures during seismic wave transit. These two effects are important during an 
experimental evaluation of the microphone with a vibration shaker. The relation 


shown for open tube air column effects is valid only during vibration shaker 
type experiments at low frequencies. 
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column effects is valid for an open-ended tube at low 
frequencies. The closed tube relation is half of this value. 
Microphone response to vertical displacements was tested 
at the National Bureau of Standards at Gaithersburg, 
Maryland, where Mr. James R. Houghton made available 
a shaker capable of producing approximately sinusoidal, 
vertical motions. The tests were performed at an ambient 
temperature of 22° C, a relative humidity of 20 percent, 
and an ambient pressure of 757 mm of Hg. The frequency 
of vertical motion was adjustable over the range 1 Hz 
to 0.01 Hz. The double amplitude was set to 4 in. 

The microbarograph system consisted of a pressure 
transducer [Cordero et al., 1957] equipped with a 7.5 s 
high-pass acoustical filter. The system pass-band including 
Discriminator N4 and Filter Amplifier N10 had 3 dB 
points at 1 and 20 s. Figure 4 shows a schematic view of 
the system microphone mechanical configuration. The 








microbarographic system was calibrated with a piston- 
phone, and a portable fixed-frequency pressure calibrator 
was operated before and after the seismic response tests. 
The pistonphone was a stainless steel wool-filled 50-gal 
volume equipped with a brass bellows variable-volume 
element to produce dynamic pressure changes. Theoretical 
considerations made it desirable to vary the key para- 
meters while subjecting the infrasonic microphone to 
vertical displacements on a shaker at a number of fre- 
quencies. Since the experiment was performed in the 
presence of local atmospheric pressure variations, a time 
of low atmospheric pressure noise was chosen. Further, 
it was necessary to choose vibration amplitudes large 
enough for the expected pressure variations to be larger 
than the local noise level. 


3. Experimental Results 


The system response on a vibration shaker table was 
found, within experimental error, to be independent of 
the orientation of the microphone inlet to the atmosphere. 
Hence, it is inferred that inlet air blast and Bernoulli 
effects did not contribute to the electrical output of the 
microphone. 

At longer periods than 5 s per cycle, the observed 
change in pressure closely approached the value to be 
expected from the hydrostatic term in the pressure equa- 
tion. There is no direct frequency dependence involved in 
this term. Figure 5 shows the pressure measured as a 
function of frequency. 

If the value for AP determined from the hydrostatic 
equation is subtracted from the observed pressure changes 
and the resulting data plotted as a function of frequency 
squared, the result is approximately a straight line in 
agreement with the predicted /* dependence. See figure 6 
for a log-log plot of AP versus frequency squared. 

The vertical external tubing length was varied in 5-in 
increments and the output of the system measured with 
the frequency of vibration kept constant at 1 Hz. The 
variation of pressure as a function of length of tubing 
below this valve is linear and changes as expected 
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Ficure 5. Experimental response of microphone as a function 
of frequency of vertical motion. 
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Ficure 6. Response of microphone as a function 
of frequency squared. 
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Ficure 7. Effect of air column length. 


(AP = 1.2 h, at 1 Hz, where h; is in inches; fig. 7). This 
is approximately the slope of the line in figure 7. 

Table 1 presents the contributions of the various seismic 
response mechanisms to the total pressure change od- 
served. The outputs have been computed at frequencies 
of 1 Hz and 0.1 Hz and for displacement amplitudes of 
0.001 in, 0.01 in, 0.1 in, and 1 in. Typical values are 
included in this table to show the relative importance of 
the various effective elements. The hydrostatic relation 
is also included for comparison purposes, although 
normally it would not be a factor in determining Rayleigh 
wave seismic response. The diaphragm compliance and 
effect of any backing volume sensitivity will reduce the 
system response, but these effects are not included in this 
worst-case consideration. 

Because local background seismic noise amplitudes can 
be over 20 wm, the displacement amplitude of 0.001 in 
is a reasonable reference to use. Table 1 shows that the 
“loudspeaker action” effect from seismic waves will domi- 

















TABLE 1. Seismic response contributions 
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Displacement X, inches 





mechanism 


X = 0.001 


X¥=0.01 X¥=0.1 x=1 





pgxX 0.003 uwbar 

1Hz wcpaird .638 wbar 
@'porassh2X  .002 ubar 
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pgx 0.003 ubar 
0.1 Hz wepairdX .0638 ubar 
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0.03 ubar 0.3 wbar 3 wbar 


.638 wbar 6.38 ubar 63.8 ubar 


w@*porassheX .2 X 107~*ubar 2 XK 107 ‘*ubar .002 wbar 02 wbar 
w*partX 3X 10~‘ubar 3 x 107‘*ubar .003 wbar .03 wbar 





hi = 
he = 
Porass — 


s=— 


nate. The effects due to transducer seismic response are 
typically less than 1 percent of the “loudspeaker effect” 
at the pertinent frequencies. Thus, for this infrasonic 
microphone configuration, it is quite reasonable to deduce 
ground motions from observed pressure changes without 
correcting for instrument seismic response. Also, with the 
proper length of air column termination, one can null out 
the seismic response of the microphone at infrasonic fre- 
quencies. In the NOAA infrasonic microphone, the return 
hose from the microphone to the noise-reducing array 
accomplishes this. The open tube relation used here for 
computing this form of seismic compensation is valid only 
for shaker-type experiments with negligible local radi- 
ation. A seismically compensated microbarograph has 
been described [Ewing and Press, 1953] for use at 
periods longer than 20 s. A different method of seismic 
compensation is used. 

Consideration of air column effects would suggest that 
caution should be exercised, for example, in installing 
an infrasonic microphone in a deep well and routing a 
tube to the surface. Since the air column seismic response 
is directly proportional to the length of the air column, 
a person might, in fact, be constructing a system sensitive 
to vibration. For 0.001-in tube displacements at 1 Hz, the 
seismically induced pressure amplitudes for an air column 
10 ft long would be 0.03 wbar. This is over three times the 
internal noise level of the microphone now used by 
NOAA. Noise due to local vibrations would be increased. 


10 in 

0.001 in 

8.7 g/cm* 
980.1 cm/s? 


I am grateful to Mr. Frank Baldwin, who did the art 
work on the illustrations, and to Rev. McConnell of 
Georgetown University, who provided the seismic records. 
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This report gives findings and recommendations developed during 
an investigation of performance characteristics for sanitary 
plumbing fixtures, conducted at the request of the Building Re- 
search Advisory Board of the National Academy of Sciences— 
National Research Council. The report describes the test methods 
that are recommended for the evaluation of 16 performance char- 
acteristics, and the nature of further work required to complete 
the development of four or five additional test procedures. 

The suitability of various existing test methods for evaluating the 
functional and performance characteristics of sanitary plumbing 
fixtures was investigated in the laboratory. In addition, new or 
modified tests for certain characteristics were developed. The 
laboratory work was performed only on bathtubs and flat speci- 
mens provided by industry through appropriate arrangements with 
the Building Research Advisory Board. Field inspection trips were 
made to provide the NBS project staff with up-to-date information 
on certain manufacturing processes and on installation and use 
problems. The complexities involved in the selection of valid 
performance levels are discussed, as well as the elements of judg- 
ment involved. A format that might be used in specifying perfor- 
mance is suggested for each test procedure, and the rationale 
underlying each suggested format is given. 
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This paper describes a mathematical analysis for determining the 
microwave current distribution in a barretter in a rectangular 
waveguide. This distribution, when used with another analysis 
which calculates substitution error for any given current distribu- 
tion will provide a previously missing step necessary for calibrat- 
ing microwave and millimeterwave barretters for absolute power 
measurements. 
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After successful acquisition in August of reflected ruby laser 
pulses from the Apollo 11 laser ranging retro-reflector (LRRR) 
with the telescopes at the Lick and McDonald observatories, re- 
peated measurements of the round-trip travel time of light have 
been made from the McDonald Observatory in September with an 
equivalent range precision of +2.5 meters. These acquisition 
period observations demonstrated the performance of the LRRR 
through lunar night and during sunlit conditions on the moon. 
Instrumentation activated at the McDonald Observatory in October 
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has yielded a precision of +0.3 meter, and improvement to +0.15 
meter is expected shortly. Continued monitoring of the changes in 
the earth-moon distance as measured by the round-trip time of 
light from suitably distributed earth stations is expected to con- 
tribute to our knowledge of the earth-moon system. 


Balcom, M. M., Influence of red and blue pre-adaptation on 
hue matching of purple samples, J. Opt. Soc. Am. 60, No. 1, 
118-121 (Jan. 1970). 


Key words: Adaptation; color; hue; vision. 


This experiment was designed to show the effect, if any, of chro- 
matic adaptation on hue matching of purple painted samples of 
different saturations. Five experienced observers were adapted for 
one minute to a red field, then asked to find the hue match for 
each of three test samples from among a set of samples represent- 
ing the complete hue circuit in 100 steps at middle saturation. 
The three test samples were of slightly greater saturation and 
representative of a range of purple samples. The entire procedure 
was repeated with blue and neutral pre-adapting fields, all under 
Source C illumination. Results indicated that pre-adaptation did 
influence the hue-match selections, the average red-blue adaptive 
shift being about one Munsell hue step. Observers made system- 
atically different hue matches for the same test sample in accord- 
ance with their ages. A method for determining graphically the 
state of adaptation at the time the hue judgments were made 
showed the chromaticity of the test samples to be at least as 
influential as either the neutral surround used or the pre-adapta- 
tion stimuli. 


Beatty, R. W., Editor, Commission 1 progress in radio meas- 
urement methods and standards (Triennial report of prog- 
ress), Radio Sci. 4, No. 7, 579-590 (July 1969). 


Key words: International scientific radio union report; measure- 
ment methods; radio measurements; radio standards; URSI XVIth 
General Assembly. 


A review of progress within the USA in Radio Standards and 
Measurement Methods, covering ihe period 1966-1968, inclusive, 
is presented. This review is a part of the Triennial Report of the 
U. S. National Committee of URSI (International Scientific Radio 
Union) to the XVIth General Assembly of URSI. 

The following topics dealing with the U. S. A. progress in Radio 
Standards and Measurement Methods are included: 1) Time and 
Frequency, 2) CW Power, 3) CW Voltage, 4) Noise, 5) Pulsed 
Voltage and Power, 6) Attenuation, Phase Shift, and Time Delay, 
7) Impedance, 8) Antenna Characteristics, 9) Field Strength, 10) 
RF Properties of Materials, 11) Broadband Measurement Tech- 
niques, and 12) Laser Measurements. 

The report does not include all the developments which have oc- 
curred in the U.S.A. in the period indicated, but rather has at- 
tempted to include the most significant developments in each field. 


Becker, E. D., Ferretti, J. A., Farrar, T. C., Driven equilibrium 
fourier transform spectroscopy. A new method for nuclear 
magnetic Resonance signal enhancement, J. Am. Chem. Soc. 


91, 7784-7785 (1969). 


Key words: Driven-equilibrium; Fourier-transform; high resolu- 
tion; nuclear magnetic resonance. 


A new method using a high-resolution pulsed NMR spectrometer 
is described by which the equilibrium magnetization of nuclei with 
long 7: may be restored rapidly in order to take full advantage of 
the Fourier transform technique. 








Bender, P. L., Alley, C. O., Currie, D. G., Dicke, R. J., Faller, J. B., 
Some implications for physics and geophysics of laser range 
measurements from earth to a lunar retro-reflector, (Proc. 
Conf. NATO Advanced Study Institute, University of Newcastle, 
Upon Tyne, England, Mar.-Apr. 1967), Chapter in The Application 
of Modern Physics to the Earth and Planetary Interiors, S. K. 
Runcorn, Ed., pp. 523-530 (John Wiley and Sons Inc., London, 
England, 1969). 


Key words: Astronomy; continental drift; geophysics; gravity; 
laser; lunar; moon. 


The technique of high-accuracy laser range measurements to an 
optical retro-reflector package on the moon is discussed. The 
information which can be obtained by this method about the 
hypotheses of continental drift and ocean floor spreading and 
about the motion of the pole is described. A method for checking 
on the constancy of the gravitational constant is also given. The 
accuracy with which information on all of these questions can be 
obtained from lunar range measurements is higher than for any 
other methods which have been suggested. 


Berger, M. J., Seltzer, S. M., Caleulation of energy and charge 
deposition and of the electron flux in a water medium 
bombarded with 20-MeV electrons, (Proc. Conf. High Energy 
Radiation Therapy Dosimetry under the auspices of the American 
Association of Physicists in Medicine, New York, N. Y., June 15- 
17, 1967), Ann. N. Y. Acad. Sci. 161, Art. 1, 8-23 (July 3, 1969). 


Key words: Charge deposition; electrons; energy deposition; flux; 
Monte Carlo; multiple scattering. 


A combination of Monte Carlo and analytical techniques is used 
to compute the penetration of fast electrons in water, taking into 
account multiple elastic scattering by atoms, multiple inelastic 
scattering by atomic electrons, and the production of bremsstrah- 
lung. Secondary electrons are also considered which result from 
electron-electron collisions, or from the absorption or scattering 
of bremsstrahlung. The calculation provides information about the 
spatial distribution of energy deposition, charge deposition and 
the electron flux (primary and secondary) down to a spectral 
energy of 400 cV. Representative results are presented, mainly 
for 20-MeV electrons incident perpendicularly on a semi-infinite 
water medium. 


Birky, M. M., Simultaneous recording of near-field and far- 


field patterns of lasers, Appl. Opt. 8, No. 11, 2249-2253 (Nov. 
1969). 


Key words: Laser; near-field and far-field diffraction; ruby and 
neodymium. 


A technique for simultaneous recording of near-field and far-field 
diffraction patterns at several exposures for a single laser pulse 
has been developed. A ruby laser and a neodymium doped glass 
laser have been investigated. The Nd** laser shows striking high 
order cylindrical mode operation as a result of thermal stress. The 
time resolved output of this laser shows the usual spikes about 
500 nenoseconds wide. When the cylindrical mode operation takes 
place some of the spikes consist of a large number of ultrashort 
pulses. 


Blandford, J. M., In-plant testing for wash-and-wear apparel 
performance and labeling, Third Annual Tech. Conf. The Ap- 
parel Research Foundation, Inc., Washington, D. C., Oct. 14-16, 
1969, 12 pages (Oct. 14, 1969). 


Key words: Apparel; care instructions; labeling, apparel; perfor- 
mance, apparel; properties, apparel; testing textiles; wash-and- 
wear. 


A distinction is made between washable garments, in general, and 
wash-and-wear apparel, in particular, with reference to the prop- 
erties associated with appearance, performance, and ease of care. 
The testing equipment, testing procedures, and evaluation materials 
required by the apparel manufacturer in his assessment of these 
properties and formulation of ease-of-care instructions for labels 
and hang tags, are discussed. Tabulations and charts are included. 
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Brady, E. L., The National Standard Reference Data System, 
Mater. Res. Std. 9, No. 10, 19-21 (Oct. 1969). 


Key words: Critical evaluation; information centers; information 
services; NSRDS; physical properties; standard reference data. 


The National Standard Reference Data System is a government- 
wide effort to give to the technical community of the United States 
optimum access to the quantitative data of physical science, criti- 
cally evaluated and compiled for convenience. This program was 
established in 1963 by the President’s Office of Science and Tech- 
nology, acting upon the recommendation of The Federal Council 
for Science and Technology. The NBS has been assigned respon- 
sibility for administering the effort. The general objective of the 
NSRDS is to coordinate and integrate existing data evaluation 
and compilation activities into a systematic, comprehensive pro- 
gram, supplementing and expanding technical coverage when nec- 
essary, establishing and maintaining standards for the output of 
the participating groups and providing mechanisms for the dis- 
semination of the output as required. The NSRDS is conducted 
as a decentralized operation of nation-wide scope with central 
coordination by NBS. It comprises a complex of data centers and 
other activities carried on in government agencies, academic in- 
stitutions and non-governmental laboratories. The independent 
operational status of existing critical data projects is maintained 
and encouraged. Data centers that are components of the NSRDS 
produce compilations of critically evaluated data, critical reviews 
of the state of quantitative knowledge in specialized areas and 
computations of useful functions derived from standard reference 
data. 


Branscomb, L. M., Truth in packaging of scientific informa- 
tion, Meas. Data 3, No. 5, 104-105 (Sept.-Oct. 1969). 


Key words: Reviews; science information; scientific literature. 


This is an edited version of a talk presented to the Precision 
Measurements Association meeting in Boulder, Colorado, August 
28, 1968. 


No new technical information is contained in the manuscript. 


Brauer, G. M., Termini, D. J., Burns, C. L., Characterization of 
components of dental materials and components of tooth 
structure by differential thermal analysis, J. Dental Res. 49, 
No. 1, 100-110 (Feb. 1970). 


Key words: Characterization of dental materials; components of 
tooth structure; dental materials materials; differential thermal 
analysis (DTA); differential thermal analysis of calcium phos- 
phates. 


Differential Thermal analysis (DTA) has been employed to char- 
acterize a variety of dental materials, components of tooth struc- 
ture, as well as certain phosphates related to the mineral com- 
ponents of tooth structure. The technique was applied to the rapid 
comparison of similar materials, quality control of ingredients 
and products used in restorative dentistry, determination of melt- 
ing ranges especially those of dental waxes and gold alloys, and 
studies of reactions occurring at elevated temperatures. The DTA 
technique is simple and requires only a few milligram sample. 
Thus it should prove especially useful for specification testing of 
melting ranges of dental materials, for establishing proper curing 
procedures, determining catalyst effectiveness and measuring stor- 
age stability in various environments. 


Braun, W., Carrington, T., Line emission sources for concen- 
tration measurement and photochemistry, J. Quant. Spectr. 
Radiative Transfer 9, 1133-1143 (1969). 


Key words: Absorption; analysis; gas; photochemistry; spec- 
troscopy; vacuum ultraviolet. 


When a resonance line is absorbed in its own gas, measurements 
of this absorption can be used to derive a value of the product, 
nf, of concentration of the absorber and f value of the transition. 
This line absorption process is also important in the study of 
fluorescence and photochemistry, where it is possible to produce 
atoms or molecules in a single quantum state without otherwise 








disturbing the system.. The interpretation and success of these 
applications of line absorption depend critically on the shape of 
the line emitted by the light source. This is strongly influenced by 
the optical depth in the emitting region, and by the inevitable 
presence of a reversing layer through which the light must travel 
on its way out of the lamp. This paper presents a simple model 
which can be used to estimate the effects of these properties of 
the light source on the threshold and sensitivity of measurements 
of nf, and on the power which can be delivered to absorbing atoms 
outside the lamp. Emphasis is on the general principles of lamp 
design and diagnosis, rather than on accurate description of a 
particular lamp. 


Brenner, F. C., Research for a uniform quality grading sys- 
tem for tires, Rubber Chem. Technol. 42, No. 5, 1446-1449 
(Dec. 1969). 


Key words: Bias ply; grading; passenger car; radial; speed capa- 
bility; tires. 


The National Traffic and Motor Vehicle Safety Act of 1966 re- 
quires the establishment of a uniform quality grading system. This 
paper discusses the needs for the system and the complexity of the 
problem. The system proposed here will depend on five properties; 
tread wear, traction, impact resistance, endurance, and wheel speed 
capability. 


Brenner, F. C., Barton, F. W., Cuts and cut growth in tires, 
Rubber Chem. Technol. 42, No. 5, 1462-1465 (Dec. 1969). 


Key words: Cut growth; cutting; road test; tires. 


Forty tires were subjected to severe road use on a course (approxi- 
mately 3500 miles) that included Belgian block, over highways at 
turnpike speeds while overloaded. The course was repeated until 
the tires were worn smooth. At the end of each course, the number 
of new cuts and the length of all cuts in grooves was determined. 
Ten miles of Caliche (sharp stones) road was included between 
the first and second measuring period. 

The tires experienced a total of 169 cuts. Forty seven (47) of these 
grew, 34 of which grew less than 0.11 inches. No cut grew con- 
tinuously throughout the test. The cuts grew but not enough to 
a tire failure except in one case where the cut exposed the 
abric. 


Brenner, F. C., Mandel, J., Simson, B. G., Research for a uni- 
form quality grading system for tires. Il. Wheel speed capa- 
pon test, Rubber Chem. Technol. 42, No. 5, 1450-1461 (Dec. 
1969). 


Keywords: Bias ply; grading; passenger car; radial; speed capa- 
bility; tires; winter. 


The wheel speed capability of a tire is defined by a test method 
which determines the speed at which the tire fails on a laboratory 
test wheel. Data is reported on over 100 different passenger car 
tires of all grades and types over a range of sizes. 

A scaling system is devised for this property. It is found that the 
system produces consistent results for tires of given manufacturer’s 
nominal grade across the size range tested and for samples pro- 
duced several months apart. 


Cassidy, E. C., Cones, H. N., Electro-optical observations and 
measurements of distorted high-intensity electric fields, 
Proc. Conf. Electrical Insulation and Dielectric Phenomena, Buck 
Hill Falls, Pa., Oct. 20-22, 1969, pp. 77-86 (National Academy of 
Sciences, Washington, D. C., 1970). 


Key words: Dielectric liquids; electro-optical measurements; elec- 
trostatic field measurements; high voltage measurements; kerr cell. 


The electric field between electrodes immersed in a dielectric liq- 
uid is often distorted by the presence of space change, thus making 
analytical determination of the field distribution extremely difficult. 
This work utilizes laser light for Kerr electro-optical analysis of 
distorted electric fields imposed when high direct voltages are ap- 
plied to nitrobenzene-filled Kerr cells. Two-dimensional fringe 
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images, similar to those achieved in photoelastic mechanical stress 
analysis, are observed. Regions of high electrical stress are de- 
tected directly from concentrations of a greater number of fringes. 
Field profiles observed with voltages as high as 90 kV applied to 
cells with various interelectrode distances are presented. Space- 
resolved meaurements of relative field strength, actual field 
strength (in V/cm) and potential are also derived. The measure- 
ment techniques employed are believed to be accurate within 1%. 


Cassidy, E. C., Cones, H. N., Use of expanded laser beam to 
analyze high-intensity electric fields, J. Soc. Motion Picture 
Television Engrs. 


Key words: Electric field measurements; electrooptical measure- 
ments; high voltage measurements; Kerr cell; laser applications; 
pulse measurements and dielectric liquids. 


A Kerr electrooptical technique that permits observation and map- 
ping of distorted high-intensity electric fields has been devised at 
NBS. Analysis of the field profile from two-dimensional fringe 
images, similar to those observed in photoelastic mechanical stress 
analysis, permits calibration of the Kerr System for measurement 
of high-voltage pulses, and determination of relative field strength, 
actual field strength, and potential. 


Catanzaro, E. J., Champion, C. E., Garner, E. L., Marinenko, G., 
Sappenfield, K. M., Shields, W. R., Standard reference mate- 
rials: Borie acid; isotope and assay standard reference mate- 
rials, Nat. Bur. Stand. (U.S.), Spec. Publ. 260-17, 70 pages (Feb. 
1970), 65 cents, SD Catalog No. C13.10:260-17. 


Key words: Absolute ratios; assay standards; boron; coulometry; 
isotopic standards; mass spectrometry. 


A precise coulometric titration method has been developed for the 
assay of boric acid. The method is capable of providing boric acid 
assays precise to within 0.0025 percent (1 sigma) which is an 
order of magnitude better than existing techniques. Using this 
method, a lot of high purity boric acid has been assayed and 
found to be virtually stoichiometric HsBOs. Its acidimetric assay 
value of 100.00 + 0.01 percent characterizes this material as a 
primary chemical standard for boron determinations. 

Multiple samples of boric acid were studied and characterized 
with respect to homogeneity and the effect of relative humidity on 
stoichiometry. Procedures for precise preparation, aliquoting and 
storage of boric acid solutions were devised. 

A mass spectrometric technique for the precise measurement of 
boron isotope ratios was developed. Single-filament tantalum-ribbon 
sources are used, and “B/*°B ratios are determined by measuring 
the relative abundances of Naz“BO.* and Nas°BO:* ions, at 
masses 89 and 88, respectively. The effects of various parameters 
such as sample mounting procedure, filament material, sample 
size, total sample composition and sodium concentration were 
studied and alternative procedures were evaluated. 

Absolute values were obtained for the isotopic abundance ratios 
of two boron isotopic standards, using surface emission mass 
spectrometry. Samples of known isotopic composition, prepared 
from nearly pure separated boron isotopes, were used to calibrate 
the mass spectrometers. The resulting absolute values are: SRM 
951, Boric Acid, “B/*B = 4.04362 + 0.00137; and SRM 952, 
Boric Acid, “B/*°B = 0.053199 + 0.000032. The indicated un- 
certainties are overall limits of error based on 95 percent con- 
fidence limits for the means and allowances for effects of known 
sources of possible systematic error. 


Cezairliyan, A., High-speed methods of measuring specific 
heat of electrical conductors at high temperatures (A re- 
view), High Temperatures—High Pressures 1, 517-529 (1969). 


Key words: High-speed measurements; high temperature; specific 
heat; thermodynamics; thermophysical properties. 


Needs and advantages of high-speed methods (subsecond duration) 
for the measurement of specific heat of electrical conductors at 
high temperatures are presented. The high-speed methods are 
described and classified according to ambient and operational con- 
ditions. Measurement and recording of experimental quantities 
are discussed. A historical development and the major features of 








high-speed experiments are summarized. A discussion and conclu- 
sions regarding the advantages-of various high-speed methods and 
their application are given. 


Chandler, R. F., Christian, R. A., Crash testing of humans in 
automobile seats, Proc. Society of Automotive Engineers Safety 
Conf. May 13-15, 1970, Detroit, Mich., and June 8-11, 1970, Brus- 
sels, Belgium, Intern. Safety Conf. Compendium No. P-30, 112-132 
(Society of Automotive Engineers, Inc., New York, N.Y., June 
1970). 


Key words: Automobile; automobile interior; automobile seats; 
crash testing. 


The design of the automobile interior is largely based upon design 
data gathered under static conditions, where the inertial properties 
of the human body do not enter into consideration. Anthropometric 
data based on static measurements are valid for design of clear- 
ance and access for normal operation of the vehicle, and have been 
widely collected and used for these purposes. The requirements of 
the human operator or passenger during a crash are not as well 
documented. Clearances, paths of travel, velocities along the path 
and other dynamic data are of vital interest to the designer. Be- 
cause of the difficulties in obtaining information about humans 
under crash conditions, most of the available ,uidelines are based 
upon dummy or cadaver tests. Until adequate human data are 
available, the designer must use the information obtained from 
human simulators and consider it only an unconfirmed approxi- 
mation. 

To alleviate this shortage of human impact data, the National 
Bureau of Standards initiated a research program in cooperation 
with the 657lst Aeromedical Research Laboratory. Tests were 
conducted on the Daisy Decelerator at Holloman AFB to gather 
comparative data about human and dummy test subjects under 
crash conditions using automotive seating and restraint systems. 
The data presented in this report include 32 human tests of con- 
trolled impacts of approximately the same magnitude. Variables 
are the type of restraint, either lap belt or lap belt plus single 
diagonal, and subject size. All tests were conducted in the —Gx 
orientation using a production automobile bucket seat. Data pre- 
sented include the displacement paths of the subjects during the 
impact, the maximum velocity along the path, anthropometric 
data, loads generated in the restraint system, and medical and 
subjective evaluation of the impacts. 


Christ, B. W., On the mechanism of interstitial-impurity- 
induced cross-slip in iron deformed near 175°K, Acta Met. 
17, No. 10, 1317-1321 (Oct. 1969). 


Key words: Cross slip; ductility; interstitial impurity; iron; nitro- 
gen; yield stress. 


A mechanism is proposed whereby the elastic dipole distortion of 
stationary interstitial impurity atoms in lattice solution can pro- 
mote nucleation and stabilization of double-kink-size dislocation 
segments on cross-slip planes. Stress relaxation accompanying 
interstitial-impurity-induced cross-slip can account for the experi- 
mental observations that small atom fractions of nitrogen in lattice 
solution in iron deformed around 175°K decrease yield and flow 
stresses and increase ductility. 


Christ, B. W., Gamble, R. P., Smith, G. V., On the distinction 
between alloy softening due to nitrogen and nickel in dilute 
lattice solution in iron, Script. Met. 3, No. 8, 521-530 (Aug. 
1969). 


Key words: Alloy softening; elastic; electronic; Fe-N alloys; 
Fe-Ni alloys; flow stress; lower yield stress; mobile dislocation 
density; Peierls stress. 


Alloy softening (AS) is the decrease in yield and flow stresses of 
iron and other b.c.c. metals which accompanies increasing quanti- 
ties of solute in lattice solution, within critical ranges of solute 
concentration and temperature. This note is a comparison of exist- 
ing alloy softening data on Fe-N and Fe-Ni. Whereas AS in each 
alloy system is due to a reduction in the thermal component of 
stress, the basic cause of the reduction is different in each system. 
This is inferred from differences in the concentration and tem- 


perature dependence of AS for each alloy system. AS in Fe-N is 
attributed to an increase in mobile dislocation density during 
yielding and flow, as compared with that of hydrogen-purified 
iron, arising from the mechanism of interstitial-impurity-induced 
cross-slip. AS in Fe-Ni is attributed to a systematic reduction in 
the Peierls stress of iron with increasing Ni. Whereas AS in Fe-N 
is primarily of elastic origin, AS in Fe-Ni is primarily of elec- 
tronic origin. 


Christensen, R. G., Hoeve, C. A. J., Comparison between theo- 
retical and experimental values of the volume changes ac- 
companying rubber extension, J. Polymer Sci. 8, Part A-l, 
1503-1512 (1970). 


Key words: Dilatometer; free energy; intramolecular interactions; 
intermolecular interactions; rubber elasticity; thermodynamics; 
volume change. 


The molecular theory of rubber elasticity assumes the free energy 
to consist of two parts: a liquid-like free energy that is governed 
by intermolecular interactions and is independent of strain at con- 
stant volume and an intramolecular interaction free energy equal 
to the sum of the free energies of the chains making up the net- 
work. The volume increases of rubber samples as a function of 
their length were found to be considerably larger than predicted 
by the molecular theory. Therefore, contrary to common belief, the 
values of (@E/AL)v.r may not be related solely to changes in 
intramolecular interactions with extension. Also, the usual proce- 
dure to obtain values of (@9E/@L)v,.r from measurements of 
(Of/0T)>».. with the aid of the molecular theory is not strictly 
correct. 


Clark, A. F., Childs, G. E., Wallace, G. H., Low-temperature 
electrical resistivity of some engineering alloys, (Proc. 1969 
Cryogenic Engineering Conf., June 16-18, 1969, Univ. of Cali- 
fornia, Los Angeles, Calif.), Chapter in Advances in Cryogenic 
Engineering 15, Paper No. C—5, 85-90 (Plenum Press, Inc., New 
York, N.Y. 1970). 


Key words: Aluminum alloys; copper alloys; electrical resistivity ; 
iron alloys; nickel alloys; titanium alloys. 


The electrical resistivities of several engineering alloys have been 
measured at 4K, 20K, 76K, 192K, and 273K. Different iron, nickel, 
copper, titanium and aluminum alloys have been measured, some of 
them in different heat treated conditions. Most of the alloys have 
a smooth decrease in resistivity with decreasing temperatures, but 
some unusual results are reported. Comparisons are made between 
different alloys and between different heat treatments of the 
same alloys. 


Currie, L. A., Indirect estimation of component variability 
in chemical and physical systems, Anal. Chem. 41, No. 14, 
2051-2054 (Dec. 1969). 


Key words: Chemical microstandards; correlation coefficient; in- 
direct variability estimation; stoichiometry; two-component system. 


An indirect method is presented for assessing the variability of a 
chemical component or a physical property, when that property 
may not be directly observed. The method is based upon a simple 
physical model incorporating stoichiometry and the observation of 
the variability of a second, related property and/or the ratio of 
the two properties. A principal result of the investigation is .an 
expression giving limits for the relative standard deviation of 
component-A, ga: | ¢74 2» — o's | = o's = G74 2 + O's where 
ga » and = represent the relative standard deviations of the ratio 
of components —A and —B and of component-B, respectively. 
Information is also obtained on bounds for the correlation co- 
efficient. 


Deslattes, R. D., Optical and x-ray interferometry of a silicon 
lattice spacing, Appl. Phys. Letters 15, No. 11, 386-388 (Dec. 
1969). 

Key words: Interferometer; lattice parameter; x-ray. 


A device permitting simultaneous x-ray and optical interferometry 
over traverses in excess of 20 um is reported. Results obtained to 











date suggest that such devices will permit measurements of certain 
crystal lattice parameters with accuracies better than one part 
per million. 


Durst, R. A., Fluoride microelectrode—fabrication and char- 
acteristics, Anal. Chem. 41, No. 14, 2089-2090 (Dec. 1969). 


Key words: Fluoride analysis; fluoride microanalysis; fluoride 
microelectrode; ion-selective microelectrode; microelectrode; 
potentiometry. 


A fluoride ion-selective microelectrode has been constructed with 
a volume requirement of less than 2 ul. This microelectrode has 
the advantage of being suitable for in situ, in vivo, and continuous 
analysis of microliter volumes of solution. Details of construction 
and the response characteristics are given. Nernstian response is 
observed down to less than 10°M fluoride. Further miniaturization 
is feasible by extension and refinement of this design. 


Ederer, D. L., Computer analysis of resonance profiles by the 
og - least squares, Appl. Opt. 8, No. 11, 2315-2325 
ov. F 


Key words: Beutler-Fano resonance profiles; least squares analysis; 
parameter correlation; random errors; slit function; slit width- 
resonance ratio; systematic errors. 


The method of least squares has been applied to the analysis of 
resonance profiles described by the Beutler-Fano (BF) absorption 
cross section. To illustrate the method, parameters were determined 
for experimental data having BF profiles of (a) the asymmetric, 
enhanced absorption type, (b) the “window” type, and (c) the 
Lorentzian absorption type. Systematic and random errors, and 
correlation among the parameters were studied as a function of the 
slit width to resonance width ratio by applying the least squares 
fitting method to data computed by folding the BF transmission 
profile with a gaussian slit function. High correlation among the 
parameters and parameter fractional standard deviations of 10% to 
20% for data with a standard deviation of 0.01 are an incentive to 
keep the slit width to resonance width ratio less than two for 
“window” type resonances. In addition it was found that the 
parameters describing a resonance could be determined with the 
greatest precision in a given time when the slit width to resonance 
width ratio was of order one. It was found that systematic errors 
greater than about one percent resulted if the number of mesh 
points used to approximate the convolution integral in the fitting 
procedure was less than at least two per resonance width. Also, for 
resonances whose width was equal to the slit width, systematic 
errors in the parameter values of order ten percent resulted if the 
width of slit function used in the fitting procedure differed from 
the “actual” slit function by ten percent. 


Eick, J. D., Caul, H. J., Hegdahl, T., Dickson, G., Chemical 
composition of dental gold casting alloy and dental wrought 
= alloys, J. Dental Res. 48, No. 6, 1284-1289 (Nov.-Dec. 


Key words: Casting gold alloys; chemical composition; dental gold 
alloys; gold; platinum; wrought gold wire; x-ray emission. 


The chemical composition of 136 dental casting gold alloys and 
21 dental wrought gold wire alloys was determined by an x-ray 
emission method. The changes in composition from 1928 and 1932 
to the present time are small. There has been a decrease in plati- 
num content in both castings and wires. About 1% nickel is found 
in the present-day wire which was not present in 1928 and 1932. 


Freeman, D. H., Currie, L. A., Kuehner, E. C., Dixon, H. D., 
Paulson, R. A., Development and characterization of ion- 
exchang: microstandards, Anal. Chem. 42, No. 2, 


e 
203-209 (Feb. 1970). 


Key words: Ion exchange; microstandards; nanogram. 


The encapsulation of Na or Ca by ion exchange beads is investi- 
gated in the counterion mass range below one nanogram. Prepar- 
ative methods are described for obtaining the exchanger in a 
pure counterion form and for subsequent prevention of ion ex- 
change so that the isolation of a single bead is possible in the 
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absence of significant contamination. Activation Analysis indi- 
cated that the network heterogeneity was characterized by a 5% 
relative standard deviation for Na in measurements extending to 
3 X 10°% grams. A bead diameter measurement is equivalent to a 
counterion mass determination and this opens the door to the use 
of the beads as microstandards in the mass region where no other 
methods are applicable. 


Freeman, D. H., Zielinski, W. L., Jr., Editors, Separation and 
Purification Section: Summary of activities July 1968 to 
June 1969, Nat. Bur. Stand. (U.S.), Tech. Note 509, 73 pages 
(Feb. 1970), 70 cents, SD Catalog No. C13.46:509. 


Key words: Gas chromatography; infrared analysis; ion exchange; 
liquid chromatography; purification; separation; styrene/divinyl- 
benzene; ultra-pure reagents. 


This is the annual progress report on the Separation and Purifica- 
tion Section activities. The major task of developing a certified ion 
exchange microstandard as a Standard Reference Material is pre- 
sented in terms of the preparative and characterizational work 
involved. Fundamental studies of ion exchange substrates includes 
the application of quantitative analytical infrared spectrophotom- 
etry to measure crosslinking in the copolymer network, and to 
determine the degree of sulfonation. Optical microscopy is applied 
to the study of swelling kinetics for single copolymer particles. 
Analytical gas chromatography is applied to the isomers of divinyl- 
benzene. The recently acquired quadrupole mass spectrometer is 
decribed. Liquid chromatography at high pressures and with high 
resolution has been begun recently. The activities of a project 
dealing exclusively with ultrapure reagents is described including 
the problems of contamination free storage. 


Gadzuk, J. W., Resonance tunneling through impurity states 
in metal-insulator-metal junctions, J. Appl. Phys. 41, No. 1, 
286-291 (Jan. 1970). 


Key words: Electron tunneling; junctions; impurities; metal- 
insulator-metal structures. 


Motivated by the recent work of Parker and Mead on the effects of 
impurity trapping states in Schottky barriers, a theory of resonance 
tunneling through absorbed atoms in field emission is modified to 
treat resonance tunneling through isolated and uniformly dis- 
tributed impurities in the insulating layer of a metal-insulator- 
metal tunnel junction. In analogy with the magnetic impurity 
tunneling theory of Appelbaum, the effects of resonant tunneling 
through impurity states opens up new tunneling channels with a 
concomitant change of slope in the I-V characteristics. These 
changes should be apparent as definite structure in curves of 
da/dV versus voltage although depending upon the impurity con- 
centrations and impurity level positions relative to the Fermi levels 
of the two metals, the height of the structure reflecting the im- 
purity energy level spectrum may not rise above the noise level. 
Implications of the theory are discussed and experimental sug- 
gestions made. 


Green, J. A. S., Mengelberg, H. D., Yolken, H. T., Oxide growth 
on copper and alpha-brasses in aqueous ammonia, 
J. Electrochem. Soc. 117, No. 4, 433-437 (Apr. 1970). 


Key words: Alpha-brass; copper; ellipsometry; oxide-growth; 
season-cracking; stress-corrosion. 


An ellipsometric study of the kinetics of oxide growth for pure 
copper and a series of alpha-brasses immersed in tarnishing am- 
moniacal environments has established that the growth of the oxide, 
shown to be cuprous oxide, obeys a linear law in each case, and 
that the rate of growth increases significantly with (a) increasing 
zinc content of the solid, (b) increasing temperature and (c) the 
application of anodic potentials. On the basis of these results it is 
suggested that the oxide is porous and that oxide growth involves 
the dissolution of metal ions at the base of the pore (anodic 
reaction) and their deposition at the oxide surface (cathodic 
reaction), electrons flowing through the oxide film. The influence 
of zinc content on the oxide-growth rate is attributed to its effect 
on the dissolution kinetics. The relevance of these observations to 
season-cracking is considered. 








Haller, W., Blackburn, D. H., Wagstaff, F. E., Charles, R. J., 
Metastable immiscibility surface in the system Na,O-B,0-— 
SiO., J. Am. Ceram. Soc. 53, No. 1, 34-39 (1970). 


Key words: Glass; immiscibility; melts; microheterogeneities; 


Na:O-B:0;-Si0,; B:0;-NasO-Si0s, Si0--B,O.—Na,0. 


Opalescence and clearing techniques have been used to determine 
the metastable immiscibility surface for sodium borosilicate solu- 
tions. These results indicate that a “three-liquid region, which may 
or may not be metastable to two-liquid regions, underlies the 
immiscibility surface. 


Hamer, W. J., DeWane, H. J., Electrolytic conductance and the 
conductances of the halogen acids in water, Nat. Stand. Ref. 
Data Ser., Nat. Bur. Stand. (U.S.), 33, 37 pages (May 1970) 50 
cents, SD Catalog No. C13.48:33. 


Key words: Conductances of HF, HCl, HBr, and HI; electrolytic 
conductance; theories of electrolytic conductance. 


Definitions, symbols, general principles, and general laws related 
to the electrolytic conductance of aqueous solutions are presented. 
The general laws considered are Coulomb’s law for charged bodies, 
Poisson’s equation relating the electrostatic potential to charge 
distribution, and the Stokes and Oseen laws for the velocity of a 
sphere in a fluid medium. The relations between electrical re- 
sistance, electrical conductance, specific resistance, specific con- 
ductance, and equivalent conductance are set forth. Theoretical 
expressions for the equivalent conductance as derived by Debye, 
Onsager, and Fuess are given in general form and in a somewhat 
more detailed fashion in an appendix. The general methods of 
treating the equivalent conductances of ionophores and ionogens, 
especially in regard to the determination of the limiting equivalent 
conductance, the degree of ionic association, and the degree of 
ionic dissociation are discussed. Data on the equivalent conduct- 
ances of the halogen acids, hydrofluoric, hydrochloric, hydrebromic, 
and hydriodic acids in water are given for a wide range of 
concentration and temperature. 


Heinrich, K. F. J., Advances in the metallurgical application 
of electron probe microanalysis, Proc. Conf. 5th International 
Congress on X-ray Optics and Microanalysis, Tubingen, Germany, 
Sept. 9-14, 1968, G. Mollebstedt and K. H. Gaukler, Eds., pp. 415- 
423 (Springer-Verlag Publ. Co., Berlin, Germany, Dec. 1969). 


Key words: Applications; electron probe; microanalysis; metal- 
lurgy; x-rays; spectrometry. 


Electron probe microanalysis is widely used in the investigation 
of phenomena of general metallurgical interest as well as in the 
study of specific materials and processes. New devices for x-ray 
detection, and automatic data collection, increases the usefulness 
of the microprobe for such applications. Besides the well-known 
methods for qualitative and quantitative electron probe analysis, 
the instrument can also be used advantageously for quantitative 
metallography and related techniques. 


Heinrich, K. F. J., Fiori, Cc Yakowitz, H., a 
technique for scanning electron mi y and electron 
probe microanalysis, Science 167, 1129-1131 ( tr eb. 20, 1970). 


Key words: Derivative signal; electron probe; image; scanning 
electron microscope; secondary electrons; target current. 


A technique is described for producing improved topographic 
images on the scanning electron microscope and the scanning 
electron probe microanalyzer. In this technique, the brightness of 
the oscilloscope is modulated by a signal obtained by mixing the 
signal (from secondary electrons or target current) with its first 
derivative. This enhances minor topographic features which are 
poorly reproduced in the conventional technique. 


Heinrich, K. F. J., Yakowitz, H., Propagation of errors in 
correction dels for q titative electron probe micro- 
analysis, Proc. Conf. 5th International Congress on X-ray Optics 
and Microanalysis, Tubingen, Germany, Sept. 9-14, 1968, G. Mol- 
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lebstedt and K. H. Gaukler, Eds., pp. 151-159 (Springer-Verlag 
Publ. Co., Berlin, Germany, Dec. 1969) 


Key words: 
effects; 
analysis. 


Absorption of x-rays fluorescence; 
error propagation; 


atomic number 
microprobe analysis; quantitative 


Lack of appropriate standards frequently forces the analysis to use 
elemental standards. The usefulness of correction models is limited 
by the accuracy to which the input parameters are known. Un- 
certainties in presumably known quantities (mass absorption, co- 
efficients, fluorescence yield, mean ionization potentials, etc.) are 
in many cases the limiting factors. The resulting analytical errors 
can be minimized by judicious choice of experimental conditons. 
This paper will give examples involving the corrections for absorp- 
tion, fluorescence by characteristic lines, and atomic number 
effects. 


Heydemann, P. L., Ultrasonic measurements at pressures up 
to 50 kbar, Proc. Intern, Collog. Solids at High Pressures, 
Grenoble, France, Sept. 8, 1969, No. 188, pp. 461-467 (Centre 
National de la Recherche Scientifique, Paris, France, 1970). 


Key words: Bulk modulus; equation of state; 


high pressure; 
phase transitions; potassium chloride; 


tellurium; ultrasonics. 


Ultrasonic measurements as a function of pressure or temperature 
provide one of the most accurate means to determine the pressure 
or temperature dependence of the bulk modulus and the density 
of both liquids and solids. Such measurements are currently car- 
ried out with very high accuracy at pressures up to about 4 kbar 
(1 kbar = 10 Nm’) in several laboratories. 

If the material under investigation undergoes transitions at higher 
pressures, or if the bulk modulus is noticeably non-linear with 
pressure, ultrasonic measurements at pressures higher than 4 kbar 
are needed. This report describes our methods to carry out such 
measurements at pressures close to 50 kbar, both under hydro- 
static and non-hydrostatic conditions. 

For most materials the accuracy of the density and bulk modulus 
data obtained from such measurements is much higher than that 
obtained from isothermal dilatometric measurements. In the case 
of tellurium, the accuracies for the bulk modulus at 30 kbar are 
1.6 and 23%, respectively. 

Another interesting application of ultrasonics at high pressures is 
the detection of phase transitions and the measurement of the 
elastic properties as the material undergoes the transition. As an 
example for measurements through a transition range, our results 
of ultrasonic and isothermal dilatometric measurements on KC] are 
presented. We find that the ratio of specific heats also shows a 
sharp increase in the transition range. 


Hust, J. G., A compilation and historical review of tempera- 
ture scale differences, Cryogenics 9, No. 6, 443-455 (Dec. 1969). 


Key words: Temperature; temperature scales; thermometry. 


A brief review is given of temperature, its definition, and measure- 
ment. Methods of practical temperature measurement are de- 
scribed. Differences between internationally accepted temperature 
scales are presented graphically and in tabular form; comparisons 
are also given for several national temperature scales. Emphasis 
has been placed on temperatures below 0°C, but some informa- 
tion on temperatures above 0°C has also been included. Suggested 
methods of conversion between different temperature scales are 


described. 


Ireland, C. T., Ku, H. H. Kullback, S., Symmetry and marginal 
homogeneity of an rXr contingency table, J. Am. Statist. 
Assoc. 64, 1323-1341 (Dec. 1969). 


Key words: Contingency table analysis; homogeneity of margi- 
nals; minimum discrimination information estimation; test for 
symmetry. 


The principle of minimum discrimination information estimation is 
employed to obtain RBAN (Regular Best Asymptotically Normal) 
estimates of the cell frequencies of an r Xr contingency table under 
hypotheses of symmetry and marginal homogeneity. For the latter 

















a convergent iterative procedure is given to compute the estimates. 
The associated minimum discrimination information statistics are 
distributed asymptotically as x* under the null hypothesis. The 
procedures differ from those previously presented in the literature 
and permit of extension to multidimensional contingency tables. 
An example is given. 


Jones, M. C., Giarratano, P. J., Simpson, A. U., Heat transfer to 
solid-vapor mixtures of cryogens below their triple points 
ry through heated tubes, AJChE J. 15, No. 5, 890-897 
(Nov. 1969). 


Key words: Boundary layer; cryogenic; heat transfer; hydrogen; 
low pressure; nitrogen; particle-wall interaction; solid-vapor; 
tube; two-phase. 


Data are presented for wall temperatures and heat transfer co- 
efficients for solid-vapor mixtures of parahydrogen and nitrogen 
flowing in an electrically heated straight tube of length 40 times 
its diameter. These are interpreted by the application of flat plate, 
constant property boundary layer theory to models in which the 
solid particle geometrical distribution takes on simple limiting 
forms. The observed enhancement of the heat transfer coefficient 
over that for gas alone traveling at the same velocity is qualita- 
tively predicted as a function of a dimensionless heat flux, the 
Sterman parameter gw/p»UX. This parameter is also shown to be 
an index of particle-wall interaction and is used to correlate ob- 
served limiting wall heat fluxes above which free flow may be 
guaranteed. 


Kamper, R. A., The Josephson effect, JEEE Trans. Electron 
Devices ED-16, No. 10, 840-844 (Oct. 1969). 


Key words: Electron tunneling; electronics; low temperature; 
superconducting devices. 


This is a review of the Josephson effect in superconductors with 
emphasis on the electrical properties of Josephson junctions and 
their application to devices and measurement techniques. 


Kilgore, S. B., Electronic circuit breaker for motor load pro- 
tection, Rev. Sci. Instr. 41, No. 1, 40-41 (Jan. 1970). 


Key words: Electronic circuit breaker; mechanical load protec- 
tion; surge current; time delay. 


An electronic circuit breaker, designed to switch off the power to 
an electric motor when the mechanical load increases, is described. 
The circuit operates rapidly and can sense an increase in mechani- 
cal load before it reaches destructive proportions. The sensitivity 
of the electronic circuit makes it necessary to incorporate a time 
delay that prevents the initial surge current from turning the 
motor off before it gets started. 


Kostkowski, H. J., Erminy, D. C., Hattenburg, A. T., High- 
accuracy spectral radiance calibration of tungsten-strip 
lamps, Chapter 4 in Advances in Geophysics 14, 111-127 (Aca- 
demic Press Inc., New York, N.Y. 1970). 


Key words: Blackbody; high temperature measurement; radiom- 
etry; spectral radiance; standards calibration; tungsten strip lamp. 


The accuracy of calibrating tungsten strip lamps in terms of abso- 
lute spectral radiance has been significantly improved. The stand- 
ard deviation uncertainty of these calibrations is estimated to vary 
from 0.3 percent at 850 nm to 1.5 percent at 210 nm. This accu- 
racy has been realized as a result of: (1) The development of a 
highly stable 3000 °K blackbody and a detailed inevstigation of its 
quality. (2) Increasing the accuracy with which the blackbody 
temperature could be determined. (3) Reducing and evaluating 
errors resulting from scattered light, polarization, and slit-function 
effects. Highlights of the instrumentation, the temperature meas- 
urement and the various error investigations are presented. 


Ku, H. H., Varner, R. N., Kullback, S., Analysis of multi-dimen- 
sional contingency tables, (Proc. 14th Conf. Design of Experi- 
ments in Army Research Development and Testing, Edgewood 
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Arsenal, Md., October 23, 1968), ARO-D Report 69-2, 141-180, 
(U.S. Army Research Office-Durham, Durham, North Carolina, 
Sept. 1969). 


Key words: Computer programs; contingency tables; estimation 
of cell frequencies from marginals; generalized independence; 
higher-order interaction; hypothesis testing; information theory; 
interaction. 


This is an expository paper on the analysis of multi-dimensional 
contingency tables. The principle of minimum discrimination in- 
formation estimation is described and used to generate estimates 
for tests of hypotheses concerning various interactions and effects. 
All classical hypotheses for contingency tables can be generated by 
the use of this principle when certain marginals are considered 
as fixed. 

Examples are given and two available computation programs are 
described in detail. 


Kusuda, T., Algorithms for psychometric calculations (skele- 
ton tables for the thermodynamic properties of moist air), 
Nat. Bur. Stand. (U.S.), Bldg. Sci. Ser. 21, 53 pages (Jan. 1970), 
55 cents, SD Catalog No. C13.29/2:21. 


Key words: Computer algorithm; psychrometrics; saturated and 
unsaturated moist air; thermodynamic properties. 


Computer algorithms to obtain thermodynamic properties of satur- 
ated and unsaturated moist air are presented in this paper. The 
sautrated moist air properties are calculated by the methodology 
developed by J. A. Goff and S. Gratch for their ASHRAE tables 
(1967 Book of Fundamentals, The American Society of Heating, 
Refrigerating and Air Conditioning Engineers). Sample calcula- 
tions were performed using a computer program based upon the 
algorithms presented herein and the results are attached. 


Lang, S. B., Rice, L. H., Shaw, S. A., Pyroelectric effect in 
barium titanate ceramic, J. Appl. Phys. 40, No. 11, 4335-4340 
(Oct. 1969). 


Key words: Barium titanate; ceramics; ferroelectricity; phase 
transition; pyroelectricity; spontaneous. 


The pyroelectric coefficient (at constant stress) of a 95 percent 
BaTiOs, 5 percent CaTiOs ceramic composition (Clevite Ceramic 
B) was measured over the temperature range from 4.9 to 400°K. 
Extrema in the pyroelectric coefficient were observed at the rhom- 
bohedral-orthorhombic transition point (150 to 170°K), the ortho- 
rhombic-tetragonal transition point (240 to 260°K), and the 
tetragonal-cubic transition point (Curie point, at 387°K). The 
pyroelectric coefficient showed no anomalous behavior below 
150°K. A spontaneous polarization curve derived from the pyro- 
electric results is compared with published single-crystal measure- 
ments. A model for ceramic polarization suggests that domain 
switching at phase transitions is biased by electrical interactions 
arising from space charges and by mechanical interactions result- 
ing from internal strains. 


Latanision, R. M., On the dislocation distribution near the 
surface of lightly deformed copper single crystals, Scripta 
Met. 3, No. 7, 465-470 (July 1969). 


Key words: Debris layer; edge dislocation; near surface source; 
screw dislocation; single crystals; surface source. 


Two recent and essentially identical etch-pitting studies on oriented 
copper single crystals have produced contradictory results. An 
explanation for the discrepancy is not obvious, but the results of 
these investigations may be interpreted in terms of either the sur- 
face “debris” layer hypothesis or the activation of surface or near- 
surface dislocation sources. The purpose of this note is to empha- 
size the ambiguity involved in such studies. 


Lawton, R. A., New standard of electric field strength, JEEE 
Trans. Inst. Meas. IM-19, No. 1, 45-51 (Feb. 1970). 


Key words: Antenna; electric field strength; measurement; stand- 
ard field. 








A new technique is presented for establishing a standard of field 
strength using a conducting sphere. An analysis is made to deter- 
mine the current on the sphere as a function of the electric field 
strength of an incident plane wave, a method of measuring that 
current using electronic circuitry and an optical indicator within 
the sphere is described, and an intercomparison is made with an 
independent field-strength standard. 

This technique is a significant improvement over previous tech- 
niques in that it permits the absolute determination of field 
strength with a maximum uncertainty of 1 percent or less at 30 
MHz. This technique is applicable to a broad range of frequen- 
cies and field strengths. 


McCrackin, F. L., Analyses and corrections of instrumental 


errors in ellipsometry, J. Opt. Soc. Am. 60, No. 1, 57-63 (Jan. 
1970). 


Key words: Birefringence; ellipsometer; thin films. 


The errors produced in an ellipsometer due to an imperfect com- 
pensator, due to birefringence in windows of a cell containing the 
sample, and due to tilting of the sample are derived. Methods of 
analyzing the ellipsometer readings to eliminate these errors are 
given. 


McDonald, D. G., Kose, V. E., Evenson, K. M., Wells, J. S., Cupp, 
J. D., Harmonic generation and submillimeter wave mixing 
with the Josephson effect, Appi. Phys. Letters 15, No. 4, 121- 
122 ( Aug. 15, 1969). 


Key words: Infrared; josephson effect; lasers; submillimeter 
waves; superconductivity. 


By observing constant-voltage steps from Josephson junctions at 
voltages as high as 17 mV we deduce that junctions can generate 
harmonics up to frequencies as high as 8200 GHz. In consonance 
with this, submillimeter wave laser detection, harmonic generation 
and mixing are demonstrated. These results suggest a model for 
the upper frequency limit of the Josephson effect. 


McLaughlin, W. L., Hussmann, E. K., The measurement of 
electron and gamma-ray dose distributions in various media 
(Proc. Symp. Utilization of Large Radiation Sources and Accelera- 
tors in Industrial Processing, Munich, Germany, Aug. 18—22, 1969), 
Chapter in Large Radiation Sources for Industrial Processes, 
IAEA-SM-123/43, pp. 579-590 (International Atomic Energy 
Agency, Vienna, Austria, 1969). 


Key words: Depth dose; dose distributions; dosimetry; dyes; elec- 
tron beams; films; gamma rays; microdosimetry; radiographic 
imaging. 


New radiochromic dye films and gels are proposed as simple and 
relatively accurate devices for imaging and determining the energy 
deposited locally at different depths and across interfaces of vari- 
ous irradiated media. Depth-dose measurements show remarkable 
agreement with theoretical computations in various plastics and 
aluminum, even near the end of the electron range, indicating that 
the dosimeter response is representative of the energy imparted 
to the medium over a wide range of electron energies. 


Margoshes, M., Rasberry, S. D., Application of digital com- 
puters in spectrochemical analysis—computational methods 
in photographic microphot try, Spectrochim. Acta 24B, 
497-513 (1969). 





Key words: Computers; emulsion; microphotometry; photographic 
emulsion; photographic photometry; photometry; programs; spec- 
trochemical. 


A new method is described for calibration of photographic emul- 
sions and conversion of microphotometer readings to relative in- 
tensities on a digital computer with special application to spectro- 
chemical analysis. The method for emulsion calibration replaces 
graphical procedures by a numerical method which is well suited 
for digital computations. The method of computation and its under- 
lying assumptions are described, and results are given of tests of 
the procedure. It is recognized that the entire experimental 





arrangement is being calibrated, not merely the photographic 
emulsion, and it is shown how malfunctions of the microphotometer 
which affect the calibration can be recognized from the output. of 
the computer. The program for conversion of microphotometer 
readings to relative intensities provides for several alternate calcu- 
lations, including selection of the proper calibration parameters 
according to the wavelength of the line when these data are sup- 
plied for more than one wavelength region, as well as correction 
for step and background where required. The program is written 
to provide for automatic selection of the required alternative 
calculation, based on preliminary analysis of the input data. 


Meinke, W. W., The NBS standard reference materials pro- 
gram: Past, present, and future, Mater. Res. Std. 9, No. 10, 
15-18 (Oct. 1969). 


Key words: Primary standards; standard reference material; 
standards. 


The past 60 years have been ones of rapid and sweeping techno- 
logical change in the industrial world. During this time the Na- 
tional Bureau of Standards (NBS), has played a pivotal role in 
the orderly progress and quality control of the nation’s productive 
process. One of the means by which this progress has been aided 
is through the production, certification, and issuance of Standard 
Reference Materials (SRM) from NBS. 

The development of our high-technology society is based in large 
part on the mass production of goods and services. In turn, mass 
production implies strict quality control of the means of produc- 
tion, an ever-increasing accuracy of measurement, interchange- 
ability of parts, and stringent performance criteria for the wide 
variety of materials now used. The base of this structure rests 
firmly on a well-established standards program, national in scope, 
in which the NBS-SRM program plays a significant role. It is 
interesting to note that the national standards bodies, the ASTM 
being among the foremost, have interacted from the earliest days 
with the NBS-SRM program. 


Meyerson, M. R., Friedman, L., Giles, P. M., Dimensional sta- 
bility of a 12% nickel maraging steel at ambient tempera- 
tures, Trans. ASM 62, No. 3, 809-812 (Sept. 1969). 


Key words: Gage blocks; maraging steel; nitriding; precision 
equipment; surface hardened; temporal dimensional stability; 
ultrastability. 


The maraging steels attain high strength with a relatively simple 
low temperature heat treatment. The nature of the heat treatment 
and the low carbon content of the steels indicated low residual 
stress levels, small lattice distortions, nitridability and other char- 
acteristics which are known to promote dimensional stability at 
high strength levels. A 12% nickel maraging steel was treated to 
produce a surface hardness of greater than R- 65 and its dimen- 
sional stability at submicroinch levels was observed over a period 
of about 2% years. The steel as treated was found to be ultrastable, 
with an average decrease in length of only 0.05 to 0.10 micrometers 
per meter per year. Its performance meets or exceeds that of any 
of the many other materials previously investigated at NBS. 


Mosburg, E. R., Jr.. A new probe configuration for plasma 
diagnostics, J. Appl. Phys. 40, No. 13, 5290-5300 (Dec. 1969). 


Key words: Plasma diagnostics; plasma physics; probes. 


A new probe geometry using a large array of small, closely spaced 
electrodes, alternately connected to two different potentials has 
been developed. The basic characteristics of the current-voltage 
curve of this probe differ markedly from that of the usual double 
probe. 


Mullen, L. O., On efficient use of selective vacuum pumping 
modes, J. Environ. Sci. 12, No. 5, 26-30 (Oct. 1969). 


Key words: Cryopumps; getters; ion pumps; pumping efficiency; 
sublimation pumps; vacuum pumps; vapor pressure. 


Most vacuum pumps are selective in the sense that they pump 
different molecular species at different speeds. Some pumps, such 





as sublimation pumps and cryopumps, are extremely selective. 
Efficient use of selective pumps requires that they be chosen to 
match the gas load. We present experimental data on some com- 
monly encountered gas loads and discuss the choice of compatible 
pumps. 


Ohashi, M., Paffenbarger, G. C., Some flow characteristics at 
37°C of ternary wax mixtures that may have possible dental 
uses, J. Nihon University School Dentistry 11, No. 3, 109-208 
(Sept. 1969). 


Key words: Dental; flow; wax. 


One hundred and eighty two ternary diagrams of the flow at 37°C 
of mixtures of commercial waxes indicate some of these mixtures 
may have dental use. 


Paulsen, P. J., Alvarez, R., Kelleher, D. E., Determination of 
trace elements in zinc by isotope dilution spark source mass 
spectrometry, Spectrochim. Acta 24B, 535-544 (April 1969). 


Key words: Electrodeposition; isotope dilution; preconcentration; 
spark source mass spectrometry; trace elements; zinc. 


The applicability of the isotope dilution technique in spark source 
mass spectrometry has been investigated for the determination of 
Pb, Cu, Cd, Ag, Tl, and Sn in zinc metal. In the procedure de- 
veloped, a solution of the zinc is spiked with known amounts of 
stable isotopes of the element to be determined, and the elements 
are electrodeposited onto gold wires that subsequently serve as 
electrodes for the mass spectrographic isotopic analysis. Chemical 
operations, including isothermal distillation of reagents, dissolu- 
tion of the zinc, and the electrodeposition step are conducted in a 
closed system to minimize contamination. The results of the analy- 
sis of two NBS zinc standard reference materials show determina- 
tion ranging from 11 ppm (11 ug/g) for Pb to as low as 0.02 ppm 
(20 ng/g) for Sn. It is concluded that this method is of general 
utility for simultaneous multi-element trace analysis and is espe- 
cially applicable in providing much needed standardization in 
spark source mass spectrometry. 


Pizer, R. S., Spinner, S., An apparatus for determining the 
actual footprint area of tires, Mater. Res. Std. 10, No. 2, 20-23 
(Feb. 1970). 


Key words: Contact area; footprint; optical and electrical sys- 
tem; silicon solar cell; tire; tires. 


An apparatus is described for measuring the actual contact area of 
tires regardless of their tread design. The procedure is both rapid 
and accurate. The method utilizes a bank of lights, an opal glass 
plate, a simple light focusing system, and a silicon solar cell which 
linearly transforms the incident light to an electrical current. 
The reduction in current is the measure of the true contact area 
of the tire. 


Potzick, J., Low-frequency sine-wave oscillator, EEE 18, No. 
3, 130 (Mactier Publishing Corp., Mar. 1970). 


Key words: Low frequency oscillator; oscillator. 


A simple, inexpensive, low-frequency oscillator which generates a 
fairly low distortion sine wave at a frequency of about one hertz 
is described. 


Randall, E. B., Jr., Lashof, T. W., An interlaboratory study of 
the measurement of the bursting strength of paper, Tappi 
53, No. 5, 799-809 (May 1970). 


Key words: Bursting strength; evaluation; interlaboratory study; 
paper; precision; reference materials; standards; test method. 


This paper reports on the results of an extensive interlaboratory 
study, involving more than 175 laboratories in North America and 
overseas using several different models of bursting strength testers 
to evaluate 34 papers and 6 foils. Significant differences, both in a 
statistical and in an economic sense, were found between results 
obtained on different models and in different laboratories using 
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the same model, confirming previous reports. There were also 
significant day-to-day variations within a laboratory, and aluminum 
foil specimens gave results relative to those of paper specimens 
that depended on the instrument model used. A convenient rule 
for the number of replications was obtained from the study; 
namely, make 10 replications unless the range of the first 10 is 
greater than 20% of their average, then make an additional 10. 
Two test results, each with 10 or 20 replications as above, are 
expected to agree with each other within the following values 95% 
of the time: 5.4% for repeatability within a laboratory, 9.5% for 
comparability between materials, and 14.3% for reproducibility 
between laboratories. Between laboratory comparisons may be 
standardized through the use of a composite reference material of 
paper, aluminum foils being of only marginal use for this 
standardization. 


Reed, R. P., Durcholz, R. L., Cryostat and strain measurement 
for tensile tests to 1.5°K, (Proc. 1969 Cryogenic Engineering 
Conf., June 16—18, 1969, Univ. of California, Los Angeles, Calif.), 
Chapter in Advances in Cryogenic Engineering 15, Paper No. D-3, 
109-116 (Plenum Press, Inc., New York, N.Y. 1970). 


Key words: Capacitance extensometer; cryostat; strain measure- 
ment instruments; tensile properties. 


A new tensile cryostat has been developed and tested at tempera- 
tures between 300 and 1.5K. Better design, insulation, and vacuum- 
sealed environment reduces the heat leak to about 0.2 watts dur- 
ing continuous pumping above liquid helium. New capacitance 
strain-sensing equipment has been developed to provide ranges of 
strain sensitivities between 10~° and 10~* inches per inch. Details 
are presented of both devices. 


Robbins, R. F., Cryogenic properties of a polyurethane ad- 
hesive, J. Macro-mol. Sci. A3, No. 7, 1367-1380 (Nov. 1969). 


Key words: Adhesive; low temperature; mechanical properties; 
polymer; polyurethane; resilience; thermal properties. 


Differential Thermal Analysis (DTA), Rebound Resilience, and 
fensile Properties of a polyurethane adhesive were measured at 
cryogenic temperatures. The experimental methods are described, 
and test results which aid in evaluating the polyurethane for use 
at low temperatures are discussed. The DTA thermogram reveals 
that the glass transition temperature (T,) is 235°K. The resilience 
profile indicates a resilience minimum (7m) at 270°K and a fre- 
quency of 3800 Hz, which is consistent with the T, measured by 
DTA. The low resilience below Tm, caused by secondary low tem- 
perature transitions, shows the high energy absorption capabilities 
of the polyurethane. Considerable plastic flow at 195°K (40°K 
below T,) is evidenced in the results of the tensile tests. The re- 
sults of the three tests indicate that the polyurethane adhesive will 
perform well at low temperatures. The test methods should also 
be useful for evaluating the low temperature performance of new 
experimental polymers. 


Robinson, H. E., Editor, Durability of insulating glass, Pro- 
ceedings of a seminar held at the National Bureau of Standards, 
Gaithersburg, Maryland, November 14-15, 1968, Nat. Bur. Stand. 
(U.S.), Bldg. Sci. Ser. 20, 84 pages (Feb. 1970), 75 cents, SD 
Catalog No. C13.29/2:20. 


Key words: Accelerated laboratory tests; double-glazed window 
units; factory-sealed insulating glass units; field performance tests, 
correlation with laboratory tests; sealant performance; standard- 
ized testing; test methods. 


A two-day seminar on the Durability of Insulating Glass was at- 
tended by some 130 persons on November 14 and 15, 1968. The 
seminar was held at the Gaithersburg, Maryland, facilities of the 
National Bureau of Standards and featured fourteen speakers who 
participated in panel discussions or delivered individual papers. 
Numerous agencies interested in design, manufacture, specifica- 
tion, purchase, installation or maintenance of windows were 
represented at the seminar. 

Among the topics considered in the panel discussions were: (1) 
The need for reliability and durability of insulating glass; (2) 
manufacturers’ test methods; (3) proposals for future action. The 








Canadian experience with an accelerated test method and accept- 
ance program was presented and discussed, as were the Norwegian 
accelerated test methods and their correlation with field experi- 
ence. A review of current practices leading to new test methods 
and standards was also presented, and a “round robin” program 
that would compare various test methods now employed in the 
industry was proposed. Affirmative interest in participating in the 
proposed round robin was expressed by about a score of manu- 
facturers present at the seminar. 


Silverman, S., The National Measurement System—a concept 


to assist the private sector, Mater. Res. Std. 9, No. 10, 11-14 
(Oct. 1969). 


Key words: Measurement; National Bureau of Standards; Na- 
tional Measurement System; national standards; self-calibration; 
standards, national; systems concept. 


The concept of : National Measurement System, involving the 
entire complex of measurement activities within the United States, 
is developed. The central element of this Measurement System is 
the National Bureau of Standards. The role of NBS and of other 


functional elements of the System, such as private laboratories, 
is discussed. 


Simpson, A. U., Timmerhaus, K. D., Kreith, F., Jones, M. C., Heat 
and mass transfer in dispersed, two-phase, single-com 
nent flow, Intern. J. Heat Mass Transfer 12, No. 9, 1141-1155 
(Sept. 1969). 


Key words: Boundary layer; computation; dispersed flow; heat 
transfer; solid-vapor; two-phase. 


The process of heat transfer to a two-phase mixture of well- 
dispersed subliming particles and vapor, flowing over a heated 
surface, is analyzed. It is shown by a laminary boundary analysis 
that, when the surface-area per unit volume of the particle (or 
solid) phase is large enough, the phase change dominates the heat 
transfer process and hastens the development of the thermal 
boundary layer. Under these conditions, the thermal boundary- 
layer thickness not only becomes uniform a short distance down- 
stream from the starting point, but also is substantially less than 
it would be were the particle phase absent. For such systems, the 
equations describing the heat transfer process can be considerably 
simplified and, if the physical properties of both phases are uni- 
form, a remarkably simple solution results. For systems in which 
the physical properties are not uniform, a solution involving inte- 
gration across the boundary layer is developed. The solutions are 
applicable to developing, as well as fully developed, laminary 

undary layers over a flat plate; the solutions also approximate 
conditions in flow through a tube, provided that the tube radius 
is large compared to the thermal boundary-layer thickness. The 
predictions of this theoretical analysis agree satisfactorily with 
experimental results. With slight modification, the same approach 


— possibly be applicable for turbulent flow in the boundary 
ayer. 


Simpson, J. A., Electron impact spectroscopy, (Proc. 10th 
Symp. on Electron, lon & Laser Beam Technology, Gaithersburg, 
Md., May 21-23, 1969), Chapter in Record of 10th Symposium, 
Electron, Ion, and Laser Beam Technology, L. Marton, Ed., pp. 
345-352 (San Francisco Press, San Francisco, Calif., 1969). 


Key words: Chemical analysis; electron impact spectrometer; uv 
absorption. 


Recent developments in electron spectrometer design have made 
the application of inelastic electron scattering measurements to 
gas analysis competitive with other techniques. The energy distri- 
bution of electrons of an initially monoenergetic electron beam 
after an encounter with a gas target, the loss spectrum, contains 
the optical absorption spectrum of the gas. This absorption data 
extends from the extreme vacuum ultraviolet into the visible with 
an energy resolution comparable to a % meter optical grating 
spectrometer. Its response is linear with concentration over a very 
wide range and contains between 10° and 10‘ resolution widths. 
This spectrum by revealing the valence energy states of the gas is 
the most intrinsic possible “fingerprint” of the atom or molecule. 
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An instrument has been built to explore the potentialities of this 
method. The instrument described and its performance as a trace 
analyzer for air pollution studies discussed. 


Smith, R. V., The influence of surface characteristics on the 


boiling of cryogenic fluids, J. Eng. Ind. 91, No. 4, 1217-1221 
(Nov. 1969). 


Key words: Boiling; boiling hysteresis; boiling inception; boiling- 
site spreading; radiation. 


This paper reviews the influence of the characteristics of the solid 
surface on the boiling of cryogenic fluids. Particular emphasis is 
placed on papers presented in the cryogenics section of the 1968 
Annual ASME meeting. 

Hysteresis, boiling-site spreading and the influence of radiation on 
the surfaces are discussed in detail. 


Snelleman, W., Menis, O., Rains, T. C., Yee, K. W., Cook, H. D., 
Flame emission spectrometry with repetitive optical scan- 
ning in the derivative mode, Anal. Chem. 42, No. 3, 394-398 
(Mar. 1970). 


Key words: A-c scan; derivative mode; design; detection limits; 
flame emission spectrometry; interferences; matrix; microsamples. 


A flame emission spectrometer using a rapid repetitive scan of a 
narrow wavelength region has been developed. By this method of 
wavelength scanning the second derivative of the output intensity 
is measured. This approach to matrix problems and micro samples 
greatly enhances the potentialities of flame emission spectrometry, 
and minimizes the need for a monochromator of high power. A 
quartz plate, which is vibrating at 145 Hz, is mounted in back of 
the entrance slit of the monochromator, and a voltage at its fre- 
quency of oscillation is used as the reference signal for a lock-in 
amplifier. When the ac amplifier is tuned to twice the frequency 
of vibration the second derivative of the intensity is obtained. 
This permits the measurement of weak line spectra nested in or 
on a broad band or continuum. It is demonstrated that spectral 
interference due to CaOH bands and/or a continuum are mini- 
mized in the measurement of barium, thus leading to an increase 
in detection limits of alkali and alkali earth elements in the 
presence of many matrix ions. An analysis can be performed on 
50 ul of sample which makes it applicable to biochemical and air 
pollution studies. 


Strobridge, T. R., Refrigeration techniques, U.S. Bureau Mines 
Info. Circ. No. 8417, pp. 39-56 (1969). 


Key words: Applications; cryogenic; refrigeration; thermodynamic 
cycles. 


The supply of helium accumulated through the conservation pro- 
gram has allowed technology in many disciplines to freely make 
use of the unique behavior of helium. Helium is the only fluid that 
can be conveniently used as a refrigerant below about 14°K and 
is often used at higher temperatures. The thermodynamic cycles 
that modern helium refrigerators utilize are reviewed. The Joule- 
Thomson process is important because it forms the lowest tem- 
perature stage of virtually all helium temperature refrigerators. 
The various processes and cycles for precooling the helium stream 
prior to the Joule-Thomson process form the essential differences 
between the closed cycle refrigerators operating today. Typical 
refrigeration and liquefaction facilities, both large and small, are 
discussed. 


Strohlein, J. N., Bagg, T. C., A selective roll-to-roll printer for 
producing duplicate microfilm copies, Nat. Bur. Stand. (U.S.), 
Tech. Note 516, 9 pages (Feb. 1970), 25 cents, SD Catalog No. 
C13.46:516. 


Key words: Information retrieval device; microcopier; microfilm 
duplicator; reader-printer; selective copier. 


The “Selectaframe” printer, a convenient device for copying onto 
roll film selected frames from other rolls of film for subsequent 
automatic enlargement printing is described in detail. 














Tighe, N. J., Kreglo, J. R., Jr., Electron microscopy of peri- 
clase brick, J. Am. Ceram. Soc. 49, No. 2, 188-192 (Feb. 7, 1970). 


Key words: Ceramics; electron microscopy; MgO; microstruc- 
ture; sea-water periclase. 


Transmission electron microscopy was used to study microtructure 
of sea-water periclase brick. The material contained 96.7% MgO 
and had a porosity of 19%. Specimens were cut from the brick 
and thinned for microscopy by ionic bombardment. MgO was 
present as single grains (20 to 90 um) and as a fine precipitate 
(~100A) in an amorphous matrix. Dislocations in MgO grains 
were decorated with impurity precipitates. Second phase material 
was also present as single grains. Forsterite and monticellite were 
identified as second phases. 


Vidal, C. R., Cooper, Jr., Heat pipe oven: a new, well-defined 
metal vapor device for spectroscopic measurements, J. Appl. 


Phys. 40, No. 8, 3370-3374 (July 1969). 


Key words: Absorption spectroscopy; heat pipe; metal vapor 
device; oven. 


A new well defined metal vapor device called the heat pipe oven 
has been developed on the basis of the heat pipe, a heat conduc- 
tive element designed by G. Grover and his co-workers in Los 
Alamos. It generates continuously homogeneous vapors of well 
defined temperature, pressure and optical path length. The vapor 
is confined by inert gas boundaries which remove the window 
problem and allow a direct pressure measurement without relying 
on vapor pressure curves. Due to the continuous evaporation and 
condensation the vapor purifies itself during operation. 


Viezbicke, P. P., Edward, NBS Frequency and time broadcast 
services, Radio stations WWV, WWVH, WWVB, and 
WWVL, Nat. Bur. Stand. (U.S.), Spec. Publ. 236, 16 pages 
(1970) 25 cents, SD Catalog No. C13.11:236 (Supersedes NBS 
Special Publication 236, 1969 and previous editions). 


Key words: Broadcast of standard frequencies; high frequency; 
low frequency; standard frequencies; time signals; very low 
frequency. 


Detailed descriptions are given of the technical services provided 
by the National Bureau of Standards radio stations WWV, WWVH, 
WWVB, and WWVL. These services are: 1. Standard radio fre- 
quencies; 2. Standard audio frequencies; 3. Standard musical 
pitch; 4. Standard time intervals; 5. Time signals; 6. UT2 cor- 
rections; 7. Radio propagation forecasts; and 8. Geophysical 
alerts. In order to provide users with the best possible services, 
occasional changes in broadcasting schedules are required. This 
publication shows the schedules in effect on February 1, 1970. 
Annual revisions will be made. Current data relating to standard 
frequencies and time signals are available monthly in the Time 
and Frequency Services Bulletin. Advance notices of changes 
occurring between revisions will be sent to users of NBS broad- 
cast services who request such notice on the basis of need. 


Waclawski, B. J., Hughey, L. R., Adsorption of water vapor on 
polycrystalline tungsten, Surface Sci. 19, 464-468 (1970). 


Key words: Adsorption; hydrogen evolution; polycrystalline tung- 
sten; room temperature; ultra-high vacuum; water vapor; work 
function. 


Following an induction period, a spontaneous evolution of hydro- 
gen was observed to occur in the initial stages of adsorption at 
room temperature when polycrystalline tungsten was exposed to 
water vapor. Concurrent photoemission measurements indicated 
that a change in the work function of the sample occurred during 
the hydrogen evolution. Calculations of the rate of evolution, and 
the amount of hydrogen evolved are presented. 


Wall, L. A., Flynn, J. H., Straus, S., Rates of molecular vapor- 
ization of organic plasticizers, Polymer Eng. Sci. 10, No. 1, 
19-23 (Jan. 1970). 
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Key words: Energy of vaporization; internal energies of vaporiza- 
tion; molecular vaporization; organic plasticizers; plasticizers; 
solubility parameters. 


The kinetics of the molecular vaporization process of 21 plasti- 
cizers were investigated in detail. By both isothermal and noniso- 
thermal weight loss methods it was evident that 11 were quite pure 
single compounds while 10 were clearly mixtures of compounds. 
For the single component species internal energies, for vaporiza- 
tion and rates of volatilization are obtained from the kinetic data. 
The internal energies of vaporization are about ‘4 or less of values 
one can estimate from solubility parameters estimated by the 
additive factor method of Small. Thus solubility parameters based 
on our experimental energies of vaporization are low by about 
30%. From this and other work on linear n-alkanes it is con- 
cluded that in the molecular vaporization process, the large organic 
molecules studied evaporate approximately as spheres and hence 
low values for the energy of vaporization are obtained. Conse- 
quently the difference between our experimental energy and that 
estimated from solubility parameters is the energy for extending 
the molecule in a vacuum environment. 


Walton, W. W., Cadoff, B. C., Editors, Performance of build- 
ings—concept and measurement, Proceeding of the Ist con- 
ference in a series of conferences on Man and His Shelter held 
at the National Bureau of Standards, Gaithersburg, Maryland, 
September 23-25, 1968, Nat. Bur. Stand. (U.S.), Bldg. Sci. Ser. 1, 
132 pages (Jan. 1970), $2.75, SD Catalog No. C13.29/2:1. 


Key words: Building systems; performance of buildings; stand- 
ards; test methods; urban planning; user needs. 


The Conference entitled “Performance of Buildings—Concept and 
Measurement” was held at the National Bureau of Standards, 
Gaithersburg, Md. on September 23-25, 1968. This was the first 
in a planned series of conferences on “Man and His Shelter.” The 
purpose of these conferences is to bring together those people 
from various disciplines who may contribute to improving the 
quality of man’s shelter. At the present conference, papers were 
presented by nineteen authors representing government and in- 
dustry in such diverse disciplines as architecture, engineering, 
science, urban planning, and standards. These papers emphasize 
the prime importance of considering user needs in the develop- 
ment of performance criteria, the necessity of test methods to 
determine whether the desired performance has been achieved, and 
the development of performance specifications and standards. Ap- 
plication of these ideas to building systems, and to the planning 
and design of entire communities, is also discussed. 


Wampler, R. H., An evaluation of linear least squares com- 
puter programs: A summary report, (Proc. 14th Conf. Design 
of Experiments in Army Research Development and Testing, Edge- 
wood Arsenal, Md., October 23, 1968), ARO-—D Report 69-2, 103-— 
126 (U.S. Army Research Office-Durham, Durham, North Caro- 
lina, Sept. 1969). 


Key words: Computer programs; curve fitting; Gram-Schmidt 
orthogonalization; Householder transformations; iterative refine- 
ment; least squares; linear equations; orthogonal polynomials; 
orthogonalization; regression; rounding error; stepwise regression. 


Two linear least squares test problems, based on fifth degree 
polynomials, have been run on more than twenty different com- 
puter programs in order to assess their numerical accuracy. Among 
the programs tested were representatives from various statistical 
packages as well as some from the SHARE library. Essentially 
four different algorithms were used in the various programs to 
obtain the coefficients of the least squares fits. The tests were run 
on several different computers, in double precision as well as 
single precision. By comparing the coefficients reported, it was 
found that those programs using orthogonal Householder trans- 
formations or Gram-Schmidt orthonormalization were much more 
accurate than those using elimination algorithms. Programs using 
orthogonal polynomials (suitable only for polynomial fits) also 
proved to be superior to those using elimination algorithms. The 
most successful programs accumulated inner products in double 
precision and made use of iterative refinement procedures. In a 
number of programs, the coefficients reported in one test problem 








were sometimes completely erronous, containing not even one 
correct significant digit. 


White, H. J., Jr., Absorption of cationic surfactants by cellu- 
losic substrates, Chapter 9 in Cationic Surfactants 4, 311-340 
(Marcel Dekker, Inc., New York, N.Y. 1970). 


Key words: Absorption; adsorption; cationic surfactants; cellu- 
losic materials; review. 


The physical chemistry of the interaction of cationic surfactants 
with cellulosic substrates is reviewed. 


Wiederhorn, S. M., Crack propagation in polycrystalline 
ceramics, (Proc. 15th Sagamore Army Materials Research Conf., 
Sagamore Conference Center, Raquette Lake, New York, Aug. 
20-23, 1968), Chapter in Ultrafine-Grain Ceramics 15, 317-338 
(Syracuse University Press, Syracuse, N.Y., 1970). 


Key words: Ceramics; fracture; fracture surface energy; strength; 
thermal shock. 


In this paper, the strength and shock resistance of ceramic mate- 
rials will be related to microstructural details and the energy 
necessary to form fracture surfaces. The influence of grain size, 
grain orientation, crystal anisotropy and residual stresses on crack 
propagation will be discussed. It will be argued that the strength 
of ceramic materials is closely related to the grain boundary or 
single crystal fracture surface energy while shock resistance is 
related to the surface energy required to propagate large cracks 
through polycrystalline materials. The use of fracture surface 
energies as a research technique and as a design parameter will 
be discussed. 


Wilson, W. K., Development of specifications for archival 
record materials, Am. Archivist 33, No. 2, 219-223 (April 1970). 


Key words: Acidity; aging; alum sorption; laboratory aging; 
paper; records; stability; stability of ce’’viose; stability of paper; 
thermal analysis. 


Work in progress on the development of information on the sta- 
bility of paper and modified cellulose is described in relation to 
the objectives of the program. Plans for calendar year 1970 are 
outlined. 


Wilson, W. K., Reflections on the stability of paper, Restau- 
rator 1, No. 2, 79-86 (1969). 


Key words: Degradation; paper; preservation of records; records; 
stability. 


Variables in the composition of paper are discussed in relation to 
stability, and new approaches to research on preservation of rec- 
ords are described. 


Wood, R. E., Hamer, W. J., Theoretical electromotive forces 
of metal-halogen cells: Some recalculations based on recent 
data, J. Electrochem. Soc. 117, No. 1, 82-83 (Jan. 1970). 


Key words: Cells; emfs; free energies of formation; fused salts; 
halides; thermodynamics. 


Comparison of the theoretical emfs for cells containing a single 
solid or molten halide published by Hamer, Malmberg and Rubin 
in 1956 and 1965 with emfs calculated from tabulations in the 
JANAF tables published in 1965, 1966, 1967, and 1968 reveals 
some rather large discrepancies. This note tabulates the 34 cases 
in which such dicrepancies were found to exceed 100 millivolts, 
discrepancies equivalent to 2.3 kcal — one thermocalorie (de- 
fined) = 4.1840 joules. Equiv~* or more. 


Yakowitz, H., The divergent beam x-ray technique, Chapter 
in Advances in Electronics and Electron Physics, Suppl. 6, Elec- 
tron Probe Microanalysis, A. J. Toksimis and L. L. Marton, ed., 
pp. 361-431 (Academic Press Inc., New York, N.Y., 1969). 
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Key words: Kossel method; x-ray microdiffraction. 


This is a review paper discussing the technique of divergent beam 
(Kossel) x-ray diffraction. The Kossel method is studied with 
regard to line intensities, photographic contrast, data reduction, 
specimen preparation, and contribution to scientific research. 
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